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PREFACE. 


The  preparation  of  this  volume  was  undertaken  with  the 
view  of  presenting  in  as  simple  and  concise  a  manner  an 
possible  the  fundamental  principles  of  the  Method  of  Least 
Squares.  While  it  is  believed  that  everything  essential  to 
the  solution  of  all  ordinary  problems  has  been  included,  no 
attempt  has  been  made  to  develop  at  length  those  special 
methods  and  forms  that  are  so  useful  and  almost  necessary  in 
case  large  numbers  of  observations  of  certain  kinds,  such, 
for  instance,  as  those  met  with  in  geodetic  and  astronomical 
measurements,  are  to  be  adjusted. 

Frequent  references  throughout  the  text,  and  more  particu- 
larly the  list  01  works  given  on  page  v  of  the  Appendix,  will, 
however,  enable  the  student  to  extend  his  studies  in  what- 
ever special  direction  his  profession  may  require ;  it  being 
expected  that  this  book  will  in  such  cases  be  looked  upon 
merely  as  an  introductory  treatise.  All  of  the  works  men- 
tioned have  been  freely  consulted  in  the  preparation  of  these 
pages,  and  the  author  desires  in  particular  to  acknowledge 
his  indebtedness  for  many  of  the  examples. 

DANA  P.  BARTLETT. 
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THE  METHOD  OF  LEAST  SQUARES. 


CHAPTER   I. 

GENERAL   PRINCIPLES. 

1.  In  scientific  investigations  of  all  kinds  it  is  frequently 
necessary  to  determine  the  values  of  certain  quantities  by 
means  of  actual  measurements  either  with  or  without  the  aid 
of  instruments.  The  observations  may  be  made  directly  upon 
the  values  of  the  unknown  quantities  or  upon  certain  functions 
of  the  unknowns.  In  the  latter  case  the  values  of  the  required 
quantities  must  be  obtained  by  computation  from  the  observed 
values  of  the  functions.  In  order  to  obtain  more  accurate 
values  of  the  unknowns  than  would  be  given  by  a  single 
measurement,  or  set  of  measurements,  the  observations  are 
usually  repeated  either  in  the  same  way  and  under  the  same 
conditions  or  in  a  variety  of  different  ways  and  under  vary- 
ing conditions. 

Under  these  circumstances  it  will  invariably  be  found  that 
the  different  measurements  give  discordant  results,  the  amount 
of  the  discrepancies  varying  with  the  character  of  the  observa- 
tions ;  and  the  question  that  now  presents  itself  is  how  to 
determine  from  these  discordant  observations  the  true  values 
of  the  required  quantities.  From  the  nature  of  the  case, 
however,  we  can  not  expect  to  obtain  our  values  with  absolute 
accuracy;  all  that  we  can  hope  for  is  to  obtain  those  values 
which  are  rendered  most  probable  after  all  the  observations 
are  taken  into  account,  and,  further,  to  determine  the  degree 
of  confidence  that  can  be  placed  in  those  values. 
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2.  The  attainment  of   the  above   results  constitutes  the 
primary  object  of  the  Method  of  Least  Squares.     The  method 
is  also  employed  in  comparing  the  relative  worth  of  different 
measurements  of  the  same  quantity,  and  in  determining  the 
equation  of  a  curve  which  shall  suitably  represent  the  relation 
between  two  variables  in  cases  where  the  exact  law  connecting 
them  is  not  known. 

Also,  before  making  any  observations,  we  may  employ  the 
method  to  determine  how  precise  the  component  measurements 
of  a  series  must  be  in  order  to  yield  a  required  degree  of  pre- 
cision in  the  final  result;  or,  conversely,  to  determine  what 
the  precision  of  the  final  result  will  be,  knowing  the  precision 
attainable  in  the  component  measurements.  This  latter  appli- 
cation of  the  method  will  be  treated  at  length  in  the  course 
on  "  The  Precision  of  Measurements." 

3.  Errors.    The  cause  of  the  discrepancies  between  the 
results  of  our  different  observations  is  that  every  observation 
that  is  a  measure  is  subject  to  error.     These  errors  are  of  two 
kinds, — Constant  or  Systematic  Errors  and  Accidental  Errors. 

4.  Constant  Errors  are  errors  which  in  all  measures  of  the 
same  quantity,  made  with  the  same  care  and  under  the  same 
conditions,  have  the  same  magnitude,  or  whose  presence  and 
magnitude   are  due   to   some   fixed   cause.     These   constant 
errors  may  be   of   several  classes,  which  are  designated  as 
follows: — 

First.  Theoretical  Errors,  such  as  those  due  to  the  refrac- 
tion or  aberration  of  light,  the  effect  of  a  definite  change  in 
temperature  or  moisture  on  our  standards  of  measurement, 
etc.  As  soon  as  their  causes  are  known  the  magnitude  of 
errors  may  be  calculated  and  their  effect  eliminated 
from  the  observations. 

Second.  Instrumental  Errors,  such  as  errors  of  division  of 
graduated  scales,  defects  in  micrometer  screws,  eccentricity  of 
rirrh-s,  j-tc.  Thrsc  errors  will  he  dUrovcrnl  l.y  ;ui  r\;unin:i- 
tion  of  the  instruments  ami  tln-ir  rfl'rrts  rliminutril  t'rom  tin* 
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observations,   either  by  a   particular  method   of   using    the 
instruments  or  by  subsequent  computation. 

Third.  Personal  Errors.  These  are  due  to  personal  pecu- 
liarities of  an  observer,  who  always  answers  a  signal  too  soon 
or  too  late,  always  estimates  a  quantity  smaller  than  it  is,  etc. 
The  character  and  magnitude  of  these  errors  may  be  deter- 
mined by  a  study  of  the  observer,  his  "Personal  Equation" 
may  be  obtained,  and  his  observations  thus  corrected  for  thin 
source  of  error. 

5.  Mistakes.    Although  of  a  somewhat  different  character, 
these  should  be  considered  in  connection  with  constant  errors. 
A  mistake  is  made  when  a  figure  3  is  read  for  a  figure  8,  or 
when  in  reading  a  graduated  circle  which  is  numbered  in  both 
directions  the  angle  is  read  43°  instead  of  the  complementary 
angle  47°,  etc.     These  mistakes  are  usually  of  such  a  charac- 
ter that  they  may  be  detected  by  an  inspection  of  the  observa- 
tions and  a  proper  correction  made. 

6.  Accidental  Errors  are  errors  due  to  irregular  causes, 
whose  effect  upon  the  observations  is  not  determined  by  any 
circumstances  peculiar  to  that  particular  set  of  measurements, 
and  which  cannot  therefore  be   computed  and  allowed  for 
beforehand.     Such  errors  are  those  due  to  sudden  changes  in 
refraction  owing  to  sudden  and  unobserved  changes  in  tem- 
perature; unequal  expansion  of  different  parts  of  an  instru- 
ment with  change  in  temperature;  shaking  of  an  instrument 
in  the  wind,  etc.     But  most  important  of  all  are  those  errors 
which  arise  from  imperfections  in  the  sight,  hearing,  and  other 
senses  of  the  observer,  which  render  it  impossible  for  him  to 
adjust  and  use  his  instruments  with  absolute  accuracy. 

After  a  full  investigation  of  the  constant  errors,  the  observer 
should  diminish  the  accidental  errors  as  much  as  possible,  both 
in  number  and  magnitude,  by  taking  every  precaution  and  care 
in  the  measurements  themselves.  The  problem  now  remains 
to  combine  the  observations  so  that  the  remaining  accidental 
errors  shall  have  the  least  probable  effect  upon  the  results, 
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and  it  is  to  bring  about  this  combination  of  observations  that 
we  employ  the  Method  of  Least  Squares. 

When  no  more  observations  are  made  than  are  sufficient  to 
determine  one  value  for  each  of  the  unknown  quantities,  we 
must  accept  these  values  as  the  most  probable  ones.  But  if 
additional  observations  are  made  leading  to  discordant  results, 
we  can  not  take  any  one  of  them  as  the  correct  value,  and  in 
fact,  as  already  stated,  we  shall  probably  not  be  able  to  obtain 
the  true  values  of  the  unknowns.  All  that  we  can  do  is  to  find 
values  of  the  unknowns  which  shall  remove  the  discrepancies 
between  the  different  observations  and  which  shall  be  those 
values  that  are  rendered  most  probable  by  the  existence  of 
the  observations  themselves. 

On  first  thoughts  it  may  seem  that  these  accidental  errors, 
being  due  to  so  many  different  and  unknown  causes,  will  be 
beyond  the  scope  of  mathematical  investigation.  Neverthe- 
less, the  theory  of  probability  requires  that  these  errors  shall 
follow  in  magnitude  and  frequency  a  law  that  is  capable  of 
exact  mathematical  expression,  and  experience  confirms  the 
correctness  of  this  law. 

For  more  extended  remarks  on  these  subjects  see  — 

Holman,  "  Discussion  of  the  Precision  of  Measurements,"  pp.  1-14. 
Merriman,  "  Text-Book  of  Least  Squares,"  pp.  1-6. 
Chauvcnet,  "  Spherical  and  Practical  Astronomy,"  pp.  469-473. 
Wright,  "Treatise  on  the  Adjustment  of  Observations,"  pp.  11-18. 


LAWS  OF  ERRORS  OF  OBSERVATION. 

7.  The  derivation  of  the  general  laws  of  the  occurrence  of 
errors  of  observation,  and  of  the  processes  for  determining 
the  most  probable  values  of  the  unknown  quantities,  will  l>e 
based  upon  the  following 

Axiom.  If  a  series  of  n  direct  observations,  J/i,  M^  . . .  .Vn, 
U  made  upon  the  value  of  a  quantity  M,  all  the  observations 
b«-in«:  made  with  tin-  Kirne  care  and  under  the  same  eircmii- 
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stances,  the  most  probable  value  M0  of  that  quantity  is  the 
arithmetical  mean  of  the  observations.     Or 


MQ=-  -2T  .^»  =  ~«-.  (1) 

8.  The  Eeal  Error  (x)  of  an  observation  is  the  difference 
between  the  observed  value  of  the  measured  quantity  and  the 
real  value. 

9.  The  Residual  (v)  of  an  observation  is  the  difference 
between  the  observed  value  of   the  measured  quantity  and 
the  value  rendered  most  probable  by  the  existence  of  the 
observations. 

10.  Example.     Eight   observations  are    made    upon    the 
resistance  of   a  coil  of  wire,  the  true  resistance  being  512. 
Find  from  these  observations  the  most  probable  resistance, 
and  also  the  real  errors  and  residuals. 

Observations.          Real  Errors.  Residuals. 

x  v 

+  .4  +.30 

+  .2  +.10 

-  .1  -  .20 
+  .3  +  .20 

-  .2  -  .30 
+  .3  +  .20 

-  .1  -  .20 
.0                     -  .10 


From  the  observations,  then,  we  should  say  that  the  most 
probable  resistance  of  the  coil  is  512.10.  It  will  also  be 
noticed  that  the  sum  of  the  residuals  is  zero.  That  this  is  a 
general  result  following  from  the  assumption  of  the  arith- 
metical mean  as  the  most  probable  value  may  be  proved  as 
follows:  If  the  observations  are  M^  M^  .  .  .  Mn,  the 
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arithmetical  mean  M0,  and  the  residuals  vl9  u2,  .  .  .  vw,  then 
we  have 


.  .  .  v 


since  M0=  - 
n 

.-.    Sv  =  0  (2) 

11.  Weighted  Observations.    The  weight  of  an  observa- 
tion expresses  its  relative  worth  compared  with  other  obser- 
vations.    Thus,  if  six  observations  are  made  upon  the  value 
of  a  quantity,  five  of  which  give  the  same  result,  while  the 
sixth   differs,  in   combining  these   two   different   results    to 
obtain  the  most  probable  value    of  the   unknown,  the  first 
value  ought  to  have  five  times  the  influence  upon  the  final 
result  that  the  second  has,  since  it  has  taken  five  times  as 
much  labor  and  time  to  obtain  it.     Hence  in  general  we  may 
say  — 

12.  The  Weight  (p)  of  an  observation  may  be  considered 
as  representing  the  number  of  times  the  observation  has  been 
repeated  and  the  same  result  obtained. 

The  weights  assigned  to  observations  may  be  due  to  a 
variety  of  causes,  as  difference  in  skill  of  observers,  difference 
in  the  instruments  used  or  the  circumstances  under  which 
the  observations  are  made,  etc.  But  whatever  the  cause,  the 
effect  on  the  final  values  of  weighting  an  observation  will  be 
tin-  Karne  as  indicated  in  the  preceding  paragraph. 

13.  Example.     Suppose  n  observations,  J/j,  Mt,  .  .  .  Mn, 
of  weights   />!,  7>2,  .  .  .  />n,    are  made  upon  the  value  of  a 
quantity  M.     To  find  the  most  probable  value   M0   of  the 
quantity. 

From  the  above  interpretation  of  the  meaning  of  weight, 
we  may  consider  that  the  whole  number  of  observations  is 
/>i+/>2+  •  •  •  Pw  or  2j»,  and  that  the  result  3/j  has 
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obtained   in    pl   observations,    M2    in  p2   observations,  etc., 
Therefore,  by  (1) 

ffiJ 


3f0   is  called  the  General  Mean. 

If  the  residuals  are  v»  v2,  .  .  .  vn,  we  have 

Vi  =  Mv  —  MQ,     v2  =  Jf2  —  3f0,  .  .  .  vn  =  Mn  — 


=  0  (4) 

Which  shows  that  in  the  case  of  direct  observations  of  differ- 
ent weights  the  sum  of  the  weighted  residuals  is  zero. 

14.  If  the  observations  are   not  made  directly  upon  the 
values  of  the    required  quantities,  the  method  of  adjusting 
the  results  so  as  to   obtain  the   best  possible  values  of  the 
unknowns  will  depend  upon  the  laws  which  govern  the  dis- 
tribution of  the  errors  of  these  observations.     It  is  found  in 
practice  that  the  accidental  errors  of  observations  follow  cer- 
tain well  defined  laws,  and  what  these  are  may  best  be  seen 
by  taking  an  actual  example. 

15.  Example.     One  thousand  shots  are  fired  at  a  target 
which  is  divided  into  a  number  of  horizontal  sections  by  lines 
one  foot  apart,  the  centre   line  of  the  target  being  in   the 
middle  of  one  of  these  spaces.     The  shots  were  distributed 
as  follows :  — 

In  Space.          Shots.  In  Space.  Shots.  In  Space.  Shot*. 

to  +  44  1  4.14104-4  190  —  24  to  —  3J  79 

«  +34  4  +  i  "  --  i  212  -34  "  -44  16 

"  -j-24  10  -  4  "  -14  204  —44  "  —54  2 

"  +14  89  —14  «  -24  193 


Tn  this  case  the  errors  are  evidently  the  distances  of  the 
shots  from  the  centre  of  the  target.     Further,  as  far  as  can 
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be  judged  from  these  one  thousand  shots,  if  another  shot  is 
fired  the   probability  that  this  shot  will   fall   between    the 

lines 

is  .001 
«  .004 


-   i  and  -liis.204 


+  3* 


"  .010 
"  .089 

u 


+ 

-       "  .212 


-H 

-2* 
-3* 
-4* 


— 2£  «  .193 
-3$  "  .079 
— 4J  "  .016 
-5  "  .002 


The  sum  of  the  above  probabilities  is  unity,  and,  therefore, 
as  far  as  the  preceding  shots  show  the  1001st  shot  will  cer- 
tainly hit  the  target. 

16.  Now  using  as  abscissas  the  distances  of  the  horizontal 
lines  from  the  centre  of  the  target,  and  as  ordinates  the  num- 
ber of  shots  falling  in  the  corresponding  spaces,  we  may 
construct  the  following  figure  :  — 

T 


Figure  1. 


r  o 

And  if  the  entire  area  of  this  figure  is  taken  as  unity  thru 
the  area  of  each  rectangle  will  <lnn»te  the  probability  of  a 
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shot,  if  fired,  falling  within  the  corresponding  space  of  the 
target. 

The  graphical  representation  of  the  accidental  errors  of 
observation  will  always  give  a  figure  similar  to  the  above. 
Hence  denoting  errors  by  abscissas,  and  their  frequency  by 
ordinates,  the  law  of  error  of  any  series  of  observations  may 
be  represented  by  a  curve  whose  general  form  is  determined 
by  Figure  1.  This  curve  is  called  the  "  Curve  of  Error,"  and 
is  shown  in  Figure  2. 


I 


Figure  2. 


PDM 


In  order  that  this  curve  may  represent  exactly  the  distri- 
bution of  the  errors  in  any  given  series  of  observations  it 
ought  to  meet  the  axis  of  X  at  some  definite  distance  to  the 
right  and  left  of  the  origin  and  coincide  with  the  axis  from 
there  on,  for  in  all  actual  observations  there  is  a  limit  beyond 
which  no  errors  occur.  But  as  the  exact  point  of  meeting 
could  not  be  determined  for  any  given  case,  and  as  it  would 
not  be  possible  to  obtain  the  equation  of  such  a  curve,  we 
make  it  asymptotic  to  the  axis  of  A",  t-iking  care  that  the 
error  thus  introduced  shall  in  any  set  of  observations  be  so 
small  as  to  be  negligible. 
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17.  An  inspection  of  Figures  1  and  2  will  now  exhibit 
some  of  the  general  lawt  of  errors  of  observation  and  the 
corresponding  properties  of  the  curve  of  error. 


Laws  of  Error  derivfd  from  an 
inspection  of  Figure  1. 


Representation  of  these  laws  by  the 
Curve  of  Error. 


First.  Small  errors  are  more 
frequent  than  large  ones. 

Second.  Positive  and  negative 
errors  of  the  same  absolute  mag- 
nitude are  equally  likely  to  occur. 

Third  The  probability  of  the 
occurrence  of  very  large  errors  is 
very  small. 

Fourth.  The  frequency  of  any 
error  depends  upon  the  magnitude 
of  that  error. 


The  maximum  point  of  the 
curve  is  on  the  axis  of  Y. 

The  curve  is  symmetrical  with 
respect  to  the  axis  of  F. 

The  curve  is  asymptotic  to  the 
axis  of  X. 

The  equation  of  the  curve  will 
be  of  the  form 


y=*(«0  '       (6 


18.  If,  now,  the  total  area  between  the  curve  and  the  axis 
of  X  be  denoted  by  unity  the  probability  that  the  error  of 
any  given  observation  will  fall  between  the  magnitudes  x  and 
x-\-  dx  will  be  represented  by  the  area  included  between  the 
curve,  the  axis  of  J^  and  the  ordinates  of  the  curve  at  the 
errors  x  and  x  -f-  dx,  or  by 


y  dx=<f>(x)  dx 


(6) 


And  this  probability  will  be  known  as  soon  as  we  find  the 
form  of  the  function  <f>(x). 

19.  The  above  expressions  in  (5)  and  (6)  are  the  ones 
that  we  should  use  if  we  regard  the  curve  of  error  as  repre- 
M -n ting  the  law  of  occurrence  of  errors  of  observation.  If, 
however,  we  look  upon  the  curve  as  expressing  the  law  to 
which  we  mu>t  make  tin-  n-Milualt*  conform,  in  order  th:it  the 
values  of  the  unknown  <ju:iiititii'H  obtained  from  them  may  IK- 
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the  most  probable  values,  we  should  replace  x  by  v  and  use 
the  expressions 

y  =  *(»)  (7) 

and  ydv  =  <i>(v)dv  (8) 

for  (5)  and  (6),  respectively. 

THE  EQUATION  OF  THE  CURVE  OF  ERROR, 

20.  Let  n  observations,  all  of  the  same  weight,  with 
results  3/i,  Mz,  .  .  .  Mn,  be  made  upon  any  function  or 
functions  of  a  number  of  unknown  quantities  z1?  z2,  .  .  .  zq  ; 
and  let  the  residuals  of  M^  Mz,  .  .  .  Mn  be  w1?  v2,  .  .  .  un, 
and  the  probability  of  the  occurrence  of  these  residuals  be 
<£(tfi)  dv,  <j>(vz)  dv,  .  .  .  <£(yn)  dv,  respectively.  Then  the 
probability  of  the  simultaneous  occurrence  of  all  these 
residuals  will  be 


P=<t>(vl)  4>(v2)  .  .  .  +(vn)  (efo)»  (9) 

.  .  .  log<j>(vn)-}-?ilogdv     (10) 


Each  different  method  that  might  be  adopted  for  computing 
the  values  of  the  unknowns  z1?  z2,  .  .  .  zq  would  lead  to  a  dif- 
ferent set  of  residuals  Wj,  u2,  .  .  .  vn  ;  but  obviously  that  set 
of  values  of  z^  z2,  .  .  .  zq  should  be  considered  the  best  which 
corresponds  to  the  particular  set  of  residuals  v^  v2,  .  .  .  vnj 
the  probability  of  whose  occurrence  is  greater  than  that  of 
any  other  set. 

Therefore  the  most  probable  values  of  zx,  z2,  .  .  .  zq  are 
those  that  make  P  in  (9),  or  log  P  in  (10),  a  maximum. 

The  values  of  zl5  z2,  .  .  .  zq  corresponding  with  this  latter 
condition  are  those  that  satisfy  equations  (11).  It  maybe 
noticed  that  these  equations  also  express  the  preliminary 
conditions  leading  to  a  minimum  value  of  log  P,  but  the 
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nature  of  the  problem  is  such  that  a  maximum  value  of  P 
evidently  exists  while  a  minimum  does  not,  and  it  is  there- 
fore unnecessary  to  investigate  further  the  mathematical 
conditions  for  a  maximum.  Hence  we  have 


, 


Bzl 


i 

'<HV«  2 


, 


and  for  convenience  we  may  put 


substituting  in  (11)  we  have 

5y  5u  dv 


These  equations  contain  ;ill  the  unknowns  z^  z^  .  .  .  z^  and 

there  are  at*  many  r<|u:iti«niH  :is  unknowns,  lu-ncf  as  soon  as 
we  find  tin-  form  of  the  function      (/')    \ve  can  solve  these 
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equations  for  the  most  probable  values  of  zly  z2,  .  .  .  zq.  Since 
we  have  considered  the  general  case,  and  the  above  results 
are  to  hold  true  whatever  the  number  of  unknown  quantities 
and  the  form  of  the  functions  observed,  we  may  deduce  the 
form  of  i/f(u)  by  solving  a  special  example. 

Example.  Let  n  observations  of  equal  weight  be  made 
upon  the  value  of  a  single  unknown  z1?  with  results  M^  J/,, 
.  .  .  Mn,  and  let  the  residuals  be  v1?  w2,  .  .  .  vn.  Then  the 
most  probable  value  of  z:  is  given  by 


differentiating  with  respect  to  z1? 

auj  __          d«2  _  dvn 

-L      •  —  r\  .—       —  —       ^.  _—      •      •      •      "~™       rv  (    ol    ) 

dzl  dZi  dzl 

substituting  (a)  in  (14),  changing  all  the  signs,  we  have 

^(i»1)  +  *(t>a)+...*(0ll)=0  (b) 

But  in  this  case,  as  was  shown  in  (2), 

•  «i  +  v2  +  .  .  .  vn  =  0  (c) 

In  order  that  (b)  and  (c)  may  both  be  true  the  functional 
symbol  ^  must  indicate  multiplication  by  a  constant.  That 
is,  in  general 

{j,(v)  =  cv  (15) 

Substituting  this  in  (13)  and  (12), 


dv 

,.       =  <f>(v)cv  ^r- 
dz  ^v    '       dz 


therefore  .  --  =  -  =  cw 

<^>(y)       da 

Integrating,  log  <j>(v)  =  |c 
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Since  y=<f>(v)  is  the  equation  of  the  curve  of  error,  (7), 
we  may  therefore  write  it 


But  on  examination  of  the  curve,  y  is  seen  to  be  a  decreas- 
ing function  of  v,  and  hence  the  exponent  of  e  is  essentially 
negative.  Accordingly  we  will  write  our  equation  in  the 
form 

v  =  lee-*^  (16) 

the  values  of  k  and  h  depending  upon  the  character  of  the 
observations,  but  in  all  observations  of  the  same  kind  and 
weight  having  the  same  values. 

This  equation  represents  the  law  in  accordance  with  which 
the  residuals  must  be  distributed  in  order  that  the  best  results 
may  be  obtained  from  our  observations.  But,  as  before 
mentioned,  if  we  wish  our  curve  to  represent  the  most  prob- 
able distribution  of  the  real  errors  of  observation  we  should 
write  the  equation  in  the  form 

f     ' 

(17) 

Hereafter  we  shall  use  without  further  remark  either  form 
of  the  equation  according  to  the  aspect  in  which  we  are 
considering  our  curve. 

An  inspection  of  the  above  equation  will  show  that  it 
satisfies  all  the  conditions  noted  in  discussing  the  form  of  the 
curve  of  error  in  paragraph  17. 

21.  It  is  important  to  notice  that  in  all  discussions  in  tlu> 
Method  of  Least  Squares  the  number  of  observations  is  sup- 
posed to  be  large  and  always  greater  than  the  number  of 
unknown  quantities.  As  will  be  illustrated  later  on,  para- 
graph 'H,  whenever  this  is  the  case  there  is  a  remarkable 
agreement  bet  ue<n  the  results  obtained  in  practice  and  those 
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indicated  by  the  theory.  And  even  when  the  observations 
are  few  in  number  the  method  still  affords  the  best  means  at 
our  command  for  their  adjustment,  the  results  obtained 
merely  having  a  smaller  weight  than  they  would  have  had  if 
derived  from  a  greater  number  of  observations. 

THE   METHOD  OF  LEAST  SQUARES. 

22.  We  are  now  in  a  position  to  see  whence  comes  the 
name  "  Least  Squares." 

In  paragraph  20  it  was  pointed  out  that  whenever  we 
make  a  series  of  observations,  each  observation  of  the  set 
having  the  same  weight,  the  most  probable  system  of  values 
of  the  unknown  quantities  will  be  that  which  corresponds 
with  the  set  of  residuals  the  probability  of  whose  occurrence 
is  a  maximum.  That  is,  the  best  set  of  values  of  the 
unknowns  will  be  that  which  gives  a  maximum  value  to 

P  =  <t>(vl)  $(v,}  .  ..*(y»)   (dv)n 
But  from  equation  (16)  this  reduces  to 

P  =  kne  -*2  W  +  v**  +  •  -  •  *»2)  (  dv)n  (18) 

Since  the  exponent  of  e  in  this  expression  is  negative, 
evidently  P  will  be  a  maximum  when 


v?  -\-vf-\-  •  •  •  v»2  =  Sva  is  a  minimum.     (19) 

Hence  the  adjustment  of  observations  by  the  Method  of 
Least  Squares  is  based  upon  the  principle  that  the  most  prob- 
able system  of  values  of  the  unknowns  is  that  which  renders 
the  sum  of  the  squares  of  the  residuals  a  minimum.  Hence 
the  name. 

The  conditions  for  a  maximum  value  of  P  were  expressed 
in  equations  (14),  and  since  it  has  been  shown  that  the 


16  METHOD   OF  LEAST  SQUARES. 

function  \f/  means  multiplication  by  a  constant,  those  equa- 
tions reduce  to  the  following,  called 

NORMAL    EQUATIONS. 


#Vf      ,  3lJn_ 

An  inspection  of  these  equations  will  show  that  they  also 
express  the  conditions  that  will  make  the  sum  of  the  squares 
of  the  residuals  a  minimum. 

In  the  adjustment  of  observations  the  above  are  the  funda- 
mental equations.  In  order  to  obtain  the  most  probable 
values  of  the  unknowns  in  any  set  of  observations,  all  that  is 
necessary  is  to  form  the  Normal  Equations  for  that  set  and 
solve  them  simultaneously.  The  examples  already  solved  for 
direct  observations  are  merely  special  cases  of  the  above 
general  solution. 


CHAPTER  II. 

THE   ADJUSTMENT   OF   OBSERVATIONS. 
INDIRECT   OBSERVATIONS. 

23.  In  the  determination  of  the  values  of  quantities  by 
means  of  observations  the  functions  of  the  unknowns  that  it 
is  necessary  to  observe  may  be  of  any  form,  but  if  they  are 
not  linear  the  normal  equations  derived  from  them  are  likely 
to  be  complicated  and  difficult,  if  not  impossible,  to  solve. 
Hence  if  the  observations  are  not  upon  linear  functions  of 
the    unknowns,    the   first   step  will   be   to   reduce   them   to 
equivalent  linear  expressions  by  transformations  depending 
on  the  character  of  the  observed  functions;  see  paragraphs 
43  and  44.     It  will  be  necessary  to  consider,  therefore,  the 
method  of  adjusting  observations  on  linear  functions  alone, 
and  the  procedure  in  cases  of  this  kind  may  be  illustrated  by 
the  following  simple  example. 

24.  Example.     8^  /#2,  $s,  are  three  solids  whose  masses 
are  required.     Not  having  standard  weights  enough  to  obtain 
all  these  masses  directly,  by  varying  the  distribution  of  the 
solids  in  the  pans  of  the  balance  the  following  observations 
are  made  :  — 

AS'1=^2-f-  1.7  grams. 
£  =  2.4      « 


tf,  =  #3  --  3.0       « 

If,  now,  the  most  probable  values  of  S^  S^  £8,  are  repre- 
sented by  Zj,  z2,  2s>  and  the  corresponding  residuals  of  the 
observations  by  v^  v2»  •  •  •  vn  the  values  of  these  residuals  in 
terms  of  z1}  z2,  z8  may  be  found  from  the  above  observations 
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by  transposing  all  terms  to  the  first  members  of  the  equa- 
tions, and  we  obtain  at  once  the  following,  called 

OBSERVATION    EQUATIONS. 

21  —  22  —   1.7  =  Vl 

za  —  2.4  =  u2  A. 

_21  +  22  +  23-  1.0  =  ^ 
z2  —  23  —  3.0  =  v4 

Applying  equations  (20)  as  formulas  we  have  the 

NORMAL    EQUATIONS. 
(2l  -22-  1.7)  +  (_Zl  +  22  +  z3-1.0 


+  (z2-2s-3.0)  =  0     (B) 
(2,-  2.4)  +  (-zl  +  z,  +  z3  -1.0) 

-f(22-23-3.0)(_l)  =  0 

Simplified,  these  become 

2z1_2z2-    23  —  0.7=0  (a) 

—2.3  =  0  (b) 

+323  —  0.4  =  0  (c) 


»X(») 

(«=) 

2X(b) 
substitute 

M 

QZl  —  622  — 

3z8  —  2.1  =  0 
3z3-  0.4=0 

(d) 
(•) 

(0 

52!  —  62, 

-42,4-62, 

—  2.5=0 
-4.6=0 

«1 

(d)  in  (b) 
(d)  in  (c) 

—  7.1  =  0 
z,=  7.1 
z,  =  5.5 
2,  =  2.5 

An  inspection  of  the  work  in  this  example  will  show  that. 
for  the  adjustment  of  ol>MTvati<>n->  of  c<juul  weight  on  liiu-ar 
finirtionH  of  the  unknowns  \\  »•  may  «l«'ri\c  the  following:  — 
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25.  Rule.     For  each  observation  write  an  "  Observation 

"  /  then  for  each  imkiioirn  form  a  "Normal 
"  by  inultijtlyiity  the  first  member  of  each  observa- 
tion equation  by  the  coefficient  of  that  unknown  in  tlmt 
equation,  adding  the  results  and  placing  the  sum  equal  to 
zero.  Solve  these  equations  simultaneously  for  the  values 
of  the  unknowns. 

In  solving  for  the  most  probable  values  of  the  unknowns 
the  second  members  of  the  observation  equations  are  very 
commonly  written  zero  instead  of  vt,  «2,  .  .  .  vn.  For  this  is 
the  form  in  which  the  equations  naturally  appear,  and  if  the 
observations  were  exact  the  residuals  would  actually  all  be 
zero.  The  method  of  solution  is  the  same  in  either  case. 

26.  Observations  of  Unequal  Weight.    If  the  observations 
are    not   all  of  equal  weight   the  same   method  will  apply, 
except  that  in  the  formation  of  the  normal  equations  each 
observation  equation  will  be  used  the  number  of  times  denoted 
by  its  weight.     Thus  in  the  last  example  if  the  observations 
have  the  weights  4,  9,  1,  4,  the  normal  equations  will  have 
the  same  form  as  in    (B),  page  18,  but  each  part  of  each 
equation  will  be  multiplied  by  the  weight  of  the  observation 
equation  from  which  it  is  derived.     This  will  give  the 

NORMAL    EQUATIONS. 


—  28  —  3.0)  =  0 

9(23  -  2.4)  +  (-  2L  +  22  +  23-  1.0) 

_L4(22-28-3.0)(-l)  =  0 
or  reduced 

52!  -5z,  -       23—    5.8—0  (a) 

—  5zl  +  922—    323—    6/2  =  0  (b) 

—  2l  —  322+1423—  10.6=0  (c) 

the  solution  of  which     i 


z,  =  7.07       22  =  5.42       2;1  =  ±4-2 


20  METHOD    OF  LEAST  SQUARES. 

Further  it  will  at  once  be  seen  that  if  pl9  pz,  .  .  .  pn,  are  the 
weights  of  the  corresponding  observations,  equations  (20) 
take  the  general  form:— 

WEIGHTED   NORMAL    EQUATIONS. 

BVz    .  Bvn 


(21) 


Hence  for  the  formation  of  the  normal  equations  in  weighted 
observations  on  linear  functions  of  the  unknowns,  we  have 
the  following:  — 

27.  Rule.     For  each  observation  write  an  "  Observation 
Equation "  /     then  for   each   unknown  form,   a   "Normal 
Equation"  by  multiplying  the  first  member  of  each  observa- 
tion equation  by  the  coefficient  of  that  unknown  in  that 
equation  and  by  the  weight  of  that  equation,  adding  the 
results  and  placing  the  sum  equal  to   zero.     Solve   these 
equations  simultaneously. 

28.  The  same  result  will  be   obtained  if   we   begin   by 
multiplying  each  observation  equation  by  the  square  root  of 
its   weight   and   then    proceed   according   to   the   first   rule 
(paragraph  25). 

This  result  illustrates  the  important  principle  that  multi- 
plying a  set  of  equations  by  the  square  roots  of  their  weights 
reduces  them  all  to  equivalent  equations  of  weight  unity. 

29.  Relation  between  the  Weight  of  an  Observation  and 
the  Value  of  h.     If  in  paragraph  '20   the    n    observations 
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have   weights  />1?  />2,  .  .  .  pn,    an(l   tne   quantity    h   values 
hi,  A2»  •  «  «  ^w»  tnen  equation  (18)  becomes 

^  (dv)n 
•  •  •  h'l*Vn^  (dv)n     (22) 

Tlie  most  probable  set  of  values  of  the  unknowns  is  that 
which  makes  P  a  maximum,  and  P  is  a  maximum  when 

7iiHV2+  hfvf-\-  .  .  .  hnzvn2  is  a  minimum.    (23) 

The  conditions  for  a  minimum  value  of  this  expression  are 
the  following,  which  are  then  for  this  case  the 

NORMAL    EQUATIONS. 


But  equations  (21)  are  also  the  normal  equations  for 
this  case.  Hence  (21)  and  (24)  must  be  identical,  and 

Pi  :  P*  :  ...  Pn  =  /*i2  :  7i22  :  .  .  .  /*,/-        (25) 

That  is,  the  square  of  A  is  proportional  to  the  weight  of 
the  observation.  Accordingly,  since  h  increases  in  value  as 
the  quality  of  the  observations  is  improved,  it  is  called  "The 
Measure  of  Precision." 
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Further,  it  follows  from  (23)  and  (25)  that  the  most 
probable  system  of  values  of  the  unknowns  will  be  that  in 
which 

Pi  ^\-\-PI  v22+  •  •  •  Pn  Vn2  is  a  minimum.     (26) 

And  this  is  the  most  general  form  of  statement  of  the 
principle  of  "  Least  Squares."  The  same  principle  is  repre- 
sented in  equations  (21). 

30.  Computation  of  Corrections.    If  large  numbers  occur 
in  the  observations  it  is  better  to  compute  the  most  probable 
corrections  to  apply  to  the  observed  values  rather  than  the 
most  probable  values  of  the  unknowns  themselves.     In  this 
way  we  can  often  avoid  a  large  amount  of  numerical  work. 

31.  Example.     Pl9  P^  P8,  P4,  P6,   are  five  points  whose 
altitudes  above  the  mean  level  of  the  sea  are  to  be  determined 
from  the  following  observations  of  difference  of  level. 

Pl  =  573.08  P4-  Pz=  170.28 

PZ-P1  =  2.60  P4  —  P!>  =  425.00 

P2  =  575.27  PB  =  319.91 

P8—P2=  167.33  PI>  =  319.75 

P4-P9=  3.80 

An  inspection  of  these  observations  shows  that  we  may  put 

Pl  =  573  +  «j  P4  =  745  +  z4 

/>f  =676  +  1,  P6  =  320  +  z6  (A) 


where    Zj,  z«,  zs,  ZH  zs»   are   small   corrections   whose   most 
probable  values  are  to  be  determined.     We  now  have  for 

OBSERVATION    EQUATIONS. 

573  +  Zj  -  573.08=0        or  z,-.08  =  0 

-573  —  *!-      2.60=0       or       z2  -  zl  -.(>()  =  0 
575  _|_  Za  _  575.27  =  0       or  «,—  .27  =  0 
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742  +  z3- 575-  z2-  167.33  =  0  or  za-z2- .33  =  0 

745  _j_  Z4  —  742  _  z3  -       3.80=0  or  24-z8— .80=0 

745  _|_  z4  _  575  _  za  _  170.28  =0  or  z4  -  z2  -  .28  =  0 

745  _L.  z4  _  320  -zc-  425.00=0  or  z4  —  z6            =0 

320  +  z6-  319.91  =  0  or  z6  +  .09=0 

320  +  zc-  319.75=0  or  z6-j-.25  =  0 

From  these  we  now  form  the 

NORMAL    EQUATIONS. 

+    .52=0 
z4  .26=0 

z4  -f    .47  =  0 

z2-     z3-f3z4-    z5  —  1.08  =  0 
•    z4  +  3z5-f    .34  =  0 
and  solving, 

2l=-.19;    z2  =  .14;    z3  =  .05;    z4  =  .43;    z5  =  .03 

Substituting  these  in  equations  (A)  we  have  for  the  most 
probable  altitudes, 

^  =  572.81          P3  =  742.05          P6  =  320.03 
P2  =  575.14          P4=  745.43 

If  the  original  observation  equations  had  been  retained, 
the  independent  terms  in  the  normal  equations  would  have 
been 

570.48         240.26         163.53         599.08         214.66 

32.  Significant  Figures.  The  adjustment  by  the  Method 
of  Least  Squares  of  observations  which  occur  in  practice, 
although  not  difficult,  is  apt  to  be  long  and  laborious.  Hence 
to  reduce  this  labor  as  much  as  possible  it  is  of  great  import- 
ance that  careful  attention  should  be  given  in  the  solutions 
to  the  proper  use  of  significant  figures.  When  in  doubt, 
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however,  as  to  the  proper  number  of  figures  to  retain  it  is 
better  to  keep  too  many  rather  than  too  few,  as  the  superfluous 
figures  can  be  rejected  at  the  end  of  the  computation ;  while 
if  too  few  are  retained  the  results  obtained  from  the  compu- 
tations will  be  worthless. 

For  a  general  discussion  of  the  subject  of  significant  figures 
see  Holman's  "  Precision  of  Measurements,"  pages  76  to  84, 
but  for  the  present  the  following  rules  will  suffice  for  most 
cases. 

Rule  1.  In  casting  off  places  of  figures  increase  by  1  the 
last  figure  retained,  when  the  following  figure  is  5  or  over. 

Rule  2.  In  the  precision  measure  retain  two  significant 
figures. 

Rule  3.  In  any  quantity  retain  enough  significant  figures 
to  include  the  place  in  which  the  second  significant  figure 
of  its  precision  measure  occurs. 

Rule  4.  When  several  quantities  are  to  be  added  or 
subtracted,  apply  Rule  3  to  the  least  precise  and  keep  only 
the  corresponding  figures  in  the  other  quantities. 

Rule  5.  When  several  quantities  are  to  be  multiplied  or 
divided  into  each  other,  find  the  percentage  precision  of  the 
least  precise.  If  this  is 

1  per  cent  or  more,  use  four  significant  figures. 
.1     "        "      "       "         "    five  "  " 

jO I  u       «      ..       u         u    8{x  u  « 

in  all  the  work.  If  the  final  result  obtained  in  this  way 
conflicts  with  Rule  3,  apply  the  latter. 

Rule  6.  When  logarithms  are  used,  retain  as  many 
places  in  the  mantissce  as  there  are  significant  figures 
retained  in  the  data  under  Rule  6. 

The  application  of  these  rules  is  not  always  possible  in  the 
course  of  the  work,  since  the  precision  measures  may  not  be 
known  until  the  end  of  the  computation.  But  as  ;i  <_'*  ;u-r:il 
it  IH  aufficient  in  direct  observations  to  retain  <>nc  more 
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place  of  figures  than  is  given  by  the  individual  observations, 
and  in  indirect  observations  to  retain  two  additional  places. 

CONDITIONED  OBSERVATIONS. 

33.  Conditioned  Observations  are  those  in  which  the 
unknown  quantities  must  be  determined  not  only  so  as  to 
satisfy  as  closely  as  possible  the  observation  equations,  but 
also  so  as  to  satisfy  exactly  certain  other  conditions.  These 
conditions  must  be  less  in  number  than  the  unknown  quantities, 
otherwise  the  unknowns  could  be  determined  from  the  con- 
ditions alone. 

The  adjustment  of  observations  of  this  class  may  be  reduced 
to  the  method  already  used  for  unconditioned  observations  in 
the  following  way. 

The  observations  are  represented  by  "  Observation  Equa- 
tions," and  the  conditions  by  certain  other  equations,  called 
"  Condition  Equations." 

Between  these  two  sets  of  equations  we  will  eliminate  as 
many  unknowns  as  there  are  conditions.  From  the  resulting 
equations,  which  will  be  the  same  in  number  as  the  observa- 
tions, we  will  form  in  the  usual  manner  the  "Normal 
Equations"  for  the  remaining  unknowns.  Having  solved 
these  normal  equations  and  substituted  the  results  in  the  con- 
dition equations,  we  shall  obtain  the  values  of  the  unknowns 
first  eliminated. 

All  conditions  of  the  problem  are  now  fulfilled,  for  the 
condition  equations  are  satisfied  exactly  and,  moreover, 
according  to  the  principle  of  Least  Squares  our  results 
are  those  rendered  most  probable  by  the  existence  of  the 
observations. 

As  in  the  example  last  considered,  it  is  often  more 
advantageous  to  compute  corrections  to  the  observed  values 
of  the  unknown  quantities  rather  than  the  values  of  the 
quantities  themselves. 
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34.    Example.    Find  the  most  probable  values  of  the  angles 
of  a  quadrilateral  from  the  observations, 

^4  =  101°  13'  22"  weight  3 

B=   93  49    17        «        2  A, 

C=    87  5   39        "        2 

D=    11  52   40        "        1 


0'  58" 
The  condition  to  be  satisfied  is  in  this  problem 

(B) 


Let   gj,  z2>  zs»  24   ^e  the  most  probable  corrections  to  add 
to  the  observed  values.     This  gives  for 

OBSERVATION   EQUATIONS 

zl  =:  0  weight  3 

*2=0        «  2 

z3^0        «  2 

z4=0        »  1 
and  for  the 

CONDITION    EQUATION 
Zi  +  z2+Z8  +  *4  +  58=0  (D) 

Eliminating  z4  between  (D)  and  (C),  the  equations  from 
which  the  normal  equations  are  to  be  derived  become 


,  =  0  weight  H 

t=0        «      2 

z,  =  0        "       2 


zt=0        «      2  E 
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Applying  the  rule  in  paragraph  27,  these  give  the 

NORMAL    EQUATIONS 

4Zl-      22         23       58:=0 

0  (F) 


Solving,  and  substituting  the  results  in  equation  (D),  we 
find 

z,=  -    8.29  28=-  12.43  G, 

22=-  12.43  z<=-  24.85 

Applying  these  corrections  to  the  observations  (A),  the 
most  probable  values  of  the  angles  are 

A  =  101°  13'  13".71 

B=    93  49         4.57                            HX 

C=    87  5  26.57 

D=    77  52  15.15 

Note.  In  eliminating  unknowns  between  the  observation 
and  condition  equations  care  must  be  taken  that  the  obser- 
vation equations  are  not  combined  with  each  other  or  multi- 
plied by  any  quantity.  For  if  this  is  done  the  weights  of 
the  observation  equations  will  be  altered.  (See  §28.) 

35.  In  the  above  example  it  is  evident  that  the  corrections 
to  be  applied  to  the  different  observations  are  inversely  as 
their  weights.  And,  in  general,  when  there  is  but  one 
equation  of  condition,  the  observations  expressing  direct 
determinations  of  the  unknowns,  the  corrections  will  be  pro- 
portional to  the  coefficients  of  the  unknowns  in  the  equation 
of  condition  divided  by  the  weights  of  the  corresponding 
observations.  A  proof  of  this  is  given  in  paragraph  117. 
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The  most  common  case  is  that  in  which  these  coefficients 
are  all  unity,  as  in  the  example  just  solved,  and  we  may  then 
derive  the 

Rule.  Find  the  difference  between  the  theoretical  and 
observed  results  and  divide  this  correction  among  the 
observations  in  the  inverse  ratio  of  their  weights. 

In  the  last  example  the  sum  of  the  observed  angles  exceeds 
360°  by  58".  Therefore  the  correction  to  be  applied  to  A  is 


~  58  X  -    ^~  -  =  -  58  X       =  -  8.29 


EMPIRICAL   FORMULAS  AND  CONSTANTS. 

36.  In  the  work  so  far  considered  the  observations  are 
supposed  to  be  made  either  directly  upon  the  values  of  the 
unknown  quantities  or  upon  some  function  of  the  unknowns 
whose  form,  and  the  constants  entering  into  it,  are  definitely 
known.  But  another  sort  of  problem  frequently  occurs, 
in  which  observations  are  made  upon  the  values  of  a  certain 
variable  and  the  corresponding  values  of  some  function  of  it, 
the  exact  form  of  the  function  not  being  known.  The  object 
in  this  case  is  the  determination  of  the  most  probable  form  of 
the  function  and  the  values  of  the  constants  involved  ;  that 
is,  the  derivation  of  the  algebraic  expression  best  representing 
the  law  connecting  the  variable  and  function. 

This  expression  may  be  looked  upon  as  the  equation  of  a 
curve,  abscissas  denoting  values  of  the  variable  and  ordinates 
values  of  the  function,  and  for  all  values  of  the  variable 
within  the  range  of  the  observations  we  may  determine  from 
it  the  most  probable  values  of  the  function  corresponding. 
But  except  in  special  cases,  where  tin-  number  of  observations 
i>  large,  where  the  law  connecting  variable  and  function  is 
well  detim-.l,  and  where  the  equation  obtained  is  an  aeeurate 
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representation  of  this  law,  it  cannot  be  assumed  to  apply 
beyond  the  range  of  the  observations.  And  in  no  case  would 
it  be  safe  to  make  use  of  the  curve  very  far  beyond  the  limits 
of  the  observations. 

37.  The  Method  of  Least  Squares  will  not  assist  in  deter- 
mining the  form  of  the  function.     This  must  be  settled  upon 
beforehand,   either   from   theoretical    considerations    or    by 
constructing  a  plot,  using  values  of  the  variable  as  abscissas 
and  of   the  function  as  ordinates,  when  the  smooth  curve 
drawn  through  the  points  thus  obtained  will  indicate  the 
form  of  equation  to  be  used. 

It  is  to  be  observed  that  this  is  a  method  of  trial  and  will 
not  necessarily  give  the  most  probable  form  of  the  function  ; 
and  in  fact  we  may  not  be  able  to  obtain  the  form  that  would 
be  absolutely  best.  Further,  several  forms  of  equation  may 
be  known  which  would  represent  well  the  plotted  points.  In 
such  a  case  that  should  be  considered  the  best  in  which  the 
sum  of  the  squares  of  the  residuals  is  found  to  be  the  least. 

38.  As  soon  as  the  form  of  the  function  is  decided  upon  it 
should  be  reduced  to  the  linear  form,  and  the  determination 
of   the  values  of   the  constants  involved  is  then  a  simple 
application  of  the  preceding  methods. 

As  the  "  Observation  Equations  "  in  any  given  problem  will 
all  be  of  the  same  kind,  it  is  usually  advisable  to  write  out 
the  general  form  of  the  "Normal  Equations"  and  arrange 
the  computations  in  tabular  form,  while  the  retention  of 
tlu'  proper  number  of  significant  figures  is  of  particular 
importance  in  this  work. 

39.  A  case  that  frequently  occurs  is  that  in  which  the 
quantity    y   is  a  constantly  increasing  function  of  the  vari- 
able &,  or  where  the  plotted  curve  is  approximately  parabolic 
in  form.     Here  the  equation 


...  (27) 

may  be  taken  to  represent  the  relation  between  the  variable 
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and  the  function.  The  larger  the  number  of  terms  taken  in 
the  second  member,  the  more  accurately  may  the  equation 
obtained  be  made  to  represent  the  results  of  the  observations  ; 
but  the  labor  involved  increases  rapidly  with  increase  in  the 
number  of  terms,  and  if  the  plot  shows  a  very  nearly  straight 
line  the  first  two  terms  alone  may  suffice. 

40.    Example.     In  measuring  the  velocity  of  the  current 
of  a  river  the  following  results  were  obtained : — 

Depths.  Velocities. 

x  V 

1  4.86 

2  5.14 

3  5.15 

4  4.85 

5  4.24 

6  3.36 

7  2.16 

8  0.67 

The  velocity  at  the  surface  is  4.250.     Find  the  equation  of 
a  curve  which  will  express  the  relation  between  x  and  V. 
Plotting  the  observations  we  find  the  curve 

V 

•f 
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This  is  approximately  parabolic  in  form  and  passes  through 
the  fixed  point  (0,  4/25).  Therefore  the  relation  between  x 
and  V  may  be  expressed  by  the  equation 


V=  4.2 


6V 


(A) 


and  substituting  in  this  the  corresponding  values  of  x  and  V 
as  given  by  tlu-  observations,  we  shall  have  eight  observation 
equations  from  which  the  most  probable  values  of  B  and  C 
are  to  be  computed. 

All  of  the  observation  equations  being  of  the  form  (A)  we 
have  the 

NORMAL    EQUATIONS 


B 


.25 
.25 


-  2  Vx2  —  0 

-  2  Vx  =  0 


(a) 

(b) 


For  computing  the  coefficients  in  these  equations  it  is  most 
convenient  to  arrange  the  following  table. 


X 

V 

Vx 

x2 

Fe2 

x9 

x4 

1 

4.86 

4.86 

1 

4.86 

I 

1 

2 

5.14 

10.28 

4 

20.56 

8 

16 

3 

5.15 

15.45 

9 

46.35 

27 

81 

4 

4.85 

19.40 

16 

77.60 

64 

256 

5 

4.24 

21.20 

25 

106.00 

125 

625 

6 

3.36 

20.16 

36 

120.96 

216 

1296 

7 

2.16 

15.12 

49 

105.84 

343 

•2401 

8 

0.67 

5.36 

64 

42.88 

512 

4096 

36 

111.83 

204 

525.05 

1296 

8772 

lx 

SVx 

2*2 

2  Far1 

2*' 

2*4 
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Substituting  these  results  in  (a)  and  (b)  we  have 

8772  C+  1296  B-\-  341.95  =  0  (c) 

1296(7+    204  ./?-{-   41.17  =  0  (d) 
and  solving, 

C  =-.1493             J?  =  .7465  (e) 

Therefore  the  required  equation  is 

V  =  4.25  +  .7465x  -  .1493*'  (f  ) 

Whenever  the  quantities  in  the  observations  are  so  large 
that  the  use  of  logarithms  is  desirable,  these  can  best  be  put 
in  the  same  columns  directly  over  the  natural  numbers. 

41.  It  is  to  be  remarked  that  in  solutions  like  the  above 
we  assume,  from  our  method  of  forming  the  observation  and 
normal  equations,  that  the  observations  on  the  values  of  the 
function  are  alone  subject  to  error,  the  observations  on  the 
variable  being  supposed  to  be  exact  or  to  have  errors  so  small 
as  to  be  negligible. 

42.  Periodic  Phenomena.    If  as  the  variable  increases  the 
function  passes  through  recurring  values,  that  is,  if   y   is   a 
periodic  function  of   a?,  some  form  of  trigonometric  equation 
would  be  the  proper  one  to  select.     For  instance,  a  good  form 
to  use  is 


y  =  A  +  B  sin  C  cos^fx         (28) 

where  -4,  jB,  C  are  the  constants  whose  most  probable  values 
are  to  be  found,  and  ra  is  the  number  of  units  of  x  comprised 
in  the  entire  cycle  of  values  of  y.  The  quantity  m  is  to  be 
determined  from  an  inspection  of  the  observations,  and  if  it 
appears  that  several  values  of  m  might  be  used,  all  should  be 
tri.  <1,  and  that  which  leads  to  the  smallest  sum  for  the  squares 
of  the  residuals  is  to  be  considered  the  best.  If  the  several 
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cycles  are  not  similar  and  regular,  additional  terms  involving 
multiples  of  x  will  have  to  be  added  to  equation  (28). 

43.  Special  Treatment  of  Exponential  Equations.    If  the 
equation  selected  is  not  in  the  linear  form  as  regards  the 
unknown    constants    the   general    method    of    procedure   is 
given  in  paragraph  44,  but   in  some  special  cases   a  more 
simple  reduction  is  possible.     A  quite  common  case  is  the 
following: — 

Suppose  the  relation  between  x  and  y  is  expressed  by  the 
equation 

y  =  kxm  (29) 

the  problem  being  to  obtain  the  best  values  of  k  and  m. 
Taking  the  logarithms  of  both  members  of  (29), 

log  y  =  log  k  -\-  m  log  x 

Denoting  log  k  by  k',  this  equation  becomes 
m  log  x  -j-  k'  —  log  y  =  0 

which  is  in  the  linear  form  as  regards  the  unknowns  m  and 
k',  and  the  normal  equations  may  now  be  formed  in  the  usual 
manner. 

The  most  convenient  way  of  determining  whether  equation 
(29)  is  a  suitable  one  to  select  or  not  is  to  plot  the  corre- 
sponding values  of  x  and  y  on  logarithmic  cross-section 
paper.  If  (29)  is  a  proper  equation  to  use  the  plotted  points 
will  lie  upon  a  straight  line. 

REDUCTION  OF  OBSERVATION  EQUATIONS  TO  THE 
LINEAR  FORM. 

44.  When   the  observation  equations  are    not   linear   as 
regards  the  unknowns,  the  only  practicable   method  of  pro- 
cedure, as  already  mentioned  in  paragraph  23,  is  to  reduce 
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them   to   that   form.     This   reduction    may  be    effected    as 
follows  :  — 

Let  the  observation  equations  be 

,/UZ,,  Z2,  .  .  .  Zq)  =  M, 
t,  Z2,  .  .  .  Zq)  =  J/8 

.....  (A) 


in  which  Zj,  Z2,  .  .  .  Zg  represent  the  unknown  quantities  and 
J/u  J/2,  .  .  .  Mn  the  observations,  the  functions  being  of 
known  form. 

Let  Z/,  Z2',  .  .  .  Zq'  be  approximate  values  of  Z15  Z2,  .  .  .  Z9, 
found  by  trial  or  the  solution  of  a  sufficient  number  of  the 
observation  equations,  and  let  the  most  probable  values  of 
Zt,  Z2,  .  .  .  Zq  be 

ZS  +  ZH   Z2'  +  z2,...Zq'  +  zq  (B) 

Zj,  z2,  .  .  .  zq  being  small  corrections  whose  values  are  to  be 
determined  by  the  Method  of  Least  Squares.  The  first 
observation  equation  in  (A)  may  then  be  written, 


Expanding  the  first  member  by  Taylor's  Theorem,  denoting 
/i(Z,',  Z/,  .  .  .  Zq)  by  &!,  and  neglecting  terms  containing 
powers  of  zn  za,  .  .  .  zq  higher  than  the  first,  this  becomes 

Ok, 


If  now  we  represent  the  coefficients  of   z,,  z2,  .  .  .  ztf,    1  \ 
»  *!»•••  ^i»  also    *i  —  Jl/i   by   m,,    and    treat   the    oilier 
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equations  in  (A)  in  the  same  manner,  our  observation  equa- 
tions will  take  the  form 


The  second  members  reduce  to  the  residuals  since  the 
quantities  zl5  z2,  .  .  .  z7  represent  merely  the  most  probable 
values  of  the  corrections.  Equations  (C)  can  now  be  solved 
in  the  usual  manner  and  the  most  probable  values  of 
Z^  Z.2,  .  .  .  Zq  found  by  substituting  the  results  in  (B). 

45.  Example.  From  the  following  observations  find  the 
most  probable  values  of  x  and  y. 

sin  x  -\-  cos  2y  =  1.5 
cos  x  -f-  3  sin  y  =  1.7  (A) 

a*  +  5y  =  2.1 

By  trial  it  is  found  that  42°  and  18°  are  approximate 
values  of  x  and  y.  Then  in  accordance  with  paragraph  44 
we  put 

1  %j  \         / 

and  expanding  the  different  functions,  we  find 

&!  =  sin  42°  +  cos  36°,      -^  =  cos  42°,      ^--  =  -  2  sin  36° 

=  1.48  =.74  =-1.18 

*2  =  cos  42° +  3  sin  18°,    ^-2  =  -sin42°,    ~=  3  cos  18° 

=  1.67                                    =-.67                 =2.85 
18ir  dlc8 942?r          d  Ar, r 


=  2.11  =1.47  =5 

^  —  JT/!  =  —  .02,          &2  —  J/2  =  -  .03,          £3  —  J/3  =  .01 
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Therefore  we  have  for 

OBSERVATION    EQUATIONS 

J425J  —  1.18z2  —  .02  =  vt 
_  .67z!  -f  2.85z2  —  .03  =  vt  (C) 

1.472l  _|_       5z2  -f  .01  =  vz 

From  these  are  obtained  the 

NORMAL    EQUATIONS 

3.162l  +    4.57z2  +  .020  =  0  D. 

4.57zx  +  34.51z2  -  .012  =  0 

Solving, 

z1==-.  00845  z2  =  .  00147  (E) 

These  results  are  in  circular  measure.  Reducing  to  degree 
measure  we  have 

zl=  -  29'.0  z2  =  5'.1  (F) 

Substituting  these  in  (B), 

x  =  41°  Sl'.O  y=18°  5'.1  (G) 

46.  If  the  observations  are  conditioned,  precisely  the  same 
method  will  be  followed  in  reducing  all  the  equations  to  the 
linear  form.  The  rest  of  the  solution  will  then  be  as  usual. 

If  the  values  found  for  z^  z2,  .  .  .  zq  should  turn  out  to  be 
so  large  that  the  terms  involving  their  second  and  higher 
powers  can  not  be  neglected  as  assumed,  the  process  must  W 
repeated  using  the  values  of  Zt,  Z2,  .  .  .  Zq  first  obtained  as 
approximations. 

In  a  few  cases  the  equations  may  be  reduced  to  the  linear 
form  by  some  special  artifice  of  a  simple*  character,  as  in 
paragraph  43.  In  such  cases  the  method  of  expansion  by 
Taylor's  Theorem  should  not  be  used. 


CHAPTER  III. 

THE   PRECISION  OF   OBSERVATIONS. 

47.  The  work  so  far  considered  has  treated  solely  of  the 
methods  by  means  of  which  the  most  probable  values  of  the 
unknown    quantities  may  be   determined  from  a  series  of 
observations.     But  in  general  something  more  than  this  is 
desired.     We  wish  to  know,  if  possible,  how  much  reliance 
can  be  placed  upon  the  results  obtained,  and  how  they  com- 
pare in   precision  with   other   determinations  of   the   same 
quantities.     Preliminary  to  the  discussion  of  the  precision  of 
our  results  it  will  be  necessary  to  consider  more  fully  the 
Curve  of  Error. 

48.  The  Constant  h.     The  law  of  distribution  of  errors  of 
observation  has  been  shown  to  be   represented  by  a  curve 
whose  equation  is 

y  =  ke-™ 

In  this,  if  x  =  0,      y  =  k  (30) 

Therefore  the  constant  k  represents  the  intercept  of  the 
Curve  of  Error  on  the  axis  of  Y.  It  is  not,  however,  an 
independent  quantity  but  is  determined  by  the  value  of  A, 
as  will  now  be  shown. 

49.  To  Find  k  in  Terms  of  h.     Since,  as  shown  in  para- 
graph   18,  the   total   area   between  the  curve  and   the  axis 
of   X  is  denoted  by  unity,  we  have 


=  1     or     k  C*e-h'W  dx  =  ~          (a) 

*/0  * 

This  may  be  written 
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Let   t  =  Aic,    .-.    dt  =.  hdx.     Also,  when    x  =  GO,  t  =  oc, 
and  when   a;  =  0,   tf  =  0. 


f  "  6-«*  <ft  =  C  e-M*  h  dx  (c) 

»/o  *^o 

Multiplying  this  equation  by 


we  have 


[JT  ^  *  T  =  J]"  J!"  e~"2(1  +  X2 

=  f"-—  -  _  <fo  re-*I 
Jo      2(1  +  x2)       Jo 


1    /••     dx 


Jo  2 

Substituting  (31)   in   (c)  and  (b) 

^          —  k          — 

2  2&  "  y/^ 

Therefore  equations  (16)  and  (17)  become 


(33) 
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50.  It  has  been  shown  that  the  quantity  A    is  a  "Measure 
of  Precision"  of  the  observations,  and  hence  a  determination 
of  its  value  in  each  case  would  enable  us  to  compare   the 
relative  reliability  of  different  measurements.      In   practice, 
however,  it  is  not  found  convenient  to  compute  the  value  of 
h   directly,  and  so  this  quantity  is  used  only  in  developing 
the  theory  of  the  subject,  while  in  the  comparison  of  observa- 
tions certain  other  quantities  now  to  be  derived  are  used  as 
precision  measure's.     These  latter  quantities  are  called  the 
"Mean  of  the  Errors,"  the  "Mean  Error,"  and  the  "Probable 
Error,"  respectively,  and  maybe  computed  directly  from  the 
observations.     Further,  as  it  will  be  found  that  they  all  bear 
a  definite  relation  to   A,   the  value  of  this  quantity  can  be 
determined  from  them  if  desired. 

In  the  following  discussions  no  distinction  is  made  between 
positive  and  negative  errors  of  the  same  numerical  magnitude, 
and  unless  otherwise  stated  the  observations  are  all  of  the 
same  weight. 

THE  MEAN  OF  TKE  ERRORS  OR  AVERAGE  DEVIATION. 

51.  The  Average  Deviation    (a.d.)   of   an  observation  is 
the  arithmetical  mean  of  the  errors  all  taken  with  the  positive 
sign. 

Since  from  (6)  the  probability  that  the  error  of  a  single 
observation  will  fall  between  x  and  x-{-dx  is 

<f>(x)  dx, 

if   n   observations  are  made,   the  number  of  errors  falling 
between  these  limits  is 

n  4>(x)  dx. 
Hence  the  sum  of  all  the  errors  of  the  observations  is 

x»oo  ~» 

n  I      x<j>(x)  dx      or       2n  |      x  <j>(x)  dx. 
J-x  Jo 
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Dividing  this  last  expression  through  by    n    we  have 

/»*>  2A   /»» 

a.d.  =  2  I     x+(x)  dx    =    —  |     e~h^  'x  dx 

•Jo  VTT^O 

=  -  -^z  f  V*2*'  (-2A2jc)  <fo 
AVir^o 

1 

or  a.d.  =  —  (34) 


THE  MEAN   ERROR. 

52.  The  Mean  Error  (/A)  of  an  observation  is  the  square 
root  of  the  arithmetical  mean  of  the  squares  of  the  errors. 

The  total  number  of  errors  being  n,  the  number  falling 
between  x  and  x-\-dx  is,  as  just  shown, 

n  <j>(x)  dx, 
and  the  sum  of  the  squares  of  these  errors  is 

n  x1  $(x)  dx. 
Therefore  the  sum  of  the  squares  of  all  the  errors  is 

n  C    x2  <t>(x)  dx 

J—a» 

>  dx  a 


But  as  shown  in  paragraph  49, 

dx==l       or        r  e-hW  dx  =  ^.       (b) 

-La, 

Differentiating  (b)  with  respect  to   A, 

-2A  ^Vefr;   =   -~  (c) 
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and  replacing  the  integral  in  (a)  by  its  value  as  determined 
from  equation  (c),  we  have 

1 


THE   PROBABLE   ERROR. 

53.  The  Prohable  Error  (r)  of  an  observation  is  an  error 
such  that  one-half  the  errors  of  the  series  are  greater  than  it 
and  the  other  half  less  than  it.  Or  it  is  an  error  of  such  a 
magnitude  that  the  probability  of  making  an  error  greater 
than  it  in  any  given  observation  is  just  equal  to  the  probability 
of  making  one  less  than  it,  both  probabilities  being  one-half. 

The  probability  that  the  error  of  an  observation  will  fall 
between  x  and  x-\-dx  being  <j>(x)  dx,  the  probability 
that  the  error  will  fall  between  the  limits  r  and  —  r  is 


P  =    Cr4>(x)dx    =        r  Ce-Wdx          (36) 

J-r  yV^o 

If   r   is  the  probable  error,  P   is  one-half,  or 

»w«fe  =  L  (37) 


and  from  this  definite  integral   r   is  to  be  found. 

Let  t  =  hx,  .'.  dt  =  hdx.  Also  when  x  =  r,  we  have 
t  :=  Ar,  and  when  x  =.  0,  t  =  0.  Substituting  these  results 
in  (37),  we  have 

1  (38) 
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Denote  hr  by  p.  Then  by  interpolation  in  a  table  of 
values  of  this  integral,  the  value  of  hr  in  (38)  is  found  to 
be 

p  =  hr  =  .47694  (39) 

p          .47694 


54.    When   t  is  small  the   values   of    j     e~^  dt    may  be 

° 
found  by  expanding  e~*   into  a  series  and  integrating   the 

successive  terms.     Thus,  by  Maclaurin's  Theorem, 


/•<  /4  f« 

dt   =    I     (1  _  p  +    —  -  —  +  ..  .  )  dt 
Jo   v  H         IL 

«      I8     J!_ 


3         5[^        7[3_ 


RELATIONS  BETWEEN    |i,   r,   a.d.,   Ik,    AND    p. 

55.    From  (40)   r  =  -,  and  from  (34)   a.rf.  =  -  , 


r  =  p  a.d.  \f* 
or  r  =  ,8453a.^.  (41) 

Also  from  (35)    ^  =  -  , 


r  =  Pn\ 
or  r  =  .6745  \i  (42) 
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The  relation  between  the  values  of  /x,  r,  a.d.,  and  h  may 
be  conveniently  expressed  as  follows:  — 

u.^2   ===  "~  ==  ~  ==   a'^'  0f  (^) 

or  arranging  in  tabular  form, 

p  r  a.d. 

p         1.0000         1.4826         1.2533 

r        0.6745         1.0000         0.8453 

a.d.        0.7979         1.1829         1.0000 

From  this  it  will  be  seen  that 

p.  >  a.d.  >  r  (44) 

56.  Further,  since  in    (25)   it  was  shown  that  p  <x  h\    it 
follows  at  once  from   (43)   that 

i        JL  QC  JL  oc  fc« 

a.d.*        K-2         r2 

That  is,  the  weights  of  different  determinations  of  a 
quantity  vary  inversely  as  the  squares  of  their  Mean  Errors, 
their  Probable  Errors,  or  their  Average  Deviations. 

It  is  to  be  observed,  however,  that  the  determination  of  the 
relative  weights  of  quantities  from  a  comparison  of  their 
precision  measures  according  to  (45)  applies  only  when  the 
quantities  are  of  the  same  kind  and  subject  to  the  same  con- 
stant errors,  if  any  of  the  latter  exist.  (See  §  98.)  The 
applications  of  (45)  are  numerous  and  important. 

57.  Example  A.     Suppose   n   direct  observations,  all  of 
the  same  weight,   be   made   upon   a  quantity,   and   that   the 
probable  error  of  a  single  observation  is   r.     Then  since  the 
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weight  of  the  arithmetical  mean  is    n,    its  probable  error    r0 
will  be  given  by 

*     I     or     ,.      '  (46) 


Or  in  general,  suppose  8  is  any  precision  measure  of  an 
observation  of  weight  />,  and  suppose  p0  is  the  weight  of  a 
second  similar  quantity  or  observation,  then  the  correspond- 
ing precision  measure  S0  of  the  latter  will  be 


8°  =  8 


The  case  of  most  common  occurrence    is   that  in  which 
p  =  1,  and  then  we  have 

80  =  -=  (48) 


Example  B.  A  line  is  measured  five  times  and  the  aver- 
age deviation  of  the  mean  (A.  D.)  found  to  be  .016  feet. 
How  many  additional  measurements  are  necessary  in  order 
that  the  A.  D.  of  the  mean  may  be  reduced  to  .004  feet  ? 

Let  x  be  the  total  number  of  observations  required.    Then 

x   :    5  =  .000256    :    .000016 

x  =  80 

Consequently  the  number  of  additional  measurements 
required  is  75. 

Example  C.  In  two  determination!  of  the  quantity  L 
there  were  obtained 

1^  =  427.320  ±  0.040,         Z,  =  427.30  ±  0.16 

Find  tln-ir  relative  weights,  and  the  most  probable  v:ilue 
of  L  and  its  probable  error. 
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Note.  The  above  is  the  method  commonly  employed  to 
denote  that  the  probable  errors  of  the  observations  are  0.040 
and  0.16. 

pl         162  16 

From    (45)  g  =     —          ~ 

From    (3),  the  most  probable  value  of   L   will  be  given  by 

A  =  427   +   16  X  "g  +  -30 

=  427.319 
and  the  weight  of    X0   being    17,   by    (48), 

r0  =   -^   =  .039 

vTT 

Therefore  we  should  write  the  result 

Z0  =  427.319  ±  .039 


REPRESENTATION  OF  \L,  a.d.,  AND  r  ON  THE  CURVE  OF 

ERROR. 

58,    To  find  the  points  of  inflection  of  the  Curve  of  Error 
we  have 

y  =  ke-M# 

-  =    -  2  h*k  x 


For  a  point  of  inflection,      .  ].,    =    0,   or 

2  A2  A  e-^2  (2  Aa  sua  -  1)  =  0 

"     by(35)" 
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That  is,  the  Mean  Error  is  represented  by  the  abscissa  of 
the  point  of  inflection  of  the  Curve  of  Error.  See  OM  in 
figure  2,  page  9. 

Next,  for  the  abscissa  of  the  centre  of  gravity  of  the  area 
to  the  rifirht  of  0  Yt  we  have 


I     y  x  dx 
=  OD  =  ^r- 

l    y  dx 


7= 

VTT 


for,  (§  18),    Cy  dx=± 

«/0  L 

Integrating, 

a;0  =  —  —  =    a.d.  by  (34). 


Finally,  if  an  ordinate  PP'  be  drawn  so  as  to  bisect  the 
area  to  the  right  of  the  origin  between  the  curve  and  the 
axis  of  X,  the  Probable  Error  will  be  represented  by  the  dis- 
tance of  this  ordinate  from  the  axis  of  Y.  For  the  proba- 
bility of  the  occurrence  of  an  error  less  than  the  amount  OP 
is  then  equal  to  the  probability  of  the  occurrence  of  an  error 
greater  than  OP.  This  being  the  case,  OP  is  the  probable 
error  by  definition. 


CHAPTER  IV. 

COMPUTATION  OF   THE   PRECISION  MEASURES. 

DIRECT   OBSERVATIONS. 

59.  Observations  of  Equal  Weight.  Given  n  direct  obser- 
vations all  of  the  same  weight  on  a  single  quantity  M,  to 
find  the  Mean  and  Probable  Errors  and  Average  Deviation 
of  a  single  observation  and  of  the  Arithmetical  Mean. 

Let  the  observations  be  M^  M^  .  .  .  Mn. 

"•      "    arithmetical  mean  be  3IQ. 

"      *'    real  errors  be  xlt    ic2,  .  .  .    xn. 

"      "    residuals  be  vl5    v2,  ...    vn. 

Denote  the  mean  and  probable  errors  and  average  devia- 
tion of  a  single  observation  by  /n,  r,  and  a.d.,  respectively, 
and  the  corresponding  quantities  for  the  arithmetical  mean 
by  /MO,  7*0,  and  A.D.  Then  by  definition, 

Vj.    =    MI     —     MQ,  V2    =    MI    —     M0,    ...  Vn   =    Mn    —     MQ 

\*& 

and  ^»  =  \  — 

If  MQ  represented  the  true  value  of  M,  the  residuals 
would  be  the  same  as  the  real  errors,  and  we  should  have 


/ 

!\ 


and  if  n  is  large  this  formula  is  practically  exact.     But  when 
n   is  small    a    more   accurate   expression    is    necessary.     To 
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obtain  this  let   3f0  +  «o   be  the  true  value  of   M.     Therefore 

«<>)    =   vl  —  XQ 


ic0     =  vn  - 
Siquaring,  adding,  and  dividing  by   n 

-  2 


=  +  «o2        for  by  (2)    Sv  =  0 

The  value  of    je0    is    not  known  and  can  not  be  found 
exactly,  but  it  is  approximately  equal  to  the  mean  error  of 

Jf0,   that  is,  by   (48),    to   ^   =   Ji.. 

^n 
Substituting  this  value  in  the  above  we  have 


n          n 

(n  -  \)p?  =   2u2  (49) 


by  (42)  ,.  =  .6745       ---  (52) 

and  n  =  .6745 
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60.    In    order   to    avoid   the    use   of  the    squares   of   the 
residuals  we  may  proceed  as  follows  :    From   (49) 


On  the  average,  the  values  of  the  residuals  will  then  be 


n-  1 


Adding  and  dividing  by  n,    neglecting  the  signs  of  the 
residuals, 

2u  I  n  —  I  "S,x  In  -^~T 

n  \       n  n  \        n 


ff.ff.   = 

by  (48)  A.D.  = 

by  (41)  r  = 

and  rn  = 


Vtf«; 
(55) 
(56) 
/•S7^ 

\ln  (n—  1) 

Zv 

n  ^  n  —  1 
.8453  Zt? 

<Jn(n-  1) 
.8453  Zr 

n  Vti  —  1 


The  mean  errors  may  also  be  computed  from  the  above  by 
using  the  table  of  equivalents  in  paragraph  55,  but  this  is 
not  customary,  formulas  (50)  and  (51)  being  used  for  this 


50 
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purpose.     Results  derived  from    (50)    are  to  be  regarded 
more  accurate  than  those   obtained  from    (54),    the    latter 
being  a  second  approximation. 

61.  Example.  From  the  following  measurements  on  the 
length  of  a  base  line  find  the  most  probable  length  and  the 
values  of  the  various  precision  measures. 


M 

V 

w2 

455.35 

.02 

.0004 

.35 

.02 

4 

.20 

-  .13 

169 

.05 

-  .28 

784 

.75 

.42 

1764 

.40 

.07 

49 

.10 

-  .23 

529 

.30 

-  .03 

9 

.50 

.17 

289 

.30 

-.03 

9 

3.30 

+  .70 

.3610 

455.330 

-.70 

M. 

2v  =  0 

2tf2 

By  (50)  and  (51) 


.3610 
9 

r   =    .6745  n 
By  (54)  and  (55) 
1.40 


.20  uo     =     -E-     =  .063 


.13  r0  =  .6745  ^  =  .042 


=  -15      A-D-  -  w  =  -047 


r  =  .8453  a.rf.  =  .13          r0  =  .8453  ^.D.  =  .040 
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and  we  should  write  for  the   most    probable  length  of  the 
base  line 

Jf0   =   455.330  ±  .042 

62.    Observations   of  Unequal  Weight.     Using  the  same 
notation  as  above,  with  slight  modifications,  we  will 

Let   MO   represent  the  General  Mean. 
"    P\i  Pit  •  •  •  Pn   be  the  weights. 
"     a.d.t,  a.d.i,  .  .  .  a.d.n   be  the  average  deviations. 
"     /*!?   /*a»  •  •  •  /*n  be  the  mean  errors. 
"     rit    r<f>  •  •  •  rn   be  the  probable  errors. 
"     a.d.,  fr  r,  and  v  refer  to  observations  of  weight  unity. 

Then  by  (48) 

i/  11 

etc. 


If  the  "  Observation  Equations  "  are  formed  for  this  case 
they  will  be 

MI  —  M0  =  VD     M2  —  MQ  =  v2,  .  .  .  Mn  —  MQ  =  vn. 

And,  as  was  shown  in  paragraph  28,  if  these  equations  are 
each  multiplied  by  the  square  root  of  the  weight  of  the  corre- 
sponding observation,  they  will  all  be  reduced  to  equivalent 
equations  of  weight  unity.  On  performing  this  operation 
it  will  be  seen  that  the  residuals  of  the  new  equations 
become 


And  evidently  to  these  reduced  observations  the  formulas 
of  paragraph  59  apply. 
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Therefore 


=    V/      n-"l  r  =  l6745  •*         (58) 


^  V 


rt  =  .6745  |i.         (59) 


M.   =    V  v^lD         n==  .««m        (60) 


Also,  by  a  method  similar  to  that  used  in  paragraph  60,  it 
may  be  shown  that 

a.d.  =  — Sl?    ^  r  =  .8453 a.d.       (61) 

\  n(n—  1) 


a.d.k  =  -  rk  =  .8453  a.d.A.     (62) 


r0  =  .8453  ^i.D.     (63) 


Formulas  (61),  (62)  and  (63)  are  not  of  much  value  in 
practice  unless  all  the  weights  are  perfect  squares,  for  other- 
wise no  labor  is  saved  in  the  computations,  since  the  square 
root  of  each  weight  will  have  to  be  determined,  and  the  results 
obtained  can  not  be  considered  as  reliable  as  those  given  by 
formulas  (58),  (59)  and  (60). 
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63.  Example.  Given  a  series  of  observations  on  Jf/J 
difference  in  longitude  between  two  stations.  To  find 
/x,  /IQ,  a.d.,  A.D.,  r,  r4,  r0. 


the 


M 

p 

pM 

V 

v^p 

V2 

pv* 

1°  4'  30" 

4 

120 

1.6 

3.2 

2.6 

10.4 

41 

1 

41 

9.4 

9.4 

88.4 

88.4 

43 

1 

43 

11.4 

11.4 

130.0 

130.0 

37 

9 

333 

5.4 

16.2 

29.1 

261.9 

48 

4 

192 

16.4 

32.8 

269.0 

1076.0 

34 

16 

544 

2.4 

9.6 

5.7 

91.2 

25 

9 

225 

6.6 

19.8 

43.6 

392.4 

46 

1 

46 

14.4 

14.4 

207.4 

207.4 

28 

25 

700 

3.6 

18.0 

13.0 

325.0 

24 

4 

96 

7.6 

15.2 

57.7 

230.8 

1°  4'  31".6 

74 

2340 

150.0 

2813.5 

M0 

*P 

2pm 

S*^ 

2pv* 

=    1°  4'  31".6 


=  -£..=  2.1 


=   .6745  p.  =   11.9 
15.8 


r     = 


150 

V/90 

.6745  uft 


r.  -    =   4.0 

X.A==^==   1.8 


1.4     =     .8453  A.D.   =   1.5 
=    1°  4'  31".   6    ±    1".4 


54  METHOD  OF  LEAST  SQUARES. 

FUNCTIONS   OF    OBSERVED    QUANTITIES. 

64.  Theorem.     Given  any  number  of  quantities  and  their 
Mean  and  Probable  Errors  and  Average  Deviations,  to  find 
the  Mean  and  Probable  Errors  and  Average  Deviation  of  any 
function  of  the  quantities. 

Let  the  quantities  be        JH/j,  Mz,  .  .  .  Mq. 
"      "    mean  errors  be     j^,  /^,  .  .  .  pq. 
"      "    function  be  M=f(M»Mi,  .  ..  Mq). 

"      "    mean  error  of      M  be  E. 
"      "    probable  error  of  M  be  R. 
"      "    average  deviation  of  M  be  D. 

The  derivation  of  the  general  formula  will  be  simplified  if 
we  consider  first  a  few  special  forms  of  functions. 

65.  Case  I.     Suppose    M  =  MI  ±  M2. 

The  number  of  observations  from  which  M±  and  M^  and 
hence  ^  and  /^  have  been  determined  is  not  necessarily 
known,  but  we  may  assume  that  for  each  quantity  it  is  any 
large  number  n,  and  that  the  real  errors  of  the  observations 
are 

for         J/i,     */,  */',  V,  .  .  . 
Mz,     a,',  a*",  *,'",  •  •  . 

Then  the  real  errors  of  My  computed  from  the  separate 
observations  on  M  and  M^  will  be 


or  E>  --     m*+m«  (64) 

since  in  the  most  probable  case  the  term   Sxjg,   will  dis- 


appear, as  there   most  likely  will  be  as  many  positive  as 
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negative   products  of  the   same  absolute  magnitude   of  the 
form    xl  ic2. 

By  successive  applications  of  the  above  the  same  principle 
may  be  extended  to  cover  the  algebraic  sum  of  any  number 
of  quantities.  So  that  if 

M  =   Ml±Mz±...  Mq 
then       &  =   ^  +  n,2  +  .  .  .  }V   =  Sfi2  (ttf) 


Since  probable  errors  and  average  deviations  differ  from 
errors  merely  by  a  constant  factor,  we  shall  likewise 
have 

12*  =   ri*  +  ?V  +  .  .  .  r*   =  Sr2  (66) 

D2  =   a.d.?  +  a.d.S  +  .  .  .  =  Sa.tf.2      (67) 

66.  Example,     Given  the  telegraphic  longitude  results, 

h.  m.  sec.           sec. 

(A)  Cambridge  west  of  Greenwich,      4  44  30.99  ±  0.23 

(B)  Omaha            "      "Cambridge,       1  39  15.04  ±  0.06 

(C)  Springfield  east    "  Omaha,  25  8.69  ±  0.11 

Find   L,,   the  longitude  of  Springfield,  and  its  probable  error. 

L   =   A  -f-  B  —  C  =   5*  58m  37*.34  ±  0*.26 
for  by  (66)          7?   =    y/(.232  +  .062  +  .II2)    =   .26 

67.  Case  IT.     Suppose     Jf£"  =    ai-Mi 

Using  the  same  notation   as  in    paragraph    64,    the    real 
errors  of  J/  will  be 


n 


= 


or  JS    =    ai^i  (68) 
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68.     Case  III.    Suppose   M=  a^Mi  ±  azM2  -  - .  aqMq 
By  combining  (65)  and  (68) 


also  J22  =   Sa2r2  (69) 

and  J>2  =   S  a2  a.d.* 

69.  Example.     The  length  of  a  bar  at  20°  Centigrade  is 
found   to   be   75".0041    ±    0".0037,    and  its   expansion    per 
centigrade  degree  0".0036  ±  0".0018.     What  is  the  length 
of  the  bar  at  56°  Fahrenheit  ? 

20°  C.    =    -  X  20  +  32    =    68°  F. 
5 

5  5 

The  expansion  for  1°  F.  will  be   •-  X  .0036  ±  -  X  .0018 

9  9 

or  .0020  ±  .0010 

L  =  75.0041  -  12  X  .0020 

=  74.9801 
E  =  V.00372  +  (12  X  .001 )« 

=  .013 
L  =  74".980  ±  0".013 

70.  Case  IV.     Suppose   M  =  f  (Jtfi,   M»  .  .  .  Mq). 

Let  MI  =  al  -f-  mlt  M^  =  aa  -f-  wia,  . . .  Mq  =  aq  -\-  m^ 
where  at,  as,  .  .  .  aq  are  arbitrarily  assumed  quantities  very 
nearly  equal  to  JHJ,  M^  .  .  .  Mqt  so  that  mly  m^  .  .  .  niq  are 
so  small  that  their  second  and  higher  powers  may  be  neg- 
lected. Then  the  errors  of  Jf/i,  M^  .  .  .  Mq  may  be  con- 
sidered to  be  in  mlt  m^  .  .  .  mq,  and  hence  /it,  /x.,,  .  .  .  /xy 
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may   be    regarded    as   the    mean    errors    of   the    quantities 
/;/,,  >/>.,,  .  .  .  mq.     We  now  have 

M  =  f  (ai  +  mi»  «2  +  ™2,  •  •  •  «?  +  ™q) 

Expanding    this    expression    by   Taylor's    Theorem,    and 
denoting  /(«i,  «2,  .  .  .  «9)  by  M',  we  have 

M  —   M' 

dM'  dM1 

-  mi  •  -  ™2  •  -  '  '  '  m 


-|-  negligible  terms   in    the  second   and 
higher  powers  of  mly  m2l  .  .  .  mq. 

Then  the  mean  error  of   M  will  be  the  same  as  the  mean 
error  of  terms    (A),    and  by    (69)    this  is  given  by 


or  what  is  practically  the  same  thing,  by 


71.    Example  A.     Two  sides  of  a  right  triangle  are  meas 
ured  with  results 

a    =    49.53  ±  0.59 
b    =    50.38  ±  0.93 


58 
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Find  the  length  of  the  hypotenuse    c    and   its  probable 
error.     In  this  example 


and 


r49.53~l2  r50.38~l* 

=    .592      4-  .932 

[_70.65j  |_70.65j 

—    .78 

c    =    70.65  ±  0.78 


.Example  B.  If  the  probable  error  of  x  is  r,  what  is 
the  probable  error  of  the  common  logarithm  of  x  ?  In  this 
case 

M  =   Iog10  ^     =    Iog10  x 

d  Iog10  x  Iog10  e 

dx 


Example  C.  If  the  weight  of  x  is  />,  what  is  the  weight 
Po  of  sin  x  ? 

Denoting  the  mean  errors  of  x  and  sin  x  by  /t  and  -£*, 
respectively,  we  have,  by  (45)  and  (70), 

E*  =   ua  =   ft2  cos*  as 

c>* 
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Po    =     t^_ 
J)  j& 

P 

Po    =    — —    =    P  sec-1  x 
cos1*  x 

72.  Equation   (70),  which  expresses  the  law  of  propaga- 
tion of  error  in  functions  of  observed  quantities,  is  one  of  the 
most  important  in  the  whole  theory  of  the  Method  of  Least 
Squares.     Upon  it  in  particular  is  based  the  discussion  of  th* 
"  Precision  of  Measurements."     This   subject  treats,  in  the 
first  place,  of  the   methods  of   finding  the   precision    of   a 
quantity  obtained  by  computation  from  a  series  of  measured 
quantities  ;  snd,  in  the  second  place,  it  investigates  the  pre- 
cision with  which  the  component  measurements  of  a  series 
must  be  made  in  order  to  obtain  a  required  degree  of  pre- 
cision in  the  final  result.     The  following  simple  example  win 
illustrate  the  character  of  the  solutions  :  — 

73.  Example.     In  the   determination  of  a  current  bv  a 
tangent  galvanometer  we  have 

rr 

I  =    10  --  tan  <£ 
G 

where  I  is  the  current  in  amperes,  II  the  horizontal  com- 
ponent of  the  earth's  magnetic  force,  G  the  galvanometer 
constant,  and  <f>  the  angle  of  deflection.  Given  the  errors 
81,  §25  83,  in  -f-fy  G  and  tan  <£,  to  find  the  error  A  in  I. 

By  (70) 

A2  = 


100 


60  METHOD   OF  LEAST  SQUARES. 

Dividing  this  equation  by   I*  =   -          -  tan2  <£   we  have 


That  is,  the  square  of  the  percentage  error  in  I  is  equal 
to  the  sum  of  the  squares  of  the  percentage  errors  in  J3J  6r, 
and  tan  <£.  Hence  if 

H  is  determined  within  .4  per  cent. 

Q     «  a  a  ^2       "         " 

tan  <f>   "  "  "        .1     "       " 


then  —  =   v/  .16  -+-  .04  -f  .01   =   .46  per  cent.      (c) 

Next,  suppose  the  value  of  I  is  required  to  within  .1  per 
cent.  To  find  the  necessary  accuracy  in  the  determinations 
of  H,  (r,  and  tan  <£  when  the  error  in  each  of  these  quan- 
tities is  to  have  the  same  influence  upon  the  total  error. 

From  (b)  we  shall  now  have 


=   —  =   .000577 


8,   =   .00058^ 

Sa  =   .00058  Q  (e) 

83   =   .00058  tan  <£ 

It  is  comparatively  easy  to  obtain  the  necessary  accuracy 
in  the  measurements  of  O  and  tan  <£,  but  difficult  in  the 
case  of  77. 

For  additional  work  of  this  kind  see  Ilolman's  u  Precision 
of  Measurements." 
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74.  Combination  of  Functions  of  the  Same  Variables. 
It  is  to  be  noticed  that  equation  (70)  applies  only  when  M 
is  a  function  of  independent  quantities.  If  M^  M^  .  .  .  Mq 
are  merely  different  functions  of  the  same  quantities  we 
must  proceed  as  follows  :  — 

Let  Mi    =    <£  (zx,  z2,  .  .  .  z,t) 

Mz    =    $  (Zu  z2,  .  .  .  z*) 
M    = 


If  any  single  observations  of  z1?  z2,  .  .  .  2fc  are  subject  to 
errors  JCj,  ic2,  .  .  .  ic^  the  corresponding  errors  in  Ml  and 
Jlf  will  be 


for   Mly  XL    =    a^  -|-  azxz  -\-  .  .  .  akxk  (a) 

"     M2y          A^    =    oj'aj!  -j-  «2'z2  +  •  •  •  ak'xt  (b) 


Where  «!,  «2»  •  •  •  ak  are  the  differential  coefficients  of 
J/i,  and  a/,  a2'  .  .  .  at;  the  differential  coefficients  of  M^ 
with  respect  to  z^  z2,  .  .  .  zk.  The  corresponding  error  in 
M  will  then  be 

X  =    AX,  +  A'X,  (c) 

Where    A   and    ^1'    are  the  differential  coefficients  of   M 
with  respect  to    Jl/i    and    3f2.     Substituting  in  (c)  from  (a) 
and  (b), 

X  =    (Aai  +  A'a\)  xv  -\-  (Aaz  +  Aa\)  «2  .  .  . 
.  .  .   \Xk 


Then  if  the  number  of  observations  or  values  of   X  be 
denoted  by   w, 


u 


.  XV*2  (d) 
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since   in  the   most   probable   case   the   product   terms   will 
cancel  out. 

Expanding  (d)  we  have 

E*  =    (Aa,  +  A'a,')*  rf  +  (Aa,  +  A'a,'}2  tf  +  .  .  . 

^'2KV  +  «.'V  •  •  •) 

...)  (70 


75.    Example.     As  a  very  simple  problem  take 

J/i    =    22i,  Jf2    =    325!,  /n    =    0.1, 

and  Jf  =  Jlft  +  J/2. 

Then      ^4  =  1,         ^4'  =  1,         at  =  2,         at'  =  3. 
By  (71)  ^2  r=  4  X  .01  +  9  X  .01  +  2  X  2  X  3  X  -01  =  .25 
or  E   =    0.5 

In  this  particular  example  the  result  may  be  found  directly 
from  (68)  by  substituting  at  first  in  M  the  values  of  M^ 
and  JH/j. 

Thus  M  =   22X  -f  32!    =    52t 

E  =    5^   =    0.5 

If  J/i  and  J/2  had  been  independent  quantities,  by  (64) 
or  (69)  we  should  have  had 


E  =  V(2  X  0.1)a+  (3  X  0.1)a 
=   0.36 

INDIRECT    OBSERVATIONS. 

76.  The  determination  of  the  precision  measures  of  the 
unknown  quantities  in  case  the  observations  are  indirect 
involved  a  knowledge  of  the  weights  of  the  unknowns,  and 
consequently  the  method  of  computing  these  weights  must 
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first  be  demonstrated.     It  will  be  assumed  at  present  that  all 
the  observations  are  of  weight  unity. 

FIRST  METHOD  OF  COMPUTING  THE  WEIGHTS. 
Let  the  observations  be 


in  which  M^  M»,  .  .  .  Mn  denote  the  actual  observations, 
and  2j,  22,  ...  zq  the  most  probable  values  of  the  unknown 
quantities.  Let 

k1  —  3/t  =  m1?     &2  —  Mz  =  7??2,  ...  kn  —  Mn  =  mn 
Then  the  above  equations  give  rise  to  the 

OBSERVATION    EQUATIONS 


«s«i  +  fat  +  •  •  •   ^2^  +  ^2    =    ^ 

(A) 


By  the  rule,  paragraph  25,  we  now  form  the 

NORMAL    EQUATIONS 

zl  2  a2  -f-  22  2  ai  +  .  .  .  zq  2  aq  -f-  2  «m    =    0 
2!  2  a&  +  22  2  *2   +  .  .  .  zq  2  bq  -f  2  ^>m    =    0 

(B) 

2!  2  a?  -f-  za  2  6?  +  .  .  .  zq  2  q*    +2  </w    =    0 
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Multiply  the  first  of  (B)   by   &,  the  second  of  (B)   by 

.  .  .  and  add  the  results.     Then 


...  Qqlaq) 

+  Q^V   +  ... 


-f-  &S  am  -f-  §22  bm 
=    0 


Let  Q19  QZ,  ...  Qq  be  determined  so  that  the  coefficient 
of  zl  in  (C)  shall  be  unity,  and  the  coefficients  of  za, 
ZH  ...  zq  each  equal  to  zero.  That  is,  let 


+  §22ai  +  .  .  . 

=    0 


...    Qq*q*  =    0 

Then  (C)  becomes 

zv  +  Q^am  +  ftSJm  +  .  .  .  QqX  qm    =    0     (E) 

Equations  (D)  may  be  derived  from  the  normal  equations 
(B),  if  in  them  we  replace  zt,  z2,  ...  zq  by  Q^  §„,...  Qq, 
2  am  by  —  1,  and  2&m,  2cm,  ...  2qm  by  zero.  Hence 
the  solution  of  the  normal  equations  with  these  changes  will 
give  the  values  of  Q19  (>2,  ...  Qq. 

TO  SHOW  THAT   Qi  IS  THE  RECIPROCAL  OF  THE  WEIGHT  OF  XL 
Expanding  the  coefficients  of  (E)  we  have 
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=    0 
or  collecting  the  coefficients  of    mly  m2,  .  .  .  mn 


=    0 
For  convenience  in  writing  we  will  let 


ai    = 

a,     = 


(E') 


.  .    Qqqn 

=      0  (G) 


Multiply  the  first  of    (F)    by    al9    the  second  by    a2,  .  .  . 
and  add  the  results.     Then 


=    1  by  first  of    (D)  (H) 

Multiply  the  first  of    (F)    by    bly    the  second  by    62,  ... 
and  add  the  results.     Then 


66  METHOD   OF  LEAST  SQUARES. 


=    0  by  second  of    (D) 

Likewise         2  ca    =    0  (H') 


2?a    =    0 

Multiply  the  first  of  (F)  by  04,  the  second  by  o^,  .  ,  . 
and  add  the  results.  This  gives 

2a2     =     &S  Oa  +   ftS  Aa  +   .  .  .  Qg^  ga 

=    Q,          by  (H)  and  (H')  (I) 

Let    ft     be  the  mean  error  of  an  observation  of  weight  unity. 
Let    ftg,  be  the  mean  error  of    z^ 
Let   pZl  be  the  weight  of    z±. 

The  mean  errors  of  m^  m2,  .  .  .  mn  are  the  same  as  the 
mean  errors  of  M±,  M^  .  .  .  Mn  and  each  is  accordingly 
equal  to  ft.  Therefore  from  (G),  by  (69) 


*y  (i)  (J) 

But  by  (48)    ftzi'    =    £ 

Pzi 

Comparing  this  with  (  J)   we  see  at  once  that 

«•  =  s  (K) 

Therefore,  for  the  First  Method  of  computing  the  weights 
we  have  the  following  :  — 
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77.  Rule  I.    In  the  normal  equation  for   zl   write  —  1 
for  the  absolute  term    Sara,    and  in    the  other  equations 
zero  for  each  of  the  absolute  terms    2  bm,  2  cm,  ...  2  qm. 
The  value  of  zl  found  from  these  equations,  is  the  recipro- 
cal of  the  weight  of  the  value  of  zl    obtained  by  the  solution 
of  the  normal  equations. 

To  find  the  weights  of  z%,  z8,  .  .  .  zg,  proceed  in  a  similar 
?/wy,  forming  a  corresponding  set  of  equations  for  each 
unknown. 

SECOND  METHOD  OF  COMPUTING  THE  WEIGHTS. 

78.  Write  equations    (B)    of  paragraph  76  in  the  form 

z-i  2  «2  +  zz  2  ab  -f-  ...  *q  2  aq  -\-  2  am    =    A 

zl  2  ab  +  z2  2  b2   +  .  .  .  zq  2  bq  +  2  bm    ==    B 

(L) 


Then  in  the  solution  by  the  preceding  method,  equation 
(E)  becomes 


=    0i^  +  C.l?+  •••   0,0  (M) 

in  which,  as  was  proved  in  (K),  Q±  is  the  reciprocal  of  the 
weight  of  Zi.  Whatever  method  of  elimination  is  employed 
in  the  solution  of  the  normal  equations,  the  coefficient  of  A 
in  the  value  of  z±  must  necessarily  be  always  the  same. 
Hence  we  have 

79.  Rule  II.  Write  A,  B,  .  .  .  Q  instead  of  zero  in  the 
second  members  of  the  normal  equations  and  carry  out 
(//>  ir  solution  in  any  convenient  /'•">/.  Then  the  most 
probable  values  of  Zj,  z2,  .  .  .  z  are  given  by  those 
in  the  results  v/V^VA  are  independent  of  A,  7?.  . 
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The  weight  of  z±  is  the  reciprocal  of  the  coefficient  of  A 
in  the  value  of  z^  The  weight  of  22  is  the  reciprocal  of 
the  coefficient  of  B  in  the  value  of  z2,  etc.,  etc. 

THIRD    METHOD    OF    COMPUTING    THE    WEIGHTS. 

.  80.  From  the  second,  third,  .  .  .  equations  of  (L)  find 
the  values  of  22,  28,  ...  zq  in  terras  of  z±  and  substitute 
in  the  first  of  (L)  without  reduction.  Then  the  first  of 
(L)  becomes 

Rz^    =     T  +  A  +  terms  in  J?,  <7,  ...  Q 

Where    T    is  the  sum  of  all  the  numerical  quantities  result- 
ing from  the  substitutions.     Dividing  through  by   jR, 


T 

in  which    -—    is  the  most  probable  value  of    21?    and,  as  was 
ft 

shown  in  deriving  the  second  method, 

!-'.*• 


R  =  P»  (O) 

From  this  follows  at  once 

81.  Rule  III.  Substitute  in  the  normal  equation  for  zl 
the  values  of  za,  z8,  .  .  .  zq  in  terms  of  z^  as  found  from 
the  remaining  equations.  Then  before  freeing  of  fractions 
or  introducing  any  reduction  factor,  the  coefficient  of  ^ 
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in  this  equation  is  the  weight  of  the  value  of  zt  obtained 
in  the  solution. 

To  find  the  weights  of  22,  23,  .  .  .  zq,  proceed  in  a  similar 
way  with  the  normal  equations  for  each  of  these  unknowns. 

For  the  solution  of  an  example  by  the  three  different 
im-thods  see  paragraph  84. 

THE  MEAN  ERROR  OF  AN  OBSERVATION. 

82.  The  next  step  will  be  to  derive  /*,  the  mean  error 
of  an  observation  of  weight  unity.  In  the  following  demon- 
stration the  equations  referred  to  by  letters  are  those  in 
paragraphs  76  to  81. 

Let  the  real  values  of    zt,  22,  ...  zq    be 

2!    -f    »!,          22    +    JBj,    .   .   .          Zq    +    Xq 

and  substituting  in    (A)  we  have 


xq)  +  m^  =  A2 
(P) 


where   A1?  A2,  .  .  .  An  are  the  real  errors  of    J!/i,  J!/2,  .  .  .  J^M 

Or  of     TWj,  7712>    •  •  •    mn« 

Multiply  tlie  first  of    (P)    by    at,    the  second  by    a2,  ... 
and  add  the  results.     This  gives 


zl  2  «2  -f-  22  ^  «#  +  ...  sg  2  aq  -\-  ^  am 
-\-  x^a^  -\-  x^ab  -\-  .  .  .  xq^aq  =    2  a  A 


But  by  the  first  of     (B)     the  first  line  in  this  equation  is 
equal  to  zero,  and  therefore 
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xl  2  a2  +  x2  2  al  +  .  .  .  xq  2  aq  —  2  aA   =    0 
Also      «!  2  a*  -f  a-2  2  i2  +  .  .  .  xq  2  &?  —  2  *  A   =    0 

.........  (Q) 

«!  2  a?  +  a-,  2  fy  +  .  .  .  xq  2  ?2  —  2  ?A   =    0 

These  being  of  the  same  form  as  the  normal  equations 
(B),  the  value  of  x±  resulting  from  their  solution  will  be 
of  the  same  form  as  that  of  zl  in  the  solution  of  those  equa- 
tions, with  only  the  substitution  of  —  A  for  m.  From 
(G)  we  shall  therefore  have 

«!    —    ajAi    —    cujAjj    —    ...   On&n     =     0  (R) 

Multiply  the  first  of  (P)  by  v^  the  second  of  (P)  by 
i\j,  .  .  .  and  add  the  results.  Then 


+  .  .  .  (z 

=    2  Ay 


and  multiplying  the   first   of    (A)    by    a^    the  second   by 
a2,  .  .  .  and  adding  the  results 

Say   =   zl  2  a2  +  z2  2  ab  -f-  .  .  .  zq  2  aq  -|-  2  am 

=    0          by  the  first  of    (B). 
Also      2  fly    =  0 

2?y   =  0 
Substituting  these  values  in  the  above,  we  find  that 

2  my    =    2  AU  (S) 

Now   multiply  the   first   of    (A)    by    v»    the  second  by 
v     .  .  .  and  add  the  results.     This  shows  that 
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zx  2  av  +  z2  2  £y  +  .  .  .  zq  2  yu  +  2  mu    =    2  ua 
and  as  above, 

Say    =    2£y  ...    =    2?y    =    0 
Combining  this  result  with    (S)    we  have 

2  my    =    2y2    =    2  Ay  (T) 

Next  multiply   the    first  of    (A)    by    At,    the  second    by 
Aa,  .  .  .  ,  add  the  results  and  compare  with    (T).     This  gives 


zl  2  aA  -f-  22  2  #A  -f  ...  zg  2  ?A  -f- 

i=    2  Aw    =    2w2  (U) 

And  finally,  multiplying  the  first  of    (P)    by    A1?    the  second 
by   A2,  .  .  .  and  adding  the  results,  we  have 


zl 
+  a?!  S  «A  +  a?a  2  ^A  +  .  .  .  »92  S'A  =    2  A2 

Therefore  from    (U) 

2  v2  +  a?!  2  «A  +  a?2  2  *A  +  .  .  .  xq  2  ?A   =    2A2      (V) 


.  ... 

n  n 

We  must  now  find  the  mean  values  of  the  terms   X 
jc.j  2  #A,  ...  a*g  S  ^A.     Expanding    2  «A, 

2  aA  =  axAi  +  «2A2  +  .  .  .  a«An 
from  (R)  jcj  =  aiAt  +  a2A2  +  .  .  .  ^A^ 
Multiplying, 

«!  2  aA    =    a^Ai2  +  a2a2A22  +  .  .  .  attanA»f 
-[-    terms  in    AjA2,  AxAs,   .  .  . 
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In  the  most  probable  case  these  product  terms  vanish,  and 
substituting  for  A^,  A2a,  .  .  .  An2  the  mean  value  ^\  we 
have 


=     /x2  2  Oa 

=    M2       b 
Similarly  a;2  2  #A    =    /u2 


From    (W)    then,    q   being  the  number  of  unknowns, 

s  =  ^+«s 


=   ViT^  (»> 


V-z     = 


n 


=  V/i^r^         (») 


rz    =    .6745^    = 


Pz(n  —  q) 


83.    By  a  method  similar  to  that  used  in  paragraph  60,  we 
may  derive 

•y  „ 

a.d.     =  (7.3) 

—  q) 

Ztf 

(76) 


-  q) 

84.    Example.     In    illustration    of   the    above    processes 
take  the  example  in  paragraph  24,  where  we  found  for 


COMPUTATION  OF  THE  PRECISION  MEASURES.       73 

OBSERVATION     EQUATIONS 

2i   -  -   22  —   1.7    =    vl 

za  -   2.4    =    vt  ., 

-  Zl  +  z,  +  z8   _   1.0    =    v, 

22     --     2;{     --     3.0      =      V4 

and  for 

NORMAL    EQUATIONS 

2  zl  -  2  22  -        23  —   0.7    =    0  (a) 

-  2  2X  +  3  22  -  2.3    =    0  (b) 

zl  +  3  z3  --  0.4    ==    0  (c) 

SOLUTION  FOR  THE  WEIGHTS  BY  THE  FIRST  METHOD. 

For  finding  the  weight  of    2t    the  above  normal  equations 
would  be  written 

2  zl  —  2  22  -    28  —  1  =  0      (a') 
-  2  2i  +  3  22  =0      (b') 

-  zi       +82,      =0      (c') 

Solve  for  zl 

3  X  (a')     621-622-323-3  =  0 

(c')       zi       +  3  28      =  0 

5  zl  -  -  6  22       —3  =  0 

2  X  (V)   -  4  2l  +  6  22 =  0 

z,  -  3  =  0 

*i    =    3  At    =    y  (d') 

For  finding  the  weight  of    22    the  equations  are 

2  2l  _   2  22  -      23  =0  (a") 

—  2  2t  +  3  2,  —   1    =    0  (b") 

-  «,  +32.  =0  (c") 
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Solve  for  z, 

3  X  (a")  6  zt   -   6  z,  -  3  z,    =    0 

(c")  -      zt                 +  3  g»    =    0 


Substitute  in    (b")       za    =     -       .-.  jt>2,    =    -        (d") 

3  5 

For  finding  the  weight  of   z8   the  equations  are 

2^   --  2za  -        z4  =0  (a"') 

-  2  zx  +  3  z2  =0  (b'") 

-  z,  +  3  z8  -  1    =    0  (c'") 

Solve  for   z§ 

from  (b'")  2  sa    =    -  zx 

o 

substitute  in   (a"')  2  Zi  --   3  z8  =    0 

2  X  (c"')  -  2sx  +  6s8  --  2    =    0 

3z8  --  2    =    0 

«,  =  f     -     ^  =  |         (d-j 


SOLUTION  BY  THE  SECOND  METHOD. 

The  normal  equations  will  now  be  modified  so  as  to  appear 
in  the  following  form  :  — 

2z,  --  2zt  -        z,  --  0.7    =    A  (a) 

-  2zj  +  3  z,  -  2.3    =    B  (b) 

-      z,  +  3  z8  -  0.4    =    C  (c) 
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Solving, 

3  X  (a)          6  2i  --   6  2a  --   3  2,  --  2.1    =    3  A 
(c)_      _2i  _  +  32,  --  0.4    ==    C 
5  zl  -  6  22  --  2.5    =    3  ^  +  (7 
_  4  2t  4-  6  22  -  4.6    =    27?  _ 

«!  -7.1    =    3^  +  2^+tf     (8) 

21  =    7.1        and        pXl    =  (d) 

0 

Substituting  (S)   in   (b) 

3  22    =   16.5  -f6^  +  5J?  +  2(7 

3 

22  =    5.5        and        jt>^    =    -  (e) 

Substituting  (8)  in  (c) 

328    =    7.5  +  3JL  +  2.B  +  2<7 


28    =    2.5        and        pzs    =    -  (f) 


SOLUTION  BY  THE  THIRD  METHOD. 
The  normal  equations  are  now  taken  in  their  original  form. 

2  2t  -   2  22  —       28  —   0.7    =    0  (a) 

—  2  z1  -\-  3  22  —  2.3    =    0  (b) 

-      zl  +  3  2,  —  0.4    =    0  (c) 
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To  obtain   Zj   and  its  weight  we  proceed  as  follows  : 

from  (c)  z,    =    £   +    '- 

6  6 

from  (b)  za    =      ^    + 

o  o 

Substitute  in  (a) 

2%_*%_i?_*_:i_o.7«  o 

33  3  3 

Collecting  terms,  =    0 

3  3 

Zl    =    7.1         and         pZl    =    -  (d) 

o 

For   za 
3X  (a)  +  (c)          52l  -  6za  --  2.5    =    0 

6 

Zj    =    -  za  +    0.5 
o 

12 

Substitute  in  (b)       --  za  —   1.0  +  3  zt  —  2.3    =    0 
5 

o 

Collecting  terms,        -  z,  —  3.3    =    0 
5 

3 
za    =    5.5         and        />„    =    -  (e) 

For   z, 

3X(a)  +  2x(b)       2z,  --  3  z,  --  6.7    =    0 

3  6.7 


n  o  IT 

Substitute  in  (c)      ---  zt  —   -      +  3  z,  --  0.4    =    0 
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3  7.5 

Collecting  terms,  -  z,  —  -        =0 

2  2 

z8    =    2.5         and         pZi    =    |  (f) 

It  is  evident  that  the  three  methods  give  identically  the 
same  results  and  that  the  work  is  about  the  same  in  each  case. 
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Substituting   the  values   found   for    zly  z2  and  zs    in  the 
observations  equations  (A),  we  have 

7.1  _  5.5  —  1.7   =   vl   =    —  .1  .01  =  vf 

2.5  -  2.4   =   v2  =  +  .1  .01  =  u22 

-  7.1  +  5.5  +  2.5  -  1.0   =   va  —    •-  .1  .01  =  v32 

5.5  -  2.5  -  3.0   =   v4  =        .0  .00  =  v42 


.03   = 
In  this  example     n  =   4,     q   •=.   3. 


By  (72)  ,,  y  .17 

By  (74)  r    =    .6745;*  =    .12 


By  (73)       Mz,    =         Ll    =    .30  rw    =    .20 


5 

.03 

-       =    .U 

2 

By  (75)  </.^.    =    —    =    .16 

r    =    .8453  a.rf.    =    .13 


rza     =    .10 
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85.  Observations  of  Unequal  Weights.  If  the  observations 
are  not  all  of  the  same  weight  the  formulas  and  operations 
are  merely  modified  in  the  usual  manner  and  equations  (72) 
and  (75)  take  the  more  general  form 


(77) 


S  v  V7  p 
a.d.    =    — -  (78) 

V  n(n  —  q) 


CONDITIONED   OBSERVATIONS. 

86.  Suppose  there  are  given  n  observations,  n1  condi- 
tions and  q  unknown  quantities.  Then  by  paragraph  33, 
the  method  of  solution  is  to  eliminate  n'  unknowns  between 
the  "  Observation  "  and  "  Condition  "  equations,  leaving 
q  —  n'  independent  unknowns  in  the  "  Normal  "  equations. 
Consequently  formula  (77)  now  applies  to  this  case  and  it 
would  be  written 


=        __  (79) 

n  —  q  +  n' 


al80  ^   =    J~       *W* 

y  Pz(n  -  q  +  n') 


(80) 


The  weights  of  the  q  —  n'  unknown  quantities  can  be 
found  by  any  one  of  the  three  methods  already  given  and 
then  the  mean  error  of  each  unknown  may  be  computed  by 
using  formula  (80).  If  the  mean  errors  of  the  n'  quantities 
that  were  first  eliminated  is  wanted,  their  weights  must  be 
determined  by  eliminating  a  different  set  of  n1  quantities 
from  the  original  observation  equations  and  solving  the  neces- 
sary seta  of  equations  for  these  weights.  The  first  method  of 
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solution  for  the  weights  would  perhaps  be  best  here  as  the 
actual  values  of  the  unknowns  have  already  been  found. 
87.    Example.     Given  the 

OBSERVATION    EQUATIONS 

Zl  -}.  22  —   3.0  =  0  weight  1 

-  Z,+  1.5  =  0  «        4 

z2  -   2.2  =  0  "3 

Zl             +  28  —  3.4  =  0  "2 

and  the 

CONDITION    EQUATION 

zs  --  22  --  0.5    =    0  (B) 

To  find  the  most  probable  values  of   z1}  z2,    and   23,    and 
also  their  mean  errors. 

Eliminating   z3   between   (A)  and  (B)  there  remains 

Zl  _J_  22  _  3.0  =  0  weight  1 

z,  --   22  +  1.0  =  0  «       4 

z2  —  2.2  =  0  "3 

2l  -}-  z2  _  2.9  =  0  "2 

From  these  we  have  the 


NORMAL  EQUATIONS 

7  2l  _         22  _     4.8    —    0 

-  «»  -f-  10  «,  -  •   19.4    =    0 

Solving, 

2l    =    0.98  /?Z1     =    6.9 

z2    =    2.04  jo2l    =    9.9                     (D) 
from  (B)          z3    —    2.54 
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Now  eliminating   za   between  (A)  and  (B),  we  find 

Zl  _J-  Zg  _  3.5  =  0  weight  1 

2l  -   zs  +  1.4  =  0  «       4 

z8  -  2.7  =  0  "3 

2l  _L.  Z8  __  3.4  =  0  "2 

From  these  we  derive  the  new  set  of 

NORMAL    EQUATIONS 

7  zl  -         z3  -     4.3    =    0  m 

-  zt  +  10  z3  —  24.4    =    0 

and  solving  for   z8 

z8    =    2.54  pza    =    9.9  (G) 

Substituting  the  values  of   z^  z2   an(i   2s   in  equations  (A), 
we  have  the  residuals 

v  v2  pv* 

.02  .0004  .0004 

.06  .0036  .0144 

.16  .0256  .0768 

.12  .0144  .0288 

.1204    =    2jw2 

Since  in  this  example    n    =    4,     n'    =    1,     £    =    3, 


By  (79) 


=  -08 
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The  first  significant  figure  in  the  mean  errors  is  so  large 
that  it  is  not  worth  while  to  retain  the  second  place  as  usual. 

If  desired,  the  probable  errors  and  average  deviations  can 
now  be  computed  by  the  usual  formulas. 

88.  In  case  the  observations  are  made  directly  upon  the 
values  of  several  quantities  subject  to  certain  conditions,  we 
have  n  =  £,  and  equation  (79)  reduces  to 


v^ 

from  which  the  mean  error  of  any  observation  may  at  once 
be  computed  from  its  weight. 

89.    Example.     Taking  the  example  in  paragraph  34  on 
the  measurement  of  the  angles  of  a  quadrilateral,  we  had  for 

OBSERVATION    EQUATIONS 

z1    =    0  weight  3 

z2    =    0  "2 

z3    =    0  "2 

24      =      0  «  1 

for  the 

CONDITION    EQUATION 

«i  +  z2  +  z3  +  z*  +  58    =    0  (B) 

and  for  the 

NORMAL    EQUATIONS 

4  21  +      z2  +      28  +  58    =    0 
z,  +  3  z2  +      z3  +  58    =    0  (C) 

2l  _j_      22  +  3  28  +  58    =    0 

Solving  these  equations  for  the  values  of  the  unknown 
quantities  and  also  for  their  weights,  we  have 

Zi    —  8.3  pzi    =    3.5 

z2    =          12.4  pst    =    2.5 

23    =    -   12.4  pZi    =    2.5 
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and  forming  a  new  set  of  normals  containing   24,    we  find,  on 
solving  for  that  quantity, 

z4    =      -  24.9  Pz4    =    1.75 

These  values  of  2,,  z2,  z8,  z4  are  also  in  this  case  the 
residuals  of  the  observations,  and  therefore  to  compute  the 
precision  measures  we  have 

v  v2  pv* 

8.3  68.9  206.7 

12.4  153.8  307.6 

12.4  153.8  307.6 

24.9  620.0  620.0 

1441.9    = 
n'    =    1 


By  (81)  ,*  ==   y    "x       =   38 

=   20  rzi  =   13 

=   24  r^  =    rZ9   =    16 


-  =        =  29      r^  =  20 

We  shall  accordingly  write  for  the  most  probable  values  of 
the  angles  of  the  quadrilateral 


A 

— 

101° 

13' 

14" 

± 

13" 

B 

•  — 

93 

49 

5 

± 

16 

C 

SB 

87 

5 

27 

± 

16 

D 

•=. 

77 

52 

15 

j_ 

20 

In  the  original  solution  in  paragraph  34  it  is  evident  that 
the  results  were  carried  out  to  a  greater  number  of  places  of 
hi<_milicant  figures  than  the  character  of  the  observations 
warranted. 


CHAPTER  V. 

MISCELLANEOUS   THEOREMS. 

THE   DISTRIBUTION   OF    ERRORS. 

90.  Having  developed  the  processes  for  the  adjustment  of 
observations  according  to  the  Method  of  Least  Squares,  it 
will  now  be  interesting  to  show  how  closely  the  distribution 
of  errors  found  in  actual  practice  corresponds  to  the  theo- 
retical distribution  upon  which  our  methods  of  solution  are 
based. 

By  formula  (36),  the  probability  that  the  error  of  a  single 
observation  will  be  numerically  less  than  a  is 

P   =    ~  fV^Vx  (82) 

VTT  Jo 

Let    t    =    h  x,    .«.    dt    =    h  dx.     Also  when    x    =    0, 

t    =    0 ;    and  when    x    —    a,     t    =    ha    =    p  ®.      Sub- 

r 
stituting  in   (82), 

P   =    4=  Cp'e-*dt  (83J 


Values  of   P   for  values  of  the  argument   —   are  given  in 

Table  I.  Also  for  any  series  of  observations  this  quantity 
P  will  represent  the  fraction  of  the  entire  number  which 
should  have  errors  less  than  the  amount  a.  Hence  if  P  is 
multiplied  by  the  whole  number  of  observations  the  result 
will  be  the  number  of  errors  which  should  be  less  than  the 
limit  a. 


84 


METHOD  OF  LEAST  SQUARES. 


91.  Example.  Forty  measurements  on  the  diameter  of 
Saturn's  ring  were  made  by  Bessel,  with  the  following 
results :  — 


M    v 

M     v 

M    v 

M    v 

38".91  -.40 

39".35  +.04 

39".41  +.10 

39".02  -.29 

39  .32  +.01 

39  .25  —  .06 

39  .40  +.09 

39  .01  -  .30 

38  .93  -  .38 

39  .14  —  .17 

39  .36  +  .05 

38  .86  -  .45 

39  .31    .00 

39  .47  +.16 

39  .20  -  .11 

39  .51  +.20 

39  .17  -  .14 

39  .29  -  .02 

39  .42  +.11 

39  .21  -  .10 

39  .04  —  .27 

39  .32  +.01 

39  .30  -  .01 

39  .17  -  .14 

39  .57  +.26 

39  .40  +.09 

39  .41  +.10 

39  .60  +.29 

39  .46  +-15 

39  .33  +.02 

39  .43  +.12 

39  .54  +-23 

39  .30  -  .01 

39  .28  -  .03 

39  .43  +-12 

39  .45  +.14 

39  .03  -  .28 

39  .62  +.31 

39  .36  +.05 

39  .72  +.41 

From  these  the  most  probable  value  of  the  diameter  is 
found  to  be 

D    =    39".308  ±  0".022 

the  probable  error  of  a  single  observation  being    r  =  0".136. 
Compare  the  theoretical  and  actual  distribution  of  errors 

between  0".00  and  0".05 

"  0  .05  "  0  .10 

"  0  .10  "  0  .20 

"  0  .20  "  0  .30 

"  0  .30  "  0  .40 

over  0  .40 

In  the  following  table  the  first  column  gives  the  successive 
values  of  the  limiting  error    a,    the  second  column  the  values 

of        ,    and  the  third  column  the*  corrt'sponiliii.^  values  of    1\ 
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The  fourth  column  contains  the  differences  between  the  suc- 
cessive values  of  P,  and  by  multiplying  each  of  these 
differences  by  40,  the  number  of  observations,  we  have  the 
quantities  in  column  five,  which  are  the  numbers  of  errors 
that  according  to  the  theory  should  fall  within  the  corre- 
sponding limits.  Column  six  shows  the  actual  number  of 
residuals  occurring  between  these  limits. 


a 

a 
r 

P 

d 

n 

n1 

0.00 
0.05 
0.10 
0.20 
0.30 

0.000 
0.368 
0.735 
1.471 
2.206 

0.000 
0.196 
0.380 
0.679 
0.863 

0.196 
0.184 
0.299 
0.184 
0.090 

8 
7 
12 
7 
4 

9 
6 
12 

8 
3 

0.40 

2.942 

0.953 

0.047 

2 

2 

00 

00 

1.000 

This  is  a  close  agreement  considering  that  the  number  of 
observations  is  not  very  large.  Also  the  number  of  errors 
greater  than  the  probable  error  should  be  equal  to  the  num- 
ber less  than  it.  On  counting  the  residuals  we  find  twenty-one 
less  than  0".136  and  nineteen  greater. 


THE    REJECTION    OF   OBSERVATIONS. 

92.  After  a  series  of  measurements  have  been  made,  it  is 
frequently  found  that  one  or  two  of  the  observations  differ 
widely  from  the  others,  and  hence  it  becomes  a  matter  of 
great  importance  to  establish,  if  possible,  some  criterion  by 
which  we  may  determine  whether  such  discordant  observa- 
tions should  be  rejected  or  not.  We  are  not  concerned  here 
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with  the  question  of  the  detection  of  a  mistake  or  constant 
error,  which  a  consideration  of  the  circumstances  of  the 
observations  or  of  the  instruments  might  reveal,  but  it  is 
assumed  that  there  is  nothing  whatever  to  guide  us  except 
the  mere  fact  of  the  unusual  size  of  the  residuals  of  the 
observations  under  discussion.  To  reject  an  observation 
merely  because  it  differs  considerably  from  the  others  is 
entirely  unjustifiable,  while  to  retain  it  without  any  investi- 
gation is  a  neglect  of  the  evidence  furnished  by  the  observa- 
tions themselves. 

The  adoption  of  any  rigid  criterion  based  upon  the  magni- 
tude of  the  residuals  is  perhaps  more  satisfactory  from  a 
mathematical  standpoint  than  from  that  of  a  practical  observer, 
and  some  of  the  latter  are  of  the  opinion  that  no  observation 
should  be  rejected  entirely,  even  the  most  widely  discordant 
ones  being  given  a  certain  weight.  In  this  latter  case, 
however,  the  Theory  of  Probability  will  furnish  a  guide  as  to 
the  proper  weights  to  assign  to  the  different  observations. 

Of  the  various  criteria  that  have  been  proposed,  that 
developed  by  Pierce  (see  Chauvenet,  page  558)  is  perhaps 
the  most  complete.  The  derivation  and  application  of  this 
criterion  is,  however,  somewhat  long  and  complicated,  and 
for  all  ordinary  cases  the  following  simple  methods  will  give 
practically  as  good  results. 

93.  Criterion  for  the  Rejection  of  a  Single  Doubtful 
Observation.  It  was  shown  in  (83)  that  in  a  series  of  n 
observations  the  number  of  errors  numerically  less  than  a 
should  be  nP,  and  therefore  the  number  of  errors  greater 
than  a  should  be 

n  -  nP    =    n(l  -  P)  (84) 

If  the  value  of  the  expression  in  (84)  is  less  than  one-half, 
the  occurrence  of  an  error  of  magnitude  a  will  have  a 
greater  probability  against  it  than  for  it,  and  hence  the  obser- 
vation corresponding  may  be  rejected. 
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Accordingly  the  limit  of  rejection,    a,    of  a  single  doubtful 
observation  is  obtained  from  the  equation 


*  i  -      =  - 


or 


P    = 


2  t&  -   I 

In 


(85) 


94.  Example. 
angle  are  made, 
rejected  ? 


Fifteen  observations  on  the  value  of  an 
Ought   any   of    the   observations   to   be 


M 

V 

V2 

v' 

v'2 

v" 

v"2 

2°     23'.90 

—    .30 

.090 

-  .41 

.168 

-  .33 

.109 

23.76 

-    .44 

.194 

-.55 

.303 

-  .47 

.221 

25.21 

+  1.01 

1.020 

+  .90 

.810 

24.68 

+    .48 

.230 

+  .37 

.137 

+  .45 

.203 

23.96 

-    .24 

.058 

-  .35 

.123 

-.27 

.073 

24.26 

+    .06 

.004 

-  .05 

.003 

+  .03 

.001 

24.83 

+    .63 

.397 

+  .52 

.270 

+  .60 

.360 

24.07 

.13 

.017 

-.24 

.058 

-  .16 

.026 

23.98 

-    .22 

.048 

-  .33 

.109 

—  .25 

.063 

24.14 

-    .06 

.004 

-.17 

.029 

-.09 

.008 

24.40 

+    .20 

.040 

+  .09 

.008 

+  .17 

.029 

24.38 

+    .18 

.032 

+  .07 

.005 

+  .15 

.023 

24.59 

+    .39 

.152 

-  .28 

.078 

+  .36 

.130 

24.10 

-    .10 

.010 

+  .21 

.044 

-.13 

.017 

22.80 

-1.40 

1.960 

2°     24  .20 

4.256 

2.145 

1.263 

M0 

2ua 

2i>'2 

2y"2 
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Using  all  the  observations  we  find 


M9    =    2°    24'.20         r    =    .6745  t/4-256    =    .37 

14 

90 

By  (85),  P    =    —    =    .967 

30 

-By  Table  I,          -    =    3.17      .-.      a    =    1.17 

As  the  residual  —  1.40  is  larger  than  a,  we  reject  the 
last  observation. 

From  the  remaining  observations  we  now  compute  a  new 
mean  value  and  a  new  set  of  residuals.  And  we  find 


'0    =    2°  24'.31         r'    =    .6745  V  —     -    =    .27 

V     13 


By    (85),  P   =  =    .964 

By  Table  I,          ^    =    3.11      .-.      a    =    .84 

The  third  observation  may  accordingly  be  rejected. 
From  the  thirteen  observations  that  remain  we  find 


M  "Q    =    2°  24'.23         r"    =    .6745  y  i^.  =    .22 


P    =  =    .962 

26 


-    =  3.08      .-.      a    =    .68 

r 


Therefore  no  more  observations  are  to  be  rejected. 
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95.  The  Huge  Error.  In  cases  where  the  number  of 
observations  is  not  unusually  large,  a  simple  and  safe  criterion 
for  the  rejection  of  a  doubtful  observation  is  found  in  tin  DM 
of  the  "  Huge  Error." 

This  is  an  error  of  such  a  magnitude  that  999  out  of  every 
1000  errors  are  less  than  it  and  only  1  as  large  as  or  greater 
than  it. 

Therefore  the  probability  that  the  error  of  any  given 
observation  will  be  less  than  the  "  Huge  Error  "  is  .999,  and 
from  Table  I, 

when  P    —    .999,  —    =    4.9 

a   =    Huge  Error   =    4.9  r 

=    3,3  |i  (86) 

=    4.1  a.d. 

Then  in  any  limited  series  of  observations,  if  an  error 
greater  than  the  huge  error  is  found,  we  should  reject  the 
observation  corresponding. 

See  also,  Holman,  page  30 ;    Wright,  page  131. 


CONSTANT    ERRORS. 

96.  Throughout  our  discussion  of  the  methods  of  adjusting 
observations  so  as  to  obtain  from  them  the  most  probable 
values  of  the  unknown  quantities,  all  constant  errors  are 
supposed  to  have  been  eliminated  before  the  Method  of 
Least  Squares  is  applied  in  deducing  the  results. 

If  this  is  not  done,  and  each  observation  is  subject  to  the 
same  constant  error,  the  final  result  will  be  affected  by  an 
equal  amount,  and  in  short,  the  Method  of  Least  Squares  is 
not  capable  of  removing  or  reducing  the  effect  of  errors  of 
this  kind.  All  that  is  accomplished  l.y  the  use  of  the  nu-thod 
is  to  reduce  to  a  minimum  the  effect  of  the  Accidental  Errors. 
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Hence  it  will  be  seen  that  although  by  increasing  the 
number  of  observations  of  a  given  kind  we  may  increase  the 
precision,  that  is,  reduce  the  probable  error,  of  our  final 
result  as  much  as  we  choose,  yet  we  do  not  in  this  way 
necessarily  increase  the  accuracy  of  the  determination. 

But  if  the  unknowns  can  be  determined  in  several  ways,  or 
under  a  variety  of  different  circumstances,  with  various 
instruments,  or  by  different  observers,  then  it  is  most  proba- 
ble that  the  constant  errors  of  the  different  sets  of  measure- 
ments will  be  grouped  about  the  true  values  of  the  unknowns 
according  to  the  exponential  law  of  error.  Accordingly  a 
combination  of  such  observations  will  enable  us  to  increase 
not  only  the  precision,  but  also  the  accuracy  of  the  final 
result,  the  constant  errors  of  the  different  sets  tending  to 
cancel  each  other  in  the  same  way  that  the  accidental  errors 
of  a  single  .set  do. 

It  is  for  this  reason  that  determinations  of  a  quantity  from 
observations  made  in  a  variety  of  ways  are  more  valuable 
than  those  obtained  merely  from  different  sets  of  measure- 
ments of  the  same  kind. 

97.  The  probability  of  the  existence  of  a  constant  error 
may  often  be  expressed  in  the  following  manner. 

Example.  A  standard  100  ohm  coil  is  compared  with 
a  Wheatstone's  bridge  and  the  mean  result  found  to  be 
100.90  ±  0.20.  To  find  the  probability  that  there  is  an  error 
in  the  bridge  between  -(-  0.80  and  -j-  1.50  ohms. 

Suppose  the  result  100.90  ±  0.20  is  treated  as  a  single 
observation,  and  we  find  by  an  application  of  (83)  the  prob- 
ability that  the  error  of  this  observation  is  numerically  less 
than  0.60  ohms. 

Here  -    =    —    =    3.00      .-.      P    =    .957 

r  .20 

Hence,  as  far  as  is  shown  by  the  observations,  the  proba- 
bility that  100.90  ohms  is  within  0.60  ohms  of  the  true 
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value  is  .957.  But  since  it  is  known  that  the  true  resistance 
i>  I  00  ohms,  it  follows  that  there  is  the  same  probability 
th:it  there  is  a  constant  error  in  the  bridge  between  -|-  0.30 
and  -f-  1.50  ohms. 

98.  Combination  of  Determinations  having  Different  Con- 
stant Errors.  In  case  two  or  more  determinations  of  a 
quantity,  together  with  their  probable  errors,  are  obtained, 
tlu-  method  of  combining  them  so  as  to  secure  the  best  final 
result  was  considered  in  paragraph  56,  and  in  Example  C, 
rraph  57.  But  it  was  there  assumed  that  all  the  results 
were  subject  to  the  same  constant  errors,  while  if  this  is  not 
true  the  probable  errors  of  the  separate  determinations  bear 
no  relation  to  their  weights,  and  accordingly  in  such  cases 
another  process  must  be  adopted. 

To  determine  whether  the  different  measurements  may 
fairly  be  considered  to  have  the  same  constant  errors  we  may 
proceed  as  follows  :  — 

Let  the  determinations  of  the  quantity   M  be 

MI  ±  n  (a) 

JtT.  ±  rz  (b) 

and  let  the  difference  between  these  results  be 

d    =    J[/i  —  MI  (c) 

Then  the  probable  error  of    d    is  by  (66), 


E    =    yV  +  ra»  (d) 

If  d  is  of  such  a  magnitude  that  an  accidental  error  as 
great  as  it  may  reasonably  be  expected,  we  may  assume  that 
the  constant  errors  of  (a)  and  (b)  are  the  same,  and  pro- 
ceed as  in  Example  C,  paragraph  57. 

But  if  the  probability  of  making  two  determinations  which 
differ  by  the  amount  d  is  very  small,  we  had  best  consider 
J/!  and  MI  to  have  the  same  weight,  provided  there  are  no 


92  METHOD  OF  LEAST  SQUARES. 

special  reasons  for  regarding  one  better  than  the  other.  The 
final  value  ot  M  will  then  be  the  arithmetical  mean  of  MI 
and  MS,  and  its  probable  error  will  be  found  by  (53). 

99.    Example  A.     An  angle  is  measured  by  a  theodolite 
and  by  a  transit  with  results 

By  Theodolite,     24°    13'   36".0    ±      3".l 
By  Transit,          24°    13'   24"      ±    14" 

What  is  the  most  probable  value  of  the  angle  and  its  prob- 
able error? 

Referring  to  the  preceding  paragraph, 

rt    =    3.1,         rz    —    14,         d    =    12, 
and  the  probable  error  of   d  is 


R    —   v/  3.12  +   142    =    14 

Then  from  Table  I  the  probability  that  the  accidental 
error  of  a  determination  will  be  at  least  as  large  as  12  is 
found  from 

«    =    i?    =    .86        ,.         1  -  P    =    .57 
r  14 

That  is,  there  is  more  than  an  even  chance  that  two  such 
determinations  of  the  angle  will  differ  by  as  much  as  12. 

Hence  it  is  fair  to  assume  that  the  two  determinations  are 
not  affected  by  constant  errors  of  different  magnitudes,  and 
they  would  be  combined  as  in  Example  C,  paragraph  57. 

Example  B.  Suppose  the  zenith  distance,  JtfJ  of  a  star, 
observed  at  two  different  culminations,  is  found  to  be 

J/J    =    14°    53'   12".10    ±   0".30 
j*fa    =    14°    53'   14".30    ±   0".50 

What  is  the  best  final  value  ? 


Here          d   =    2.2,         It    =    v^9  +  .25    =    .58 

and  for       -    =    —    =    3.8,         1   -  P    =    .01 
r  .58 
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Therefore  the  chance  that  tin*  diffrn-nre  in  the  two  deter- 
minations, due  to  arridental  errors,  will  be  as  large  as  2.2  is 
only  one  in  a  hundred.  It  is  to  be  concluded  then  that  the 
constant  errors  of  observations  at  the  two  culminations  differ 
by  about  2.2,  and  as  there  is  nothing  to  show  that  one  meas- 
urement is  more  accurate  than  the  other  we  will  give  them 
both  the  same  weight  and  take  the  mean.  Then  the  best 
value  for  the  zenith  distance  is 

M0    =    14°   53'    13".20    ±    0".74 

For  JU    =    14°   53'  +    18>'10  +  14*M 

=    14°   53'   13".20 

By  (53)        ra  .6745  \J  ^^    =    .74 

For  a  more  extended  treatment  of  this  subject  see  Johnson, 
"  The  Theory  of  Errors  and  Method  of  Least  Squares," 
chap.  vii. 

THE   WEIGHTING   OF    OBSERVATIONS. 

100.  In  case  the  relative  worth  of  observations  is  not 
settled  by  methods  already  discussed,  the  proper  weight  to 
assign  to  each  quantity  in  the  final  adjustment  can  only  be 
determined  from  a  full  knowledge  of  all  the  circumstances  of 
the  measurements.  Even  then  considerable  experience  in 
the  particular  work  in  hand  is  required  before  the  best  values 
for  these  weights  can  be  assigned.  The  weight  given  to  a 
quantity  should  never  be  considered  final,  but  always  subject 
to  revision  whenever  new  information  with  regard  to  the 
quantity  is  obtained.  Thus  an  observation  which  at  first  is 
supposed  to  deserve  a  high  degree  of  confidence  is  often 
found  on  later  investigation  to  possess  very  little  value,  and 
vice  versa. 

See  also  Wright,  page  118. 
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OTHER    LAWS   OF    ERROR. 

101.  Although  in  the  great  majority  of  cases  the  distribu- 
tion of  errors  follows  the  exponential  law  thus  far  considered, 
there  are  a  few  special  cases  in  which  some  of  the  suppositions 
made  in  deriving  that  law  do  not  hold,  and  hence  for  the 
adjustment  of  such  observations  the  corresponding  special 
laws  of  error  must  be  determined. 

.  For  instance,  in  applying  the  exponential  law  we  assume  a 
large  number  of  observations,  that  each  observation  is  sub- 
ject to  the  same  law  of  error,  that  small  errors  are  more 
likely  to  occur  than  large  ones,  and  that  positive  and  negative 
errors  are  equally  probable.  Now  it  is  easy  to  conceive  of 
cases  where  only  positive  errors  can  occur,  or  where  the 
probability  of  the  occurrence  of  a  small  error  may  not  be 
greater  than  that  of  a  larger  one,  etc.  If  we  can  determine 
the  different  sources  of  error  in  any  case  and  the  relative 
effect  of  each  upon  the  quantity  sought,  we  shall  arrive  at 
the  law  of  error  for  that  particular  set  of  observations.  The 
case  of  most  common  occurrence  is  the  following. 

102.  Suppose  all  errors  between  the  limits    a   and    —  a 
are  equally  probable,  and  that  there  are  no  errors  beyond 
these  limits.      Then  if    y  =   <t>(%)    is  the  equation  of  the 
Curve  of  Error,  and  its  area   is   represented   as   in   para- 
graph 18,  we  have 

J%(«0  &    =    *  (a) 

or  2  4>(x)  Cadx    =    1  (b) 

Jo 

since  by  the  supposition  made   <£(x)    must  be  a  constant. 
Integrating  and  solving  for   <£(»),  we  have 


y    =     <|>(x)     =     --  (87) 

2  a 
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To  find  the  Mean  Error  we  have  l>y  <k-iinition  as  in  para- 
graph 52 

P*  =  r  x^(x)  dx 

J  -a 
\     f*a 

a  Jo 


n      


a 

7? 


The  Probable  Error  is  derived  from  the  equation 

1    Cr  1 

a  Jo  2 

T 

Finally,  for  the  Average  Deviation  we  have 

Ca 
a.d.    =     !     x  <j>(x)  ax 

J -a 
I     Ca        j 

=    -   I     x  dx 
a*** 

j  a 

7 

And  the  Curve  of  Error  has  the  form 


(88) 


(89) 


(90) 
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That  the  average  deviation  and  probable  error  are  in  this 
case  equal  to  one  half  of  a  may  also  be  seen  from  the  defini- 
tion of  these  quantities. 

Example.  In  taking  a  logarithm  from  a  four  place  table, 
what  is  the  probable  error  of  the  mantissa  ? 

In  this  case  the  maximum  error  is  .00005,  and  all 
errors  between  —  .00005  and  .00005  are  equally  probable. 
Therefore 

r    =    .000025 

103.  The  only  other  special  case  of  common  occurrence  is 
that  in  which  the  error  of  a  quantity  is  due  to  two  sources, 
each  of  which  can  with  the  same  probability  assume  all 
values  between  a  and  —  a.  Here  it  may  be  shown  that 
the  curve  of  error  consists  of  two  straight  lines  whose  equa- 
tions are 

2a  —  oc  2a  +  oc 

y      ~  y      ~ 


Also      |i"   =    -    «%  r   =    (2  -  \/2")  a      (92) 

I 

For  a  more  extended  discussion  of  special  laws  of  error,  see 
Wright,  paragraphs  31  to  39.  See  also  Example  151. 

104.  Contradictory  Observations.  Suppose  three  observa- 
tions of  a  quantity  give  results  55,  56,  91.  It  is  obvious  that 
to  take  the  arithmetical  mean  as  the  most  probable  value  in 
this  case  would  be  contrary  to  the  evidence  furnished  by  the 
measurements,  while  in  so  small  a  number  of  observations  it 
would  not  be  allowable  to  reject  the  value  91  entirely.  With 
such  a  series  of  observations  no  satisfactory  solution  can  be 
obtained,  but  probably  the  best  thing  to  do  would  be  to  take 
the  value  which  has  as  many  observations  less  than  it  as  it 
has  greater,  or  56. 


CHAPTER  VI. 

GAUSS'S   METHOD    OF   SUBSTITUTION. 

105.  The  most  laborious  part  of  the  application  of  the 
Method  of  Least  Squares  to  the  adjustment  of  observations 
consists  in  the  formation  and  solution  of  the  "Normal  Equa- 
tions," and  this  labor  increases  enormously  with  increase  in 
the  number  of  observations,  of  unknowns  and  of  conditions. 
It  is  not  at  all  unusual  to  find  that  the  adjustment  of  a  single 
set  of  observations  takes  several  weeks,  even  with  all  the  aid 
that  can  be  obtained  from  tables  of  logarithms,  of  squares,  of 
products,  and  of  reciprocals,  and  also  from  the  use  of  calcu- 
lating machines.     In  general  the  computations  can  be  per- 
formed more  rapidly  and  with  less  fatigue  by  using  a  machine 
or  a  table  of  products  than  by  using  logarithms,  but  a  combi- 
nation of  methods  is  often  desirable. 

106.  In  any  case  however  it  is  of  the  utmost  importance 
that  the  formation  and  solution   of   the   normal   equations 
should  be  effected  in  a  systematic  way,  and  that  as  far  as 
possible  checks  on  the  numerical  work  be  carried  along  in 
the  computations.     By  a  slight  amount  of  additional  work 
checks  upon  the  results  at  successive  stages  of  the  solution 
may  be  obtained  by  means  indicated  in  the  following  dem- 
onstrations, while  the  most  satisfactory  form  for  the  solution 
of  the  normal  equations  is  given  by  the  "Method  of  Substitu- 
tion" proposed  by  Gauss. 

This  method  will  now  be  explained ;  and  it  is  to  be  observed 
that  it  is  customary  in  this  subject  to  enclose  a  quantity  in 
brackets  when  the  sum  of  a  number  of  quantities  of  the  same 
kind  is  to  be  denoted.  Thus  [##]  means  the  same  as  ^ab. 

For  the  sake  of  simplicity  in  demonstration  it  will  be 
assumed  that  the  observations  are  all  reduced  to  weight  unity. 
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107:    Checks  on  the  Formation  of  the  Normal  Equations. 
If,  as  in  paragraph  76,  we  take  for 

OBSERVATION    EQUATIONS 


(A) 


we  shall  have  for 

NORMAL   EQUATIONS 

[aa]  «!  +  [aft]  22  +  .  .  .  [a?]  2^  +  [am]    =    0 
3l  +  [ftft]  *»+...  [ft^]  2,  +  [ftro]    =    0 

.........  (B) 


2  +  .  .  .   [qq]  zq  +  [?m]    =    0 

Let 

«i  +  *i  +  •  •  •  Si  +  mi    =    si 


(C) 


•••     [«]  +  [*]  +  •••[?]  +  [»»]    =    [•]         (D) 

Multiplying  the  first  of   (C)   by  ml9  the  second  by  m,,  .  .  . 
and  adding,  there  results 

[am]  +  [6m]  +  .  .  .  [gm]  +  [mm]  =  [«m]     (93) 

Next  multiplying  each   of   equations    (Cj   by  its   a   and 
adding,  and  then  each  by  its   ft   and  adding,  etc.,  we  have 
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[am~\  = 


[a6]  +  [66]  +  .  .  .  [60] 


Equation  (93)  will  be  satisfied  if  the  absolute  terras  in 
the  normal  equations  are  correct,  and  equations  (94)  when 
the  coefficients  of  the  unknown  quantities  are  correct.  These 
check  the  formation  of  the  normal  equations. 

108.  The  Reduced  Normal  Equations  and  the  Elimination 
Equations. 

The  value  of  zt  in  terms  of  the  remaining  unknowns, 
derived  from  the  first  of  equations  (B),  is 

[aft]  [ac]  [am] 

'   [aa]*2  "   [aa]*8  "          '   [^aj 

Substituting  this  in  the  remaining  n  —  1  equations,  they 
become 


And  letting 


,  1]    =    [M]   - 


[oa] 
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the  above  equations  take  the  following  form,  which,  being 
the  same  as  that  of  the  original  normal  equations,  they  are 
called  the 

FIRST   REDUCED   NORMAL   EQUATIONS. 

[W,  1]  2,  +  [ftc,  1]  z8  +  .  .  .  [fy,  1]  zq  +  \bm,  1]    =    0 

&clz  cclz         ...clz  Cml     =    0 


,  1]  z2  +  [c0,  1J  f.  +  .  .  .  [??,  1]  s,  +  [0m,  1]  =   0 

An  inspection  of  equations  (F)  will  render  it  easy  to  form 
a  rule  for  writing  out  any  one  of  them. 

Now  by  means  of  the  first  of  equations  (G),  eliminating 
z2  from  each  of  the  others  in  the  same  way  that  zl  was 
eliminated  from  the  normal  equations,  there  results  the 

SECOND  REDUCED  NORMAL  EQUATIONS 
[cc,  2]  z3  +  .  .  .  [c0,  2]  zq  +  [cm,  2]    =    0 


[eg,  2]  za  -f-  .  .  .  [gg,  2]  zq  -f-  [gra,  2]    =    0 
In  which 

[cc,  2J     =     [cc»  lj]    —     — T-/l        — = 

(i) 


Continuing  this  process  we  shall  finally  arrive  at  the  single 
equation 


from  which  the  value  of   zq   is  determined. 
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The  value  of  zq_l  will  then  be  obtained  by  substituting  the 
numerical  value  of  zq  in  the  first  of  the  preceding  >H  of 
equations,  and  so  on,  until  finally  zl  is  obtained  from  the 
first  of  the  original  normal  equations.  The  equations  from 
which  the  unknowns  are  actually  determined  are  then  the 
following,  called  the 

ELIMINATION   EQUATIONS. 

Itr  ft]  Zi  +  \_al>~\  «a  +  .  .  .      lag]  zq  +  [am]  =  0 
[bb,  1]  *,+...  [bq,  \]  zq  +  [bm,  1]  =  0 

.........  (9.5) 


It  may  be  seen  from  the  rule  in  paragraph  81  that 
\_qq,  q  —  1]  is  the  weight  of  zg,  and  the  weight  of  any 
unknown  might  be  found  at  the  same  time  as  its  value  by 
making  it  the  last  in  the  order  of  elimination,  but  except  in 
special  cases  the  weights  had  best  be  obtained  by  the  general 
process  of  paragraph  115. 

109.  Check  on  the  Solution  of  the  Normal  Equations. 
Multiplying  the  first  of  the  Observation  Equations  (A)  l»y 
mt,  the  second  by  ra2,  .  .  .  and  adding  the  results,  we  have 

[mo]     =     [am]  Zj_  -|-   [bm]  Z2  +  .  .  .  \_qrn]  zq  +   [nun] 

But  in  equation  (T),  paragraph  82,  it  was  shown  that 
[ray]  =  [vv].  Therefore 

[vv]    =    [am]  Zi  -\-  [bm]  z2  -f  .  .  .  [qm]  zq  +  [mm] 

Substituting  in  this  the  value  of  zl  from  the  first  of  (95), 
we  get  the  result 

[vv]    =    [bm,  1]  ga  +  [cm,  1]  z3  +  .  .  .  [,ym,  1]  zq  +  [m;//,  1] 
in  which 
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\_ab~]  [am] 

(K) 


[6m,  1]    = 


_  fam"]  [am] 

[mm,  1]    =  [mm]  -  ^ 

being  similar  in  form  to  equations  (F). 

Next  eliminating   22    in  a  like  manner,  we  get 

[wt>]    =    [cm,  2]  28  +  .  .  .  [?m,  2]  2<?  +  [mm,  2] 
and  continuing  this  process  it  finally  appears  that 

[w]    =    [mm,q]  [96] 

110.  Arrangement  of  the  Computations.  In  computing 
the  coefficients  that  appear  in  the  "Auxiliary"  or  "Reduced 
Normal  Equations"  it  is  most  convenient  to  arrange  the  work 
in  tabular  form.  The  arrangement  of  the  solution  will  be 
illustrated  for  an  example  containing  four  unknowns  but  it 
will  be  evident  that  the  process  can  be  extended  to  cover  any 
case. 

Let  the  Observation  and  Normal  Equations  be  represented 
by  equations  (A)  and  (B),  there  being  only  four  unknowns 
2,,  22,  28,  24,  and  arrange  a  table  as  on  the  next  page.  In  this 
scheme  of  solution  the  upper  lines  of  the  rows  in  the  first 
compartment  contain  all  the  quantities  that  appear  in  the 
Normal  Equations,  together  with  [mm]  and  the  quantities 
in  the  column  headed  8  which  are  used  in  checking  the 
results  in  accordance  with  equations  (93)  and  (94).  The 
other  compartments  contain  the  corresponding  quantities 
for  the  Reduced  Normal  Equations,  and  the  first  line  in 
each  compartment  gives  the  coefficients  in  the  Elimination 
Equations. 
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SCHEME  A. — SOLUTION  OF  THE  NORMAL  EQUATIONS. 


M 

log[aa] 

M 

log  [aft] 

M 

log  [or] 

[arf] 

log[arf] 

[am] 
log  [am] 

M 

log[«] 

log^> 
log  4: 
log  ^4 
log  -4. 

m 

At,  [«*] 
log  ^6  [aft] 

M 

<<UM 

log^[ae] 

[M] 

-^•[«rf] 

log^6[arf] 

M 

J6[am] 
log  ^4  [am] 

M' 

A*[?u} 
Iog^fc[a«] 

[00] 

4.M 

logJe[«e] 

[erf] 

4.C«rf] 

log^e[arf] 

w 

J.[«m] 
log  Jc  [am] 

M 
-4.M 

logJc[«] 

[«W] 

^[arf] 
logJ^arf] 

[rfm] 
^d  [am] 
log  -^rf  [<"»] 

[A] 

4-M 

log  ^rf  [a*] 

[mm] 
Am  [am] 
log  ^m  [am] 

M 
J.M 

log^m[«] 

log  B. 
log  3, 
\oSSn 

[66.  1] 
Jog  [66,1] 

[te,  1] 
log[6e,l] 

[1*1] 

log[6rf,  1] 

[6m,  1] 
log  [6m,  1] 

[6,.  1] 
log  [ft.,  1] 

[ee.l] 

*.[*.!] 
log£c[6c,  1] 

led,  1] 
5e[6rf,l] 
log  Be  [ftrf,  1] 

[em,  1] 

ft[ftm,l] 
log  ^c  [6m,  1] 

[«.  1] 
£[*,!] 
log^e[ft.,l] 

[rf^l] 
^[6rf,l] 
Iog^[W,l] 

[rfm,l] 
^rf[6m,l] 
log£d[6/»,l] 

[A.  1] 
-%[*..!] 
log  Dd  [6S,1] 

1    [mm,l] 
A.[Am,l] 
logf?m[im,  1] 

[m.,  1] 
A*  [6*.  1] 
logJBn[6,,l] 

log?* 
JogO. 

[ce,2] 
log  [re,  2] 

[erf,  2] 
log  [erf,  2] 

[em,  2] 
log  [em,  2] 

[«.  2] 
log  [«,  2] 

[rfrf,2] 
(?„[«*,  2] 
log  Cd  [erf,  2] 

[rfm,  2] 
Q[em,2] 
logCrf[«n,2] 

[A,  2] 
Ct  [e.,  2] 
log  Q  [e«,  2] 

[mm,  2] 
^[em.2] 
log<7M[cm,2] 

[«-kf] 

CU[«.«] 
log^^.2] 

log/?-    = 

log  -  *« 

[rfrf,3] 
log[rfrf,3] 

[rfm,  3] 
log  [rfm,  3] 

[A,  3] 
log  [A,  3] 

[mm,  3] 
£>„  [rfm,  3] 
'og  Dm  [rfm,  3] 

[m.,3] 

/)m  [rf^  3] 

Jog/?m[rf*,3] 

M  = 

[mm,  4] 

[m,,  4] 
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The  logarithms  of  the  quantities  in  the  first  row  of  each 
compartment  are  also  written  in,  and  from  these  by  proper 
subtractions  are  obtained  the  logarithms  in  the  margin,  where 

[«*]          A  [«c]  A  [am] 

[5aJ'  [oF'  m  ~ 

* 

• 


(L) 


Now  in  each  compartment  adding  the  logarithms  at  the 
margin  to  each  of  the  logarithms  in  the  first  row  of  that  com- 
partment we  obtain  the  corresponding  logarithms  written  in 
the  other  rows.  The  numbers  represented  by  these  logarithms 
are  next  written  above  them,  and  if  each  of  these  quantities 
is  then  subtracted  from  the  one  above  it  the  result  will  be 
the  corresponding  quantity  in  the  compartment  below.  Some 
of  the  squares  in  each  compartment  are  left  vacant  as  the 
quantities  belonging  to  them  have  already  appeared  above. 

111.  Application  of  Checks.  Also,  by  (93)  and  (94),  in 
the  first  compartment  the  quantities  in  the  first  lines  of  the 
last  column  should  be  equal  to  the  sura  of  all  the  quantities 
in  the  first  lines  of  the  corresponding  rows  plus  the  quantities 
similarly  situated  above  the  first  terms  of  the  rows.  Similar 
checks  will  apply  in  each  compartment;  for  if  from  the 
second  of  equations  (94)  we  subtract  the  product  of  the  first 
equation  multiplied  by  Ab  we  have 

[M,  1]  +  [Ac,  1]  +  .  .  .  [ftro,  1]    =    [fc,  1]     (97) 

In  the  same  manner  we  may  show  that  a  corresponding 
check  holds  throughout,  so  that  finally  we  shall  have 

[mm,  4]    —    [w«,  4]  (98) 
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The  last  compartment  of  the  table  is  added  to  give  this 
final  check  and  the  value  of  [uu]  in  accordance  witli  equa- 
tion (96). 

If  the  multiplications  and  divisions  are  simple  or  if  a  table 
of  squares  or  products  or  a  computing  machine  is  used  the 
logarithms  will  of  course  be  omitted  from  the  scheme  of 
solution. 

112.  Example.  In  order  to  illustrate  the  systematic  for- 
mation of  the  coefficients  that  appear  in  the  Normal  Equations 
as  well  as  the  solution  of  the  latter  by  the  above  method  we 
will  take  the 


OBSERVATION    EQUATIONS 


-  2   — 


=  0.1 
=  0.6 
=  0.1 


=  0.8 
=  —0.1 
=  0.4 


First  form  a  table  containing  the  coefficients  in  these  equa- 
tions and  also  the  sums  s.  As  a  first  check  the  sum  of  the 
quantities  in  column  s  should  be  equal  to  the  sum  of  all  the 
quantities  in  all  the  other  columns. 


I.  COEFFICIENTS  IN  THE  OBSERVATION  EQUATIONS. 


No. 

a 

b 

c 

d 

m 

s 

1 

-  1 

1 

I 

1 

-  .1 

1.9 

2 

1 

1 

—  1 

-  1 

-  .6 

—    .6 

3 

1 

2 

-  1 

1 

—  .1 

-  1.1 

4 

1 

—  1 

0 

-  2 

-  .3 

-  12.:; 

5 

0 

1 

I 

1 

.1 

1.1 

6 

1 

0 

-  1 

0 

-  .4 

.4 

Sum. 

3 

0 

-  3 

0 

-  1.4 

-  1.4 
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From  these  we  now  compute  the  coefficients  in  the  Normal 
Equations  (B),  and  also  the  necessary  quantities  for  the 
check  equations  (93)  and  (94). 


II.  COEFFICIENTS  IN  THE  NORMAL  EQUATIONS. 


No. 

aa 

ab 

ac 

ad 

am 

as 

bb 

be 

bd 

bm 

1 

1 

-1 

-1 

-1 

.1 

-1.9 

1 

1 

1 

-.1 

2 

1 

1 

-  1 

-  1 

-.6 

-    .6 

1 

-1 

-1 

-.6 

3 

1 

-2 

-1 

1 

-.1 

-1.1 

4 

2 

-2 

.2 

4 

1 

-1 

0 

-2 

-.3 

-2.3 

1 

0 

2 

.3 

5 

0 

0 

0 

0 

.0 

.0 

1 

-1 

1 

.1 

6 

1 

0 

-1 

0 

-.4 

-   .4 

0 

0 

0 

.0 

Sum 

5 

-3 

-4 

-3 

-1.3 

-6.3 

8 

1 

1 

-.1 

No. 

bs 

cc 

cd 

cm 

cs 

dd 

dm 

ds 

mm 

ms 

1 

1.9 

1 

1 

-.1 

1.9 

1 

-.1 

1.9 

.01 

-.19 

2 

-    .6 

1 

1 

.6 

.6 

1 

.6 

.6 

.36 

.36 

3 

2.2 

1 

-1 

.1 

1.1 

1 

-.1 

-1.1 

.01 

.11 

4 

2.3 

0 

0 

.0 

.0 

4 

.6 

4.6 

.09 

.69 

5 

1.1 

1 

-1 

-.1 

-1.1 

1 

.1 

1.1 

.01 

.11 

6 

.0 

1 

0 

.4 

.4 

0 

.0 

.0 

.16 

.16 

Sum 

6.9 

5 

0 

.9 

2.9 

8 

1.1 

7.1 

.64 

1.24 

If  the  coefficients  in  the  Observation  Equations  are  large, 
so  that  it  becomes  convenient  to  use  logarithms,  other  tables 
corresponding  to  I  and  II  would  be  formed  to  contain  these 
logarithms. 
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Substituting  these  quantities  now  in  the  general  tabular 
scheme  of  paragraph  110  we  have  the  results  on  (he  following 
page.  The  work  of  the  first  compartment  is  performed 
without  the  use  of  logarithms,  as  the  numbers  are  simple. 
The  quantities  in  the  last  column  should  all  be  zero  accord- 
ing to  the  check  equations,  what  small  differences  there  are 
I  Ting  due  to  the  rejection  of  figures  beyond  the  third  place 
in  the  decimals.  The  decimal  points  in  logarithms  to  which 
correspond  negative  numbers  have  been  replaced  by  the 
letter  n. 

The  demonstrations  that  have  been  made  now  enable  us  to 
see  at  once  from  an  inspection  of  the  results  in  this  table 
that 

2     =    0.238 


[iw]    =      .007 

Therefore  substituting  in  equations  (73)  and  (74)  we 
have 

fi^    =    .047  rZi    =    .032 

If  the  two  values  of  \_vv~]  obtained  in  the  solution  had 
differed  at  all  we  should  have  taken  the  mean  of  the  two. 

113.  If  the  Elimination  Equations  (95)  are  divided  by 
[««],  [bb,  1],  [cc,  2],  \_dd,  3],  respectively,  they  become 


Am    =  0 

#s   +   J£C23   +   ltd**   +   Em    =  0 

Cm     =  0 

Dw  =  0 


And  the  solution  for  the  unknowns  can  be  effected  most 
conveniently  by  arranging  the  computations  in  the  manner 
illustrated  on  page  109. 
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SCHEME  A.     SOLUTION  OF  THE  NORMAL  EQUATIONS. 


a 
b 

b 
C 
c 

d 

U 
1/1 

-m 
b 

C 

?c 
d 

*d 

m 

?m 
c 

d 

V 

m 

'm 

d 
m 

>m 

m 

a 

b              c              d            m            s          S 

5 

-3 

-4 

-3 

-1.3 

-6.3 

0 

-.6 
-.8 
-.6 
-.26 

8 
1.8 

1 
24 

1 
1.8 

-    .1 

.78 

6.9 

3.78 

0 

5 
3.2 

0 
2.4 

.9 
1.04 

2.9 
5.04 

0 

8 
1.8 

1.1 

.78 

7.1 

3.78 

0 
T 

.64 
.34 

1.24 
1.64 

9n3537 
9n1107 
9n1521 

6.2 
0.7924 

-1.4 
On1461 

-    .8 
9*9031 

-    .88 
9,,9445 

3.12 
0.4924 

0 

1.8 
.316 
9.4998 

-2.4 
.181 
9.2568 

-    .14 
.199 
9.2982 

-2.14 
-    .705 
9n8479 

6.2 
.103 
9.0138 

.32 
.114 
9.0552 

3.32 
-    .403 
9n6049 

.30 
.125 
9.0966 

-    .40 
-    .443 
9n6463 

On2403 
9W3587 

1.484 
0.1715 

-2.581 
On4118 

-    .339 
9n5302 

-1.435 
On1569 

1 
T 

T 

6.097 
4.489 
0.6521 

.206 
.589 
9.7705 

3.373 
2.496 
0.3972 

.175 
.077 

8.8889 

.043 
.328 
9.5156 

9W8769 

=  Iog-z4 

1.608 
0.2063 

-    .383 
9n5832 

1.227 

0.0888 

2 

.098 
.091 
8.9601 

-    .285 
-    .292 
9n4657 

0 

0.238 

=  «4            [w]  =  .007 

.007 

.007 

0 
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SCHEME  B. —  SOLUTION  OF  THE  ELIMINATION  EQUATIONS. 


-°m 

-cm 

:J: 

-  Am 

—   Ae   Z3 

24 

Iogz4 

Iogz3 

Iogz2 

«1 

log  Cd 

log  Ac 

log^ 

log  Ab  z2 

log  7?c  z3 
log  Ac  z8 

log  Cdz4 
logJ?rfz4 

114.    Filling  out  this  table  for  the  example  just  solved,  we 
have 


.238 

.228 

.142 

.260 

.414 

.031 

.143 

.145 

.514 

.191 

.238 

.642 

.318 

1.108 

9.3769 

9.8075 

9.5024 

On2403 

9n3537 

9n7782 

9,,1107 

9n9031 

9n7782 

9n2806 

9n7106 

9n6172 

8'\  Q7A 

n4o7o 
97J1551 
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115.  The  Weights  of  the  Unknowns.  In  order  to  deter- 
mine the  precision  measures  of  zt,  22,  and  23,  it  would 
next  be  necessary  to  compute  the  weights  of  the  latter  quan- 
tities. The  demonstration  of  the  processes  by  which  these 
weights  may  be  found  will  not  be  taken  up  here,  as  the  best 
method  to  adopt  varies  a  good  deal  with  the  character  of  the 
example,  but  a  statement  of  the  results  in  the  general  form 
of  solution  will  be  given. 

By  treating  the  Elimination  Equations  in  a  way  similar  to 
that  used  in  deriving  equation  (E)  of  paragraph  76  from 
equations  (B)  of  the  same  paragraph,  we  may  show  that 


Zl   +   Am   -f    JBmai   +    Cma2   +   J)mas     =     0 

2,  +  Hm       +  Cm  ft  +  7>m&  =    0  (100) 

z.  +   Cm       +  Dm  73    =    0 

where  the   a's,    /J's,    y's,    are  determined  from  the  equations 

Ad  +  J3da,  +  Cda,  +  a,  =   0       J?d  +   Cd  &  +  &   =   0 
-4,+JSe0l  +  a,  =   0  J?c+&  =   0  (101) 

Ab  +   a,  =   0  Cd   +  73  =   0 

Then  by  an  application  of  the  principles  of  Rule  1,  para- 
graph (77),  it  may  be  shown  that 


i 

X 
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These  equations  can  of  course  be  extended  to  cover  any 
number  of  unknown  quantities,  and  tabular  schemes  for  the 
computations  of  the  a's,  /?'s,  y's,  .  .  .  and  of  the  weights 
Pzj  PZJ  PZJ  ...  can  readily  be  arranged. 

For  a  general  demonstration  of  these  results,  and  also  for 
a  discussion  of  special  methods  of  solution,  consult  — 

Johnson,  "The  Theory  of  Errors  and   Method  of  Least  Squares," 

chap.  ix. 

Wright,  "  Treatise  on  the  Adjustment  of  Observations,"  chap.  iv. 
Chauvenet,  "  Spherical  and  Practical  Astronomy,"  pp.  530-649. 


THE    METHOD   OF   CORRELATIVES. 

116.  The  method  of  adjusting  "Conditioned  Observa- 
tions "  explained  in  paragraph  33  is  perfectly  general,  but 
where  there  are  many  conditions  to  be  satisfied  the  solution 
is  apt  to  be  very  laborious.  For  the  case  that  occurs  most  fre- 
quently in  practice,  in  which  the  observations  are  direct  and 
equal  in  number  to  the  number  of  unknown  quantities,  the 
process  of  solution  devised  by  Gauss  and  called  the  "  Method 
of  Correlatives  "  is  the  most  convenient.  This  method  is 
derived  as  follows  :  — 

Let  q  observations,  JtfJ,  M^  .  .  .  Mq,  of  the  respective 
weights  PU  PK  .  .  .  pq,  be  made  directly  upon  the  values 
of  q  unknown  quantities,  and  let  the  most  probable  values 
of  the  unknowns  be 


Where  v1?  v2,  .  .  .  v?,  are  the  most  probable  corrections  to 
apply  to  the  observed  values  as  well  as  in  this  case  the 
residuals  of  the  observations. 

If  the  n'  condition  equations  are  not  linear  they  may  be 
reduced  to  that  form  by  the  method  of  paragraph  44,  so  that 
we  may  assume  for  our 
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CONDITION    EQUATIONS 


.  .  .  aqvq  +  ml    =    Q 
.  .  .  bqvq  -f  m2    =    0 
.........  (A) 

ll»l   +    l*V*    +    •  •  •    IqVq    +    mn>    =      0 

In  which  the  quantities  m^  7W2,  .  .  .  mn',  would  all  be  zero 
if  the  observations  were  exact.  It  is  to  be  observed  that  the 
coefficients  a,  #,  .  .  .  I  are  not  arranged  in  the  same  order  in 
these  equations  as  they  are  in  the  observation  equations  of 
paragraph  107. 

The  values  of  w1?  v2,  ...  vq,  must  be  determined  so  as  to 
satisfy  the  above  equations  and  also  by  the  principle  of  Least 
Squares,  so  as  to  make 


2*  +  •  •  •  Pqvg*    =    a  minimum. 
Corresponding  with  a  minimum  value  of  this  last  we  have 

PlVl<t»l    +  PtVidVt   +    .  .  .  PqVqdvq     =      0  (B) 

for  all  possible  simultaneous  values  of   dv^  dvz,  .  .  .  dvq  ;  that 
is,  for  all  values  which  satisfy  the  equations, 

aldvl  -\-  a2dvz  +  .  .  .  aqdvq    =    0 
bidvi  -\-  t>zduz  -f-  .  .  .  bqdvq    =    0 

.........  (C) 

lidVi  -[-   I*dv2  -h   .  .  .    lqdvq    =    0 

obtained  by  differentiating  equations  (A). 

Therefore,  di-ii..tiiiLr  tla-  lirst  member  of  (B)  by  R  and  the 
first  members  of  (C)  by  A\,  xS'2,  .  .  .  AVn.,  it  will  be  IUT- 
CHsary  that 

=    0  (D) 


\\hcre    X*j,  /*,,   .  .  .   /•„.,    an-  uiidi'tiTiniiifd  c.,,  -ilici«-nts. 
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This  last  equation  will  be  satisfied  if  the  coefficient  of  each 
differential  in  it  is  made  equal  to  zero,  that  is,  if 


=    ^a2  -f-  &2£2  -f-  ...  kn>l2 

(103) 

=    k^aq  +  kzbq  -f  ...  Jcn,lq 


All  that  remains  therefore  is  to  find  values  of  vlt  V2,  ...  vq 
and  &!,  &2,  ...  kn,>  which  will  satisfy  simultaneously  equa- 
tions (A)  and  (103),  and  that  this  may  be  done  is  easily 
seen  from  the  fact  that  we  have  the  same  number  of  equations 
as  unknowns. 

Substituting  the  values  of  v^  vz,  ...  vq  from  equations 
(103)  inequations  (A)  we  have  the  following  : 


(id  —  (ib  _  dl 

i  S  —  +  K*  S  p -  +  .  .  .  kn  S  —  +  ml  =  0 


_  ab  ^bb  _  bl 

,  S  p    +  A-2  S  ~  +  . .  .  Av  S  —  +  m2  =  0 


(104) 


S  —   +  .  .  .  /?„'  S  — 


The  solution  of  these  equations  gives  at  once  the  values  of 
£M  /»•«,  ...  Arn',  which  are  called  the  "Correlatives"  of  the 
Condition  Equations.  The  values  of  v^  v^  .  .  .  vq  are  then 
found  by  substituting  the  values  of  the  tfs  in  equation  (103). 

117.  As  equations  (104)  are  of  the  same  general  form  as 
a  set  of  Normal  Equations,  Gauss's  Method  of  Substitution 
can  be  advantageously  employed  in  the  solution. 
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When  there  is  but  a  single  equation  of  condition  the  second 
members  of  equations  (103J  reduce  to  their  first  terms,  and 
equations  (104)  reduce  to  the  single  equation 

CLd 

k^  —  -f  ml    =    0  (105) 

and  the  values  of   vl9  «2,  ...  vq   in  (103)  become 


m" 


a, 

V*-  mi 
aa 

p 

It  is  from  these  results  that  the  rules  in  paragraph  35   are 
derived. 

118.    Example.      Suppose  we  have  given  the  observations 

Ml  =  2.02,  weights 

Mt  =  4.13,  "       2 

Jf,  r=  2.52,  «       5  (a) 

M.  =  2.67,  »       7 

Jf6  =  2.84,  «       4 

and  let  the  most  probable  values  of  the  unknowns  be  repre- 
sented by 

Z,    =    J/i  +  Wlf       2,    =    MI  +  V2,    .   .   .    26    =    MS  +  V8      (b) 

Also  suppose  that  the  unknowns  are  subject  to  the  conditions 

•i  +  »•  +  «•  +  z*  +  2.    =    14.0 

=      1.5 
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Then  expressing  these  conditions  in  terras  of  the  corrections 
by  means  of    (a)    and    (b),    we  have  the 


CONDITION    EQUATIONS 

Vi  +  v<t  +  va  +  v4  +  v.  +  .18    =    0 
v  -   v  -   .04    =    0 


Referring  to  paragraph   116,   we  see  that  in  this  example 
n'    =    2,         ml    =    0.18,         m2    =     —  0.04 

For  the  purpose  of  computing  the  coefficients  in  equations 
(104)    we  next  arrange  the  following  table. 


p 

a 

b 

aa 
P 

ab 
P 

Ib 
P 

1 

3 

1 

0 

•  — 

0 

0 

3 

1 

1 

1 

2 

1 

1 

2 

2 

2 

1 

5 

1 

0 

A 

0 

0 

5 

1 

1 

1 

7 

1 

\ 

—  — 

7 

7 

7 

i 

4 

1 

0 

X 

0 

0 

4 

599 

5 

9 

420 

14 

14 
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Substituting  these  results  in  equations   (104),  we  have 

Cl  +    14  **  +  '^ 

14  *l  "      14  *'  " 

Solving,  &j    =   —  .1647  (,. 

h    =        .1537 

Then  from  equations   (103)   we  get  at  once 

Vl    =    -  .0549,         04    ss    —  .0455, 

va    =    -  .0055,         v&    =    -  .0412.  (g) 

vs    =    -.0329, 

And  by  substituting  these  results  in    (b)   we  can  obtain  the 
values  of   zx,  z2,  z3,  z4,  z&. 

The  above  is  the  solution  of  Example  70,  page  127. 


EXAMPLES. 

1.  An  urn  contains  five  black  balls,  three  red  balls  and 
two  white  balls.     1 1'  three  balls  are  drawn  from  the  urn  what 
different  combinations  may  result,  and  what  is  the  probability 
of  each  ? 

2.  In  a   single  throw  with   a  pair  of  dice  what  is   the 

probability  that  neither  ace  nor  doublets  will  appear  ? 

9 

3.  Four   cards  are   drawn    from  a  pack.     What   is   the 
probability  of  getting  four  aces?     Of  getting  one    of   each 
suit? 

4.  From  a  lottery  of  thirty  tickets,  marked  1,  2,  ...  30, 

four  tickets  are  drawn.     What  is  the  probability  that  num- 

o 
bers  1  and  15  will  be  among  them? 

145 

5.  Find  the  odds  against  the  appearance  of  7  or  11  in  a 
single  throw  with  a  pair  of  dice.  7  :  2 

6.  I  toss  up  n  coins.     What  is  my  chance  of  getting  just 
one  head  ? 

7.  In   a  single  throw  what  are  the  relative   chances  of 
throwing  9  with  two  dice  and  with  three  dice  ?  24  :  25 

8.  From  2  n  counters  marked  with  consecutive  numbers 
two  are  drawn.     What  are  the  odds  against  having  an  even 
sum  ?  n  :  n  —  1 

9.  In  two  trials  with  a  single  die  what  is  the  probability 
of  throwing  (a)  an  ace  the  first  time  only?    (b)  at  least  one 
ace? 

10.  Find  the  probability  of   throwing   doublets    one   or 

91 

more  times  in  three  trials  with  a  pair  of  dice. 

216 
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11.  Find  the  probability  of  throwing  exactly  three  aces 

125 
in  five  trials  with  a  single  die.  

3888 

12.  A  certain  stake  is  to  be  won  by  the  first  person  who 
throws  5  with  a  die  of  twelve  faces.     What  is  the  chance  of 
the  sixth  person  ? 

13.  A  and  B  play  chess.     A  wins  on   the  average  two 
games  out  of  three.     What  is  A's  chance  of  winning  just 

80 
four  games  out  of  the  first  six  ? 

243 

14.  A  and  B  shoot  alternately  at  a  mark.     A  hits  once  in 
n  times  and  B  once  in  n  —  1  times.     Find  their  chances  of 
first  hit,  and  the  odds  in  favor  of  B  if  A  misses  on  his  first 
shot.  Even,     n  :  n  —  2 

15.  In  how  many  trials  will  it  be  a  wager  of  4  to  3  that 
double  five  will  be  thrown  with  a  pair  of  dice  ?  30 

16.  Find  the  probability  of  throwing  one  and  only  one 

ace  in  two  trials  with  a  single  die. 

17.  If  I  have  three  tickets  in  a  lottery  of  four  prizes  and 

41 

eight  blanks,  what  is  my  chance  of  drawing  a  prize  ? 

55 

18.  Find  the  probability  of  throwing  at  least  four  aces  in 

six  trials  with  a  single  die.  

23328 

19.  On  an  average  seven  ships  out  of  eight  return  to  port. 
Find  the  chance  that  out  of  five  ships  expected  at  least  three 

16121 

will  return. 

16384 

20.  In  a  lottery  containing  a  large  number  of  tickets, 
where  the  prizes  are  to  the  blanks  as  1  :  9,  find  the  chance 

81  J.R 

of  drawing  at  least  two  prizes  in  five  trials.  

100000 
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21.  In  a  purse  are  ten  coins,  all  nickels  except  one  which 
is  a  five-dollar  gold  piece ;    in    another   are    ten   coins,  all 
nickels.     Nine    coins    are    taken    from    the    first   purse    and 
placed  in  the  second,  and  then  nine  coins  are  taken  from  the 
latter  and  placed  in  the  former.     If  you  now  had  your  choice 
which  purse  would  you  take  ? 

22.  A  and  B  engage  in  a  game  in  which  A's  skill  is  to 
B's  as   2:3.     What  is  A's   chance  of  winning  at  least  two 
games  out  of  five  ? 

23.  If   A's   skill  at  a  certain  game  is  double  that  of   B, 
what  are  the  odds  against  A's  winning  four  games  before  B 
wins  two?  131  :  112 

24.  A  party  of  twenty-five  take  seats  at  a  round  table. 
What  are  the  odds  against  any  two  specified  persons  sitting 
next  to  each  other  ? 

25.  A  has  three  shares  in  a  lottery  in  which  there  are  three 
prizes  and  six  blanks.     B  has  one  share  in  another  where 
there   is   but   one   prize  and  two  blanks.     What  are   their 
relative  chances  of  getting  a  prize  ?  A  :  B  =  16  :  7 

26.  Expand  through  the  terms  involving  A8  and   &8,  the 
expression 

-  +  (y  +  *)• 

X  -\-  h 

When    a;    is   1    and    y    is    -,   does   -  +  ^*   increase   or 

diminish  when  x  and  y  begin  to  increase  at  the  same  rate  ? 

27.  Given  f  (x,  y)    =   x2(a  -f-  y)8,  expand  the  expres- 
sion   (x  +  hY(a  4-  y  +  &)8. 

28.  Find  the  value  of   Iog10  a  -\-  cos  #,  when    a  =  1001 
and    b  =  0.1°.     Give  the  result  first  to  five  places  and  then 
to  seven  places  of  significant  figures,  in  each  case  without  the 
aid  of  tables.  4.000433 

29.  Transform  to  the  new  origin,    (2,  —  3,  1),    the  equa- 
tion,  x2  +  if  +  z*  _  4  x  +  6  y  -  2  z  -  11  =  0,    the 
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axes  remaining  parallel  to  the  original  ones.  The  equations 
of  transformation  are  x  =  2  +  a;' ,  y  =  —  3  -{-  y' , 

2=1  +  2'.  iC2  +  y2  +  22   =    25 

30.  Find  the  minimum  value  of  x2  -|-  xy  -\-  y*  —  ax  —  by. 

(ab  —  a*  —  tf) 
3 

31.  Find  the  values  of  a:,  y  and  z   that  render  a  maximum 
or  a  minimum  the  function    ce2  -|-  y2  -|-  z2  -f-  a;  —  22  —  av/. 

32.  Find  the  values  of   a?   and   y   that  render  a  maximum 
or  a  minimum  the  expression   sin  x  -\-  sin  y  -f-  cos  (x  -\-  y). 

33.  Find  the  co-ordinates  of  a  point  the  sum  of  the  squares 
of  whose  distances  from  three  given  points,    (a^,  y,),  (a*2,  y2), 

( a?i  -f-  a?2  -I-  a;8 ) 
(xs,  y3),    is  a  minimum.  a;   =   - 

3 

34.  Given  the  volume,  a8,  of  a  rectangular  parallelepiped, 
find  its  shape  when  its  surface  is  a  minimum. 

35.  Find  the  volume  of  the  greatest  rectangular  parallelo- 
piped  that  can  be  inscribed  in  the  ellipsoid 

av*          y-          2* 
a2      "   fta      "   c2 

36.  In  the  Physical  Laboratory  apparatus  for  illustrating 
the  Estimation  of  Tenths,  the  reading  of  the  micrometer  head 
for  a  certain  setting  is  computed  to  be  2.3038,  which  may  be 
taken  as  exact.     Setting  the  apparatus  by  the  eye  the  follow- 
ing readings  are  obtained  :  — 

2.314         2.324         2.310         2.519         2.326 
2.320         2.302         2.313         ±305 

What  are  the  accidental  and  real  errors  and  the  residuals? 
Are  there  any  constant  errors  or  mistakes?  If  then-  arc 
< •.. n^taiit  errors  are  they  of  the  first,  second,  or  third  class? 
I!  HT  do  you  tell? 
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37.  Are  the  following  observations  such  as  to  call  for  an 
application  of  the  Method  of  Least  Squares  in   their  :i<lju>t- 
mriit? 

x  -\-    y  —     z-j-w—         5  3ic  —  y  —     z  —    u   =    I 

x  —  2y  -\-  2z  -f-  u   =    —  1  x  —  y  -\-  4z  +  6u    =    9 

38.  In  the  case  of  direct  observations,  what  other  quanti- 
ties  besides    the    arithmetical    mean    might    reasonably   be 
assumed  to  give  plausible  values  of  the  unknowns  ?     Why  is 
the  arithmetical  mean  preferred  to  these? 

39.  Find  the  most  probable  value  of  a  quantity   M  from 
the  observations 

216.27         216.16         216.04         216.19         216.44         216.58 
.29  .43         215.99  .39  .51 

.33  .09         216.23  .14         215.94 

Also  test  the  result  by  finding  the  sum  of  the  residuals. 

M  =   216.251 

40.  In  the  determination  of  a  certain  wave  length,  Row- 
land made  the  following  observations.     Find  the  most  prob- 
able wave  length. 

4.524  4.515  4.513  4.507  4.501  4.485  4.517  4.493  4.505 
.500  .508  .511  .497  .502  .519  .504  .492 

4.5055  ±  0.0017 

41.  Ten  measurements    of  the  density  of   a   body  gave 
results  as  follows.     Find  the  most  probable  density. 

9.662    9.664    9.677    9.663    9.645 
.673     .659     .662     .680     .654 

9.6639  ±  0.0022 

42.  In  a  triangulation  of  the  U.  S.  Coast  Survey  an  angle 
was  measured  twenty-four  times  with  results 
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116°  43'  44".45 
50  .55 
50  .95 


51 
51 

51 

48". 
49  . 

48  . 

".30 
.05 
.70 

90 
20 

85 

49 
50 
49 

47" 
47 
51 

".05 
.55 
.25 

.40 
.75 
.05 

46 
49 
53 

47".85 
50  .60 
48  .45 

".75 
.25 
.40 

51 
49 
52 

".75 
.00 
.35 

Find  the  best  value  for  the  angle.         116°  43'  49".64  ±  0".28 

43.  The  following  observations  were  made  with  a  sextant 
in  order  to  determine  the  latitude  of  a  place  :  — 

43°  4'  46"         7"       59"       52"       47"       36" 
24         28         39         52         15         40 

What  is  the  most  probable  latitude?  43°  4'  27"  ±  5" 

44.  Determine  the  quantity    M   from  the  observations 

M  =   81.55    .41    .68    .25    .27    .77    .13    .86    .10    .03 
Wts.  25     25     16       9       9       9       4       1       1       1 

45.  An  angle    M  is  measured  with  the  following  results  : 

M  p          M  p         M  p         M  p         M  p 

45".00  5  42".50  5  27".50  3  36".25  2  45".00  2 

31  .25  4  37  .50  3  43  .33  3  42  .50  3  40  .83  3 

45  .00  3  38  .33  3  40  .63  4  39  .17  3 

Find  the  most  probable  value  of  the  angle  and  test  the 
result  by  computing    ^po.  39".78  ±  0".94 

46.  In  one  hundred  measurements  of  angles  made  in  the 
primary  triangulation  in  Massachusetts  by  the  Coast  Survey 
there  were  found  between 

6".0   and   5".0       1  error.  1".0    and      0".0     26  errors. 


5  .0 

M 

4 

.0 

2  errors. 

0 

.0 

u 

—  1 

.0 

26 

u 

4  .0 

M 

3 

.0 

2      " 

-  1 

.0 

u 

-2 

.0 

17 

u 

3  .0 

u 

o 

.0 

3      " 

-2 

.0 

u 

-3 

.0 

8 

u 

2  .0 

u 

1 

.0 

13      « 

-3 

.0 

u 

—  4 

.0 

2 

u 

Plot  these  results  as  in  Figure  1,  page  8. 
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47.  In  sixty-six  determinations  of  the  velocity  of  light 
made  at  Washington  the  percentage  of  errors  of  different 
magnitudes  was  found  to  be  as  follows:  — 
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to  -  7.0 

.8 
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-5.0 

2.0 

5.0 

4.2 

-4.0 
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4.0 

5.8 
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7.8 

3.0 

8.6 

-2.0 

12.8 

2.0 

12.4 

—  1.0 

17.1 

1.0 

17.1 

Plot  the  results  and  draw  a  smooth  curve. 

48.  Given  the  observations 

x  -       y  4-  22    =    3  4<e  -f     y  -f  4z    =    21 

3X  _|_  2y  --   52    =    5         -  x  +  3y  -f  $z    =    14 

find  the  most  probable  values  of   a;,   y,    and    z. 

x    =    2.470  ±  .038 

49.  O  being  the  level  of  the  sea  and   P15  P2,  jP3,   three 
points  whose  altitudes  are  to  be  determined,  the  following 
observations  are  made  :  — 

Pl  above  O    =    10  ft.         P2  above  Ps    ==    9  ft. 
P2       "      PI  =      1  PI      "       P8    =    2 

P2       «      0    =    18 

Find  the  most  probable  altitudes.  Ps    =    8.50  ±  .29 

50.  The  altitudes  of   A   above    O,   B   above   A,   and  B 
above    O  are   found   by    measurements   to   be   respectively 
12.3,  14.1,  and  27.0  feet.     What  is  the  most  probable  value 
of  each  of  these  differences  in  level?         A   =.    12.50  ±  .17 

51.  Measurement  of  the  ordinates  of  points  on  a  straight 
line  corresponding  to   abscissas  4,  6,  8,  9,  are   made    with 
results  5,  8,  10,  12.     What  is  the  most  probable  equation  of 
the  line  in  the  form   y  =   mx  -f-  b?  b  =    —  0.29 
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52.  Find  the  altitudes  in  Example  49  if  the  observations 
have  weights  5,  3,  6,  2,  4,  respectively. 

53.  Solve  the  example  in  paragraph  31,  giving  the  obser- 
vations the  weights  25,  25,  4,  4,  4,  4,  4,  4,  1,  respectively. 

Elevation  of   PB  =   320.25 

54.  Find  the  most  probable  values  of   2lt  22,  and  z8  from 
the  observations 

zl  =  552.10  wt.  16  21  —  22   =         .75  wt.  1 

-  z2  4-  28  =         .15    "      9     2X  +  22  -  28   =   552.05    "    1 

28  =  551.23    "4  2l  -  28   —         .70     "    1 

22  =  551.30    «      4 

22   =   551.2345 

55.  In  the   triangulation   of  Lake  Superior   there  were 
measured  at  station  0  the  angles 
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Required  the  adjusted  values  of  the  angles. 

F  O  P  =   40".28  ±  0".34 

56.     In  the  U.  S.  Lake  Survey  the  following  angles  were 
measured  at  station  North  Base  :  — 

(1)  Crebassa  — Middle                  55°  57'  58".68  wt.    3 

(2)  Middle  — Quaquaming           48    49    13  .64  "    19 

(3)  Crebassa  — Quaquaming      104    47    12  .66  "    17 

(4)  Quaquaming  — South  Base    54     38    15.53  "    13 

(5)  Middle  — South  Base            103     27    28.99  "      6 

Find  the  adjusted  values  of  the  angles. 

(1)    =   58".965;    r  =   0".28 
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57.  Adjust  the  following   observations  of  differences  in 
level: — 

Altitude  of  A    401.3    wt.  16  C  above  B      72.5    wt.  9 

A  above  B    220.8      «    16  A     "       B    222.0      "    1 

^4      «       C    150.2      "      4        Altitude  of  B    180.7      "    1 

58.  In  "Conditioned  Observations"  can  the  number  of 
observations  required  be  less  than  the  number  of  unknown 
quantities?     Why  must  the  number  of   conditions  be  less 
than  the  number  of  unknowns? 

59.  From  the  following  measurements  of  the  angles  formed 
at  the  centre  of  a  disk  by  four  radial  lines,  find  the  most 
probable  values  of  the  angles. 

A    —    104°  25'  13"  C    =    86°  33'  20" 

B    =      98     13    47  D    =    70     48    23 

A    =    104°  25'  2".25 

Also  solve  giving  the  observations  the  weights  5,  2,  1,  4, 
respectively. 

60.  Four  observations  on  the  angle   A   of  a  triangle  gave 
a  mean  of   36°  25'  47",    two   observations   on    B    gave   a 
mean  of   90°  36'  28",    and  three  on    C  gave    52°  57'  57". 
Adjust  the  triangle.  A   =   36°  25'  44".2 ;    r  =   7".  7 

61.  Five  angles  at  a  station  are  measured,  and  also  their 
sum.     The  observed  sum  differs  from  the  sum  of   the   five 
observed  parts  by  the  amount   d.     What  are   the   adjusted 
values  of  the  angles  ? 

62.  The  three  angles  of  a  spherical  triangle  are  measured 
with  results 

A   =  46°  17'  38".32          B  =   73°  35'  16".15 
C  =   60°  7'  5".16. 

Adjust  the   triangle,  knowing  that  the   spherical    excess  is 
2".475.  A   =   39".3;    ^  =    1".6 
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63.     At  the  station  Pine  Mountain  the  following  angles 
were  observed  between  surrounding  stations  :  — 

Jocelyne  — Deepwater  65°  11'  52".500  wt.  3 

Deepwater  —  Deakyne  66  24    15  .553  "    3 

Deakyne  — Burden  87  2    24  .703  "    3 

Burden— Jocelyne  141  21    21.757  "    1 

Find  the  most  probable  values  of  the  angles. 

04.     Solve  Examples  55  and  56  by  the  method  of  "Con. 
ditioned  Observations." 

65.  .lisa  station  whose  altitude  is  known  to  be  5240.1 
feet,     JB  and  C  are  floats  on  a  lake,  and  D  is  a  signal  point. 
From  the  following  observations  determine  the  most  prob- 
able altitudes  of  JB,   C  and  D. 

C  below  A  720.1  wt.  3  B  below  A  719.7  wt.  3 
D  «  A  200.3  "5  JB  «  J)  520.9  "  2 
C  «  D  520.4  "  2 

66.  Given  the  following  observations,  subject  to  the  con- 
dition   zl  -\-  za    =    z8,    find   the   most  probable  values  of 
zl9  z,,  and  z,. 

22j  —     za  +     z,    =        3.0  2z2  —     z8    =    1.0 

2zx  -  3za    =  -  4.5  Zj  +  2z2    =    5.1 

Zl  +     z8    =        3.8 

67.  The  chemical  composition  of  a  specimen  was  found 
by  several  observers  to  be  as  follows :  — 

Pb  =  .52      Other  substances     =  .09      Au  and  Ag  =  .39 
Ay  =  .27      Pb  and  Ag  =  .78      Impurities    =  .10 

Au  =  .11      Pb  and  impurities  =  .62      Au  =  .12 

From  these  observations  find  the  most  probable  composition 
m  the  specimen. 
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68.  From  the  following  observations  what  are  the  best 
values  of  the  unknowns,  supposing  that   y   and   z   must  be 
equal? 

x  +  y    =    5.2     wt.  4         y  +  z    =    4.2     wt.  I 
x    =    3.0      "    9  2    =    2.0       "   4 

x  —  z    =    1.1      "1 

69.  In  determining  the  difference  in  longitude  between 
various  cities  the  results  obtained  were 

(1)  Cambridge  —  Washington  23™  4K041  wt.  30 

(2)  Cambridge  —  Cleveland  42     14.875  «      7 

(3)  Cambridge  — Columbus  47     27.713  "      8 

(4)  Washington  —  Columbus  23     46.816  "      7 

(5)  Cleveland  —  Columbus  5     12.929  "      5 

Adjust  these  observations. 

70.  The  capacity  of  a  condenser  is  known  to  be  14.0  m.  f. 
It  is  divided  into  five  sections,  a,  5,  c,  </,  e,  and  it  is  known 
that  the  difference  between    b   and   d  is    1.5  m.  f.     Find  the 
most  probable  capacities  of  the  sections  from  the  observa- 
tions 

a    =    2.02       wt.  3  d   =    2.67       wt.  7 

b    =    4.13        "2  e    =    2.84        "    4 

c    =    2-52        "    5  a    =    1.9651 

71.  If  the  unknowns  in    the  following  observations  are 
subject  to  the  condition   x  -f-  2y  -|-  3z    =:    36,   what  are 
their  adjusted  values? 

x  =  4.3     wt.  1,        y  =  5.7     wt.  4,       z  =  7.3     wt.  9 

x    =    3.77 

72.  A  cannon  is  discharged  horizontally  from  the  top  of 
a  bluff.     Observations  on  the  time,  and  distance  of  fall  of  the 
ball  gave  the  results 

t    =    0.5       1.0       1.5         2.0    seconds 
s    =    1.2       4.0       9.1       15.0    metres 
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What  curve,  passing  through  the  point  of  departure  of  the 
ball,  will  represent  the  above  observations  ? 

73.  An  Argand  burner   shows  the  following  efficiencies 
with  varying  rates  of  gas  consumption  : 

g    —    2.0     2.3     2.8     3.3     4.0     4.5     5.0     feet 
E  =    2.1     2.4     2.5     3.0     3.2     3.8     4.1 

Find  the  equation  of  the  straight  line  which  best  rep- 
resents the  relation  between  g  and  E.  The  measurements 
on  g  are  without  appreciable  error. 

74.  Observations  are  made  upon  the  expansion  of  Amyl 
alcohol  with  change  in  temperature  as  follows :  — 

V  =    1.04       1.12       1.19       1.24       1.27     cu.  cm. 
t    =    13.9       43.0       67.8       89.0       99.2     C.  degrees 

If  V  =  1  -f-  B  t  -\-  C  t2  expresses  the  law  connect- 
ing the  volume  and  temperature,  find  the  most  probable 
values  of  B  and  C. 

75.  In  a  Hooke's  joint  where  the  angle  between  the  axes 
is  45°,    x  being  the  angular  rotation  of  the  driver,  and    y 
that  of  the  follower,  from  the  following  measurements  find 
the   equation   of   a   curve   that  will  represent   the   relation 
between   x  and   y  —  x. 
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y  —  x  =   -  -  0.85  —  9.82  sin  2aj  -f  0.92  cos  2a 

76.  A  series  of  observations  extending  over  a  period  of 
thirty  years  was  made  by  Quetelet  to  determine  the  daily 
variation  in  temperature  at  Brussels.  The  mean  results  of 
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the  measurements  are  given  below.     From  them  derive  an 
equation  to  express  the  temperature  at  any  time  of  tin-  year. 

Jan.  4°.66  May      9°.83  Sept.  8°.16 

1-Vb.  5°.42  June  10°.09  Oct.  6°.55 

Mar.  6°.77  July      9°.71  Nov.  5°.10 

Apr.  8°.59  Aug.     9°.14  Dec.  4°.41 

y    =    7.369  +  0.9854  sin  30*  -  2.7084  cosSOa; 

+  0.0100  sin  60s  —  0.1950  cosGOa; 

-   0.0133  sin  90x+ 0.1783  cos  90s; 

In  this  answer  the  values  of   x  begin  at  the  15th  of  Janu- 
ary, and  represent  the  time  in  months. 

77.  The  law  connecting  the  time  of  vibration  of  a  pendu- 
lum with  its  length  is  assumed  to  be  of  the  form,    T  •==.  m  Ln. 
From  the  following  observations  find  the  most  probable  values 
of   m   and   n. 

T   =      12.9       11.6       10.4       9.7       5.3       4.6 
L    =    164.4     132.9     107.6     93.5     28.4     20.6 

L   is  in  centimetres,    T  in  tenths  seconds.          n   =   0.5000 

m   =   1.0044 

78.  Determine  the  equation  of  a  curve  which  will  repre- 
sent the  following  observations : 


X 

0.0 

y 

0.00 

X 

1.5 

y 

1.09 

X 

3.0 

y 

8.65 

0.5 

0.04 

2.0 

2.56 

3.5 

13.72 

1.0 

0.31 

2.5 

4.99 

4.0 

20.47 

79.     Determine  the  equation  of  a  curve  which  will  repre- 
sent the  relation  between    x   and    y   in  the  observations, 
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80.  At  a  station  P  the  angles  between  a  straight  line 
passing  through  P  parallel  to  the  axis  of  X  and  the  direc- 
tions from  P  of  four  points  P^  P2,  P8,  P4,  are  measured. 
Having  given  the  coordinates,  (a,  £),  of  the  four  points, 
find  the  coordinates  of  P. 

Point.  Coordinates. 


a 
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PI 
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39°  18' 

PI 

1.21 
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147°  54' 
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If  the  coordinates  of  the  point    P    are    (a;,  y),    and  the 
angle  is  denoted  by   -4,    we  have 


x  —  a 


81.  If   asinbx  =   Jf/j   and  values  of   M  are   observed 
for  known  values  of   a   and   #,    determine  the  most  probable 
value  of   x.     If   x'   is  an  approximate  value  of   x   found  by 
trial,  and    m  =   a  sin  bx'  —  M,    we  shall  have 

2  a  b  m  cos  bx' 
2  («  b  cos  bx')* 

82.  If  in  one  series  of  observations  the  value  of    h   is 
twice  what  it  is  in  another,  what  is  the  relative  probability  of 
the  occurrence  of  an  error  of  given  magnitude   a  in  the  two 
series  ?     Show  what  the  curves  of  error  will  be  in  the  two 
cases.    What  error  has  the  same  probability  for  its  occurrence 
in  each  series?    What  is  the  relative  probability  of  the  occur- 
rence of  an  error  not  greater  than  a  in  the  first  case  and  not 
greater  than   2a   in  the  second  case  ? 

83.  From  64  observations  the  latitude  of  a  station  was 
found  to  be   49°  10'  9".110  ±  0".051.     What  was  the  prob- 
able error  of  a  single  observation  ?  0".41 
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84.  If  twenty  measurements  of  an  angle  give  a  result  with 
an    A.D.    of   0".38,    ami  it  is  required  to  timl  the  niiLrle  so 
that  the   A.D.   shall  be  only  0".25,  how  many  more  observa- 
tions must  be  made?  27 

85.  From  the  following  determinations  of  the  area  of  a 
field  find  the  most  probable  area  and  its  probable  error. 

5674  ±  12,    5680  ±  4,    5685  ±  3,    5682  ±  1,     5678  ±  2 

A    =    5681.41   ±  0.84 

86.  From  the  following  measurements  by  Fizeau   and 
others,  find  the  most  probable  value  for  the  velocity  of  light 
together  with  its  probable  error.     Measurements  are  in  kilo- 
meters. 

298000  ±  1000         299990  ±     200         299930  ±  100 
298500  ±  1000         300100  ±  1000 

V  =    299917  ±  88 

87.  Two  different  instruments  give  for  the  value  of  an 
angle, 

34°  55'  33".0  ±  4".l,  34°  55'  36".0  ±  6".3 

What  is  the  best  value  to  take  for  the  angle  ? 

34°  55'  33".9  ±  3".4 

88.  Determinations  of  the  difference  in  longitude  between 
Washington   and  Key  West  made  on  seven  different   days 
gave  the  results 

I9m     14.42  ±  Os.044  19m     ls.60  ±  0'.046 

1  .37  ±  0  .037  1  .55  ±  0  .045 

1  .38  ±  0  .036  1  .57  ±  0  .047 

1  .45  ±  0  .036 

What  is  the  best  value  and  its  probable  error  ? 

r.4GO  ±  0'.016 
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89.  In  the  triangulation  of  Lake  Ontario  two  different 
instruments  gave  for  an  angle,  74°  25'  5".429  ±  0".29   from 
sixteen  readings,  and    74°  25'  4".611  ±  0".22    from  twenty- 
four  readings.     Find  the  most  probable  value  of  the  angle 
and  its  probable  error. 

90.  In  each  of  Examples  39-45  find  the  mean  and  prob- 
able errors  and  average  deviation  of  each  observation  and  of 
the  most  probable  value,  using  formulas  from  (50)  to  (63) 
according  as  they  apply. 

91.  In  Example  42  divide  the  observations  in  their  order 
into  six  groups  of  four  observations  each  and  compute  the 
mean  of  each  group.     Then  determine  the  probable  error  of 
the  first  of  these  means:  (1)  considered  as  a  single  measure 
of  four  times  the  weight  of  those  in  Example  42  ;  (2)  directly 
as  one  of  six  observations  of  equal  weight;    (3)  as  a  deter- 
mination from  its  four  constituents.  0".67  ;  0".72  ;  1".00 

92.  The  following  twenty-nine  measurements  on  the  den- 
sity of  the  earth,  made  by  Cavendish,  give  as  a  mean  result 
5.48.     What  is  the  probable  error  of  an  observation  ?     Solve 
by  the  usual  method  and  also  by  taking   the  residual  that 
occupies  the  middle  position.  0.14 
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93.  What  is  the  probable  error  of  the  mean  of  two  obser- 
vations which  differ  by  the  amount   a  ? 

94.  A  base-line  is  measured  five  times  with  a  steel  tape 
reading  to  hundredths  of  a  foot,  and  five  times  with  a  chain 
reading  to  tenths  of  a  foot,  with  results 

By  tape,       741.17       741.09       741.22       741.12       741.10 
By  chain,     741.2         741.4         741.0         741.3         741.1 
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Find  the  probable  errors  and  weights  for  a  single  observa- 
tion in  each  cusr,  and  also  the  adjusted  length  of  the  lino 
and  its  probable  error.  741.146  ±  0.015 

95.  Twenty-one  determinations  of  a  chronometer  correc- 
tion gave  results 

-  8.78  -  8.78  -  8.68  -  8.80  -  8.96  -  8.83  -  8.79 
.76  .51  .63  .75  .64  .70  .90 
.85  .64  .58  .78  .65  .64  .93 

Find  the  probable  error  of  the  mean  by  using  both  formula* 
(53)  and  (57),  and  also  determine  the  probable  error  of  ;i 
single  observation  by  taking  the  middle  residual. 

0.017;    0.018;    0.09 

96.  In  the  following  observations  show  that  MQ  =  49.64, 
H  —  1.95,    r  =  1.31,    /*o  =  0.40,    r0  =  0.27,    ^  =  0.87, 
rz  =  0.59. 

M   —    48.81       48.76       49.53       51.56       50.38       49.84 
/>    =        5  4  5  3  2  5 

97.  Observations  on  the  time  of  ending  of  a  transit  of 
Mercury  are  made  by  different  observers  with  a  variety  of 
instruments  and  under  more  or  less  favorable  circumstances. 
If  the  weights  assigned  by  the  computer  are  as  indicated,  find 
the  best  value  for  the  time  and  its  probable  error. 

5*  38m  23*     wt.  1  38m  26*     wt.  3  38m  19*     wt.  3 

37  55        "    0  38    21        «    2  38    21        «   2 

38  10        "    1  38    18        "    2  38    15        «   2 

tQ  =   5*  38m  19*.9 

98.  An  angle  is  measured  five  times  with  a  theodolite,  and 
seven  times  with  a  transit,  giving  results 

Theodolite,    31".7,   39".8,   40".7,   28".6,   3'J  .8 

Transit,          32  .8,    36  .7,    38  .1,    -J!)  .8,    41  .6   35".3,    36".2 
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If  the  relative  values  of  readings  by  the  two  instruments 
are  as  3  to  2,  what  is  the  most  probable  value  of  the  angle  ? 
What  is  the  mean  error  of  the  result? 

99.  Given    Jtfi    =    65.58   ±   .59,    M*    =    35.15   ±   .93, 
MI   =   49.64  ±  .27,    find  the    probable  errors  of   4  J/L  — 

3  Mt  +  2  J/3   and  of   ^  +  ^  -  ^  3.69 ;   0.43 

2*  o  4 

100.  The  three  angles  of  a  triangle  are  measured,  and  the 
probable  error  of  each  observation  is    r  .     What  is  the  prob- 
able error  of  the  triangle  error  ?  r  J~%~ 

101.  The  zenith  distance  of  a  star  on  the  meridian  is 
observed  to  be    z    =    21°  17'  20'  .3  ±  2".3.     The  declina- 
tion of   the  star  is  given  as     d   =    19°  30'  14".8  ±  0".8. 
What  is  the  latitude  of  the  place  and  its  probable  error  ? 

L    =    z  +  d    =    40°  47'  35".l  ±  2"4. 

102.  The  zenith  distance   z   of  a  star  at  upper  culmina- 
tion   is   observed    n    times,  and   its   zenith    distance    z1    at 
lower  culmination    n'   times.     If   the  latitude    is   given   by 
L   =   90°      -   |  (z  +  z'),    and   the   probable    error   of   an 
observation  is   r,   what  is  the  probable  error  of  the  latitude  ? 

103.  The  horizontal  force  necessary  to  start  a  100-pound 
weight  sliding  along  a  table  is  observed  to  be  15.5  ±  0.2 
pounds.     Find  the  probable  error  of  the  coefficient  of  friction. 

104.  If  a  line  is  measured  by  the  continued  application  of 
a  unit  of  measure,  and   r   is  the  probable  error  of  the  placing 
and  reading  of  this  measure,  what  is  the  probable  error  of 
the  length   /  ?  r  y/T 

105.  If  the  average  deviations  of   Zj,  z2,  z8,   are  a,  £,  c, 
respectively,  what  is  the  average  deviation  of  z-f  -\-  z22  -\-  z8a  ? 

106.  If  the  radius  of  a  circle   is  measured   with  result 
1000.0  ±  2.0,   how  should   the  circumference  and  area   be 
expressed  ? 

107.  Two  sides,   a   and   b,   and  the  included  angle   C  of 
a  triangle   are   measured  with    results    a  =   252.52  ±  .06 
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feet,    b   =   300.01   ±   .06   feet,    C   =    42°  13'  00"  ±  30". 
What  is  the  urea  and  its  probable  error?  25452  ±  9 

108.  Measurements  of  adjacent  sides  of  a  rectangle  gave 
a  ±  ri,    and    b  ±  rz.      What  is  the  probable  error   of  the 
an  :i,  and    for  what    kind  of  a  rectangle  will    this    probable 
error  be  tlie  least? 

109.  If  the  measured  sides  of  a  rectangle  have  the  same 
'/.'/.,     what  is  the    a.d.    of    the    diagonal    determined    from 
them  ?  Same 

110.  If   the   sides   of   a   rectangle  are  measured  in   the 
manner  indicated  in  Example  104  and  found  to  be  a   and   b, 
what  is  the  probable  error  of  the  area  ? 

111.  The   correction  to  be  applied  to  a  chronometer  is 
found  to  be   +  12m  13*.2  ±  0*.3.     Ten  days  later  the  cor- 
rection is  again  determined  and  found  to  be    12m  21*.4  ±  0*.3. 
What  is  the  mean  daily  rate  and  its  probable  error  ? 

0'.820  ±  0'.042 

112.  Measurements  of   the    compression    of    the   earth's 
meridian  have  resulted  in 

±  .000046 

294 

What  is  the  probable  error  of  the  denominator  294  ?        3.98 

113.  The   current   flowing    in    a   circuit   is   due    to    two 
sources  whose  electromotive   forces  are  determined   to   be 
ej_    =    200  ±2,    e2    =    400  ±  3.     The  resistance  of  the 
circuit  is   30  ±  1.     Find  the  current  and  its  probable  error. 

20  ±  0.68 

114.  The  side    b   and  angles   B   and    C  of  a  triangle  are 
measured   with    results     b    =    106   ±   .06   metres,     B    = 
29°   39'   ±   1',     C    =     120°  7'   ±  2'.      What   is   the   most 
probable  value  of  the  angle   A   and  of  the  side    c  ? 

A    =    30°  14'  ±  2'.2;    c    =    185.35  ±  0.15 

115.  The  distance  between  two  divisions  on  a  graduated 
is  measured  by  a  micrometer.     Show  that  the  average 
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deviation  of  the  mean  of  two  results  is  the  same  as  the  aver- 
age deviation  of  a  single  reading. 

116.  If  the  weights  of  the  determinations  of  three  angles 
A,  B,  (7,   are   3,  3,  1,   respectively,  what  is  the  weight  of 
the  sum  of  the  three  angles  ?  0.6 

117.  If  the  weight   of   x   is  />,    what  is  the  weight   of 


118.  If   x    =  -  -   and  the  weight  of   y   is  />,   what  is  the 
weight  of   x  ?  c*p 

119.  In    Example    107,  how   closely  must  the  parts   be 
measured  in  order  to  obtain  the  area  within  0.5  per  cent? 

120.  From  observations   on    I  and   t  the  value  of   g   is 
to  be  computed  by  the  pendulum  formula 


What  changes  in  g  will  be  produced  by  changes  in  I  and 
t  of  §!  and  82  units,  respectively,  and  what  are  the  allow- 
able errors  in  I  and  t  if  g  is  to  be  determined  within 
1  per  cent  ? 

121.  The  moment  of  inertia  of  a  cylindrical  bar  is  to  be 
obtained  from  measurements  on  its  mass  m,  its  length  A, 
and  its  diameter  d.  The  error  in  the  determination  of  m 
is  negligible,  the  precision  of  the  determination  of  d  is  four 
times  that  of  h.  If  the  measurements  give  m  =  48, 
A  =  8.000,  d  =  1.200  ±  0.10,  and 


12          16 


what  is  the  probable  error  of   7J   and  what  should  be  the 
ratio  of   d  to    A    to  determine    I  most  accurately  ': 

<!  :  A   '        -j.-irt  :  9 
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122.  If   observations  give  for  a  certain  quantity    x   the 
value  303,   with  a  mean  error  of   2,    what  is  tin-  mean  error 
of  the  expression    3x  -\-  Iog10  2  x  ? 

123.  The    probable    error    of   the    determination    of   the 
angle    A   is  20".     What  is  the  maximum  probable  error  of 
sin  A  4~  cos  A  ? 

124.  If  the  probable  error  of  an  observation  on  an  angle 
is    10",   is  there  any  difference  between  the  probable  error 
of  the  function    sin  A  -|-  cos  A  +  sin  C  and  of  the  func- 
tion   sin  A  -\-  cos  7?  -f~    sin  6Y,    supposing   A    and    J?   are 
of  the  same  magnitude  ? 

125.  Given  the  observations, 

z1-2z2+    z3  -  3    =    0        3z!  4-    z24-2z8—  17    =    0 
2Zl  4-  3z2  -  4z8  -  2    =    0      -  25!  -f-  4z2  +  3z8  -  10    =    0 

Find  the  most  probable  values  of   z1?  z2,  z3,    and  also  their 
weights  and  precision  measures. 

2l    =    3.541;    pzt    =    29;    rtl    =    .024 

126.  Find  the   weights   and    precision   measures   of   the 
unknowns  in  Examples  48  to  57. 

127.  Determine  the  probable  errors  of  the  constants  in 
Examples  72  to  79,  inclusive. 

128.  The  length  of  a  pendulum  which  beats  seconds  is 
given  by 

I    =    I'  -f  /  —  q  --   s\  I'  sin2Z 

where    I'    is  the  length  at  the  equator,    q    the  ratio   of 

289 

the  centrifugal  force  at  the  equator  to  the  weight,  and    s   the 
compression  of  the  meridian  regarded  as  unknown.     Putting 

I'    =    991—  +  x,         (f  ?  -   «)  *'    =    V> 
observations  in  different  latitudes  gave  in    millimetres    the 
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following  equations,  from  which  we  are  to  determine    /   and 
8   together  with  their  probable  errors  :  — 


x  +  0.969y  =  5.13 
x  +  0.749y  =  3.97 
x  +  0.426y  =  2.24 
x  +  0.095y  =  0.56 

x         =  0.19 

x  +  0.152y  =  0.77 
x  +  0.327y  =  1.70 
x  +  0.685y  =  3.62 
x  -f  0.793y  =  4.23 

/'    =    991.069  ±  .026;       0    =    —  ±  0.00046 

294 

129.  Find  the  weights  and    precision   measures    of   the 
unknowns  in  Examples  64  to  70,  and  also  in  Examples  59  to 
63,  and  in  71. 

130.  In  Example  95  the  probable  error  of  a  single  obser- 
vation is  0.08  seconds.     Find  the  number  of  errors  which 
should  fall  between  0.00  seconds  and  0.10  seconds,  between 
0.10  seconds  and  0.20  seconds,  and  also  the   number   that 
should  be  over  0.20  seconds.     Compare  the  results  with  the 
number  actually  found. 

131.  In   470   determinations  of   the  right  ascensions   of 
Sirius  and  Altair  made  by  Bradley,  the  probable  error  of  a 
single  observation  was  0".2637.    The  number  of  errors  falling 
between  specified  limits  was  as  shown  below.     Compare  this 
result  with  the  distribution  of  errors  called  for  by  the  theory. 

Limits.            Errors.  Limits.  Error*. 

0".0  to  0".l  94  0".6  to  0".7  26 

0  .1  to  0  .2  88  0  .7  to  0  .8  14 

0  .2  to  0  .3  78  0  .8  to  0  .9  10 

0  .3  to  0  .4  58  0  .9  to  1  .0  7 

0  .4  to  0  .5  51  Over  1  .0  8 

0  .5  to  0  .6  36 

132.  What  is  the  probability  that  the  error  of  a  single 
observation   will  be  as  large  as  twice  the   probable  error? 
AH  large  as  five  times  the  probable  error  ? 
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133.  On  the  average  how  many  observations  must  In-  made 
before  an  error  as  large  as  three  times  the  mean   error  will 
occur  ? 

134.  In  Example  46,  assuming  that  all  errors  between  any 
two  limits  fall  half  way  between   those  limits,   compute  the 
average  deviation  and  mean  error  of  an  observation  and  com- 
pare their  ratio  with  the  theoretical  value  given  in  the  table 
in  paragraph  55. 

135.  A  line  is  measured  500  times  and  the  probable  error 
of  each  observation   is  0.6  c  m.     Plow  many  errors  should 
occur  between  0.4  c.m.  and  0.8  c.m.  ? 

136.  Show  how  the  value    of    TT    could   be   determined 
experimentally  from  observations  such  as  those  in  Example 
131. 

137.  In  a  system  of  observations  all  equally  good,  r  being 
the  probable  error  of  a  single  observation,  if  two  observations 
are  taken  at  random,  what  quantity  is  their  difference  as 
likely  as  not  to  exceed,  and  what  is  the  probability  that  the 
difference  will  be  less  than   r?  r^T\    0.367 

138.  In  the  following  measurements  of  an  angle,  ought 
any  of  the  observations  to  be  rejected  ? 

12'  51".75         47".85         47".40         48".90         44".45 
48  .45         51  .05         48  .85         50  .95 
50  .60         47  .75         49  .20         50  .55 

139.  Determine    whether    any    of    the    observations    in 
Example  44  should  be  rejected. 

140.  A  quantity   M  is  measured  with  the  results  given 
below.     Ought  all  the  observations  to  be  retained  ? 

M  =   236,   251,   249,   252,   248,  254,   246,   257,   243,   '274 

141.  A   certain  angle  has  been  laid  out  with  such  accu- 
racy that  its  true  value  may  be  taken  as  exactly  90°.    Twenty- 
live  observations  are  made  upon  it  with  a  transit  that  it  is 
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desired  to  test,  and  the  result  obtained  is  89°  59'  57"  ±  0".8. 
What  are  the  odds  in  favor  of  a  constant  error  in  the  instru- 
ment between  —  1"  and  —  5"?  Between  0"  and  —  6"? 

908  :  92  ;     86  :  1 

142.  Repeated  measurements  of  a  standard  metre  bar 
with  a  decimetre  scale  gave  a  result  10.032  ±  0.010.  What 
are  the  odds  in  favor  of  a  constant  error  in  the  scale  between 

0.010         ,  0.054  0 

—  and  -      -  ?  43  :  7 

v/To 


143.  Two   determinations   of  the  length  of  a  line  gave 
683.4  ±  0.3  and  684.9  ±  0.3,   respectively.     Show  that  the 
best  value  for  the  length  is    684.15  ±  0.51,    and  that   the 
probable  systematic  error  of  each  determination  is   0.65. 

144.  Two  men   A   and   J3   observe  an  angle  repeatedly 
with  the  same  instrument  with  results 

A.  B. 

47°  23'  40"       23'  35"  47°  23'  30"       24'  00" 

23    45         23  40  23    40        23    20 

23    30  23    50 

Is  there  any  relative  personal  error,  and  what  is  the  best 
final  value?  47°  23'  38".2  ±  1".6 

145.  Three  independent   determinations  of  the  capacity 
of  a  condenser  made  with   three  different  instruments  gave 
results    42.22    ±    0.21,    43.40    ±    .15    and    44.20    ±    0.18. 
What  is  the  most  probable  value  of  the  capacity? 

For  extended  treatment  of  the  subject  illustrated  in  Exam- 
ples 143  to  145  see  Johnson,  "The  Theory  of  Errors  and 
Method  of  Least  Squares,"  chap.  vii. 

146.  In  an  estimation  of  tenths  what  is  the  probable  error 
of  an  observation  ?     What  is  the  average  deviation  ?         0.025 

147.  In  obtaining  the  angle  of  deflection  of  the  needle  of 
a  tangent  galvanometer  by  the  usual  method  what  is  the 
probable  error  of  the  result  ? 
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148.  If  all  the  errors  of  a  series  of  observations  must  fall 
Itctuccn    0    and    a,    ami  tlu-  fri-queiicy  of  any  error  is  pro- 
portional to   its  magnitude,   what  is  the  Curve  of   Error? 

What  are  the  values  of    r,    a.<t.,    and    /x  ?  r    =    • 

1 49.  In  Example  148  what  is  the  probability  that  the  error 
of  a  single  observation  will  be  as  large  as    0.5  a. 

150.  If  all  values  of   x  between    0    and   a   are  possible, 
and  their  probabilities  are  proportional  to  their  squares,  find 
the  mean  value  of   x   and  the  probability  that   x   will  be  as 
large  as   0.5  a.     Also  draw  the  Curve  of  Error. 

151.  What  is  the  greatest  probable  error  of  a  logarithm 
found  by  interpolation  in  a  seven  place  table?       .000000015 

152.  Given  the  following  set  of  Normal  Equations,  together 
with    [mm]    =    1.3409,    find  the  most  probable  values  of 
the  unknowns  and  their  weights  and  probable  errors.     There 
were  sixteen  observations. 

3.1217  z,  4-     .5756  z2  -     .1565  z3  -     .0067  z4 

-   1.5710      =    0 
.5756  zl  +  2.9375  z2  +     .1103  z3  -       .0015  z4 

+     .9275      =    0 

-  .1565  zl  4-  .1103  z2  4-  4.1273  z3  +  -2051  z4 

+  .0652   =  0 

-  .0067  zx  -       .0015  z2  4-     .2051  z3  -f  4.1328  z4 

4-     .0178      =    0 
zl    =       0.583  ±  0.018 
z4    =  _  0.004  ±  0.015 
pz<   =    4.12 

153.  From  ten  observation  equations,  for  which  was  found 
[mm]    =    2.6322,    there  resulted  the  normal  equations 

5.2485  zl   -   1.7472  z2  -   2.1954  zs  +  0.5399    =    0 

-   1.7472  zx  4-  1.8859  z2  -f  0.8041  z8  - .    1.4493    =    0 

-   2.1954  zl  +  0.8041  z2  +  4.0440  z3  -   1.8681 
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Find  the  most  probable  values  of   z1?  z2  and  z8   together  with 
their  probable  errors.  zt    =    0.42  ±  0.11 

154.  Find  the  most  probable  values  of  the  unknowns  in 
the  normal  equations 

459  z,  _  308  z2  -  389  z8  +  244  z4  _   507  =  0 

-  308  Zi  +  464  za  +  408  z3  —  269  z4  +  695  =  0 

-  389  zl  +  408  z2  +  676  zs  --  331  z4  -f  653  =  0 
244  Z!  -  269  z2  -  331  z8  +  469  z4  -  -  283  =  0 

[mm]    =    1129 
z4    =    -  0.488 ;  pZt    =    281 

155.  If  thirteen  observation  equations  give  rise  to  the 
result    [mm]    =    100.34   and  to  the  normal  equations 

17.50  zl  -     6.50  z2  -     6.50  z8  -     2.14    =    0 

-  6.50  Zj  +  17.50  za  -     6.50  z8  —  13.96    =    0 

-  6.50  zx  -     6.50  z2  -f  20.50  z8  +     5.40    =    0 

show  that  the  most  probable  values  of  the  unknowns  are 

zl   =   0.67  ±  0.60,     z,   =   1.17  ±  0.60,     z8   =   0.32  ±  0.55 


APPENDIX. 


ELEMENTS   OP  THE    THEORY   OF   PROBABILITY. 

200.    Definition.    If  an  event  can  happen  in    a   ways,  and 
fail  in    b   ways,  and  all  these  ways  are  equally  likely  to  occur, 

the  probability  of  the  happening  of  the  event  is  — ^-, 

and  the  probability  of  its  failure  is   ; — r. 

a  -\-  b 

Since  the  event  must  either  happen  or  fail,  the  sum  of  the 
above  probabilities  must  represent  a  certainty.     But 


That  is,  the  probability  of  a  certainty  is  expressed  by  unity. 
Also,  if  the  probability,  J3,  of  the  happening  of  an  event  is 
known,  the  probability  of  its  failure  is  given  at  once  by 
I-  P. 

201.  Example  A.  A  single  throw  is  made  with  a  pair  of 
dice.  What  is  the  probability  that  the  sum  of  the  spots 
turned  up  will  be  5  ? 

Number  of  ways  of  throwing  the  dice  is      6  X  6  =  36 
Number  of  ways  of  throwing  five  is  4 

4  1 

.•.    Probability  of  throwing  five  is  = 

36  9 

Example  B.  A  coin  is  tossed  up  six  times.  Find  the 
chance  that  three  heads  and  three  tails  will  be  the  result. 
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Number  of  ways  of  throwing  the  coin  is  28    =   64 

6x5x4 

Number  of  ways  of  throwing  three  heads  is   -  -   =   20 

1X2X3 

20  5 

Probability  of  throwing  three  heads  is  —  =  — 

64          16 

202.  Compound  Events.  A  certain  event  can  happen  in 
a  ways  and  fail  in  b  ways  :  a  second  independent  event  can 
happen  in  a'  ways,  and  fail  in  b'  ways,  all  of  these  ways 
being  equally  likely  to  occur. 

To  find  the  probability  of  the  simultaneous  occurrence  of 
the  two  events. 

The  total  number  of  ways  in  which  the  events  can  take 
place  together  is  (a  -f-  b)  (a1  -\-  b') 

(1)  Both  events  can  happen  in   a  a'   ways. 

(2)  Both  events  can  fail  in   b  b'   ways. 

(3)  First  event  can  happen  and  second  fail  in    a  b'    ways. 

(4)  First  event  can  fail  and  second  happen  in    a'  b   ways. 

The  probability  of    (1)    is 


(a  +  b)  (a'  +  b') 

h  h' 

The  probability  of    (2)    is 


The  probability  of    (3)    is 


(a  +  b)  (a1  +  b') 
ab1 


(a  +  b)   (a'  +  b') 
The  probability  of    (4)    is    (g  +  ^  (*g,  +  ^ 

But  the  probability  of  the  happening  of  the  first  event  is 

; — r,   and  of  the  second  event  is  —. — ; — r.,   etc.     Hence  it 

a  +  b  a '  -\-   b" 

will  at  once  be  sei-n  that  the  probability  of  the  simultaneous 
occurrence  of  two  independent  events  is  equal  to  the  product 
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of  the  probabilities  of  the  occurrence  of  the  component  events. 
Or,  in  LC«'iicr;il,  if    /J,,  7\,  .  .  .   /'M    arc  the    probabilities  of 

the  occurrence  of  any  number,   />,   of  Independent  event*,  the 

probability  of  tin-  simultaneous  occurrence  of  all  the  events  is 

P,    X    A    X     •  •  •    Pn  (A) 

By  independent  events  is  meant  those  such  that  the  manner 
of  occurrence  of  one  has  no  influence  upon  the  manner  of 
occurrence  of  the  others. 

203.    Example   C.     The  chance  that   A   can  solve  a  cer- 

2 

tain  problem  is    -,    and  the  chance  that    B    can  solve  it  is 
3 

Find, 
12 

(a)  The  probability  that  both  will  solve  it. 

(b)  The  probability  that  the  problem  will  be  solved. 

For  (a).  This  is  a  question  as  to  the  probability  of  the 
concurrence  of  two  independent  events.  Therefore  by  an 
application  of  (A),  the  probability  that  both  will  solve  the 
problem  is 

?     A       A 

3  '      12  18 

For  (b).     The  problem  will  be  solved  unless  both  fail. 

17  7 

The  probability  that  both  will  fail  is        -  y  —   —   — 

8  ;     12          86 

7  29 

The  probability  of  getting  a  solution  is    1 —   — 

36  36 

Example  D.     A  pack  of  cards  is  cut,  and  those  taken  off 
then  replaced.     In  how  many  trials  will  it  be  an  even  u 
that  an  ace  will  be  cut  ? 
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Let    n    be  the  number  of  trials.     Then    n    is  to  be  found 
from 

i-(™}"=  1 

^52  /  2 

where  the  first  member  of  the  equation  represents  the  proba- 
bility that  we  shall  not  fail   n   times  in  succession. 
Solving  for   n, 


n 


log  52   —   log  48 

=    8.7 

In  nine  trials  then  there  is  a  little  more  than  an  even  chance 
of  cutting  an  ace. 

204.  Dependent  Events.    If  we  have  a  number  of  events 
whose  modes  of  occurrence  are  dependent  one  upon  another, 
the  probability  of  their  concurrence  will  be  found  by  the  same 
method  as  in  paragraph  202  ;    a'   now  denoting  the  number 
of  ways  in  which  after  the  first  event  has  happened   the 
second  will  follow,  and    b'   the  number  of  ways  in  which  after 
the  first  has  happened  the  second  will  not  follow,  etc.    Accord- 
ingly, the  general  formula    (A)    of  paragraph  202  applies  to 
dependent  events  as  well  as  to  independent  ones. 

Also,  if  an  event  can  take  place  in  a  variety  of  ways,  the 
total  probability  of  its  occurrence  will  be  the  sum  of  the 
probabilities  of  its  occurrence  in  each  of  the  different  ways. 

205.  Example  E.     Suppose  two  purses  contain  respect- 
ively five  dimes  and  a  copper,  and  six  dimes.     A  coin  is  taken 
at  random  from  the  first  purse  and  placed  in  the  second,  and 
then  a  coin  is  transferred  from  the  second  to  the  first.     What 
is  the  probability  that  the  copper  will  remain  in  the  first  purse  ? 

The  probability  that  the  copper  will  be  taken  from  the  first 
purse  and  placed  in  the  second,  and  then  returned  to  the  first 
purse  is 

I     I       _L 

6   :  49 


and  the  probability  that  the  copper  will  not  be  taken  from 
tlu-  first  purse  at  all  is 


6 

Therefore  the  probability  that  the  copper  will  finally  remain 
in  the  first  purse  is 

J_         _5^  3C  6^ 

42  "  "   6~  42  7 


FUNCTIONS   OF   SEVERAL    VARIABLES. 

206.  For  the  application  of  Taylor's  Theorem  to  the 
expansion  of  a  function  of  several  independent  variables,  see 
Osbor  tie's  "  Differential  and  Integral  Calculus,"  page  145. 
And  for  the  conditions  that  lead  to  maxima  and  minima 
values  of  such  functions,  see  page  155  of  the  same  work. 
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TABLE  II. —  Common  Logarithms. 
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1.563      1.588      1.613      1.638      1.664 
1.823     1.850      1.877      1.904      1.932 
2.103     2.132     2.161      2.190     2.220 

22 
24 
2tf 
28 
30 

1.5 
.1.6 
1.7 
1.8 
1.9 

2.250     2.280     2.310     2.341      2.372 
2.560     2.592     2.624     2.657      2.690 
2.890     2.924     2.958     2.993     3.028 
3.240     3.276     3.312     3.349     3.386 
3.610      3.648     3.686     3.725     3.764 

2.403     2.434     2.465     2.496     2.528 
2.723     2.756     2.789      2.822     2.856 
3.063     3.098     3.133     3.168     3.^U4 
3.423     3.460     3.497      3.534     3.572 
3.803     3.842     3.881      3.920     3.960 

32 
31 
36 
38 
40 

2.0 
2.1 
2.2 
2.3 
2.4 

4.000      4.040     4.080     4.121      4.1fi2 
4.410     4.452     4.494     4.537     4.580 
4.840     4.881      4.928     4.973     5.018 
5.290     5.336     5.382     5.429     5.476 
5.760      5.808     6.856      5.905     5.954 

4.203     4.244      4.285     4.326     4.368 
4.623      4.666      4.709      4.752      4.796 
5.063     5.108     5.153     5.198     5.244 
5.523     5.570     5.617     6.664     5.712 
6.003     6.052     6.101      6.150     6.200 

42 
44 
46 

48 
60 

2.5 
2.6 
2.7 
2.8 
2.9 

6.250     6.300     6.350     6.401      6.452 
6.760     6.812     6.864     6.917     6.970 
7.2!K)     7.344      7.398      7.453     7.508 
7.840     7.896      7.952      8.009     8.066 
8.410     8.468      8.526     8.585     8.644 

6.503     6.554     6.605     6.656     6.708 
7.023      7.076      7.129      7.182      7.230 
7.563     7.618     7.673     7.728      7.784 
8.123     8.180     8.237      8.294     8.352 
8.703     8.762     8.821      8.880     8.940 

52 
54 
66 
58 
00 

3.0 
3.1 
32 
3.3 
3.4 

9.000     9.060     9.120     9.181      9.242 
9.610      9.672     9734      9.797     9.860 
10.24      10.30      10.37      10.43      10.50 
10.89      1096      11.02      11.09      11.16 
11.56     11.63     11.70     11.76     11.83 

9.303     9.364     9.425     9.486      9.548 
9.923     9.986      10.05      10.11      10.18 
10.56      10.63      10.69      10.76      10.82 
11.22      11.29      11.36      11.42      11.49 
11.90      11.97      12.04      12.11      12.18 

62 
6 

7 

7 

7 

3.5 
3.6 
3.7 
3.8 
3.9 

12.25      12.32      12.39      12.46      12.53 
12.96      13.03      13.10      13.18      13.25 
1369      13.76      13.84      13.91      13.99 
14.44      14.52     14.59      14.67      14.75 
15.21      15.29      15.37      15.44      15.62 

12.60      12.67      12.74      12.82      12.89 
13.32      13.40      13.47      13.54      13.62 
14.06      14.14      14.21      14.29      14.36 
14.82     14.90      14.98      15.06      15.13 
15.60      16.68      15.76      15.84      15.92 

7 
7 
8 
8 
8 

4.0 
4.1 
4.2 
43 
4.4 

16.00     16.08      16.16      16.24      16.32 
16.81      16.89      16.97      17.06      17.14 
17.64      17.72      17.81      17.89      17.98 
18.49      18.58      18.66      18.75      18.84 
19.36      19.45      19.54      19.62      ly.71 

16.40      16.48      16.56      16.65      16.73 
17.22      17.31      17.39      17.47      17.66 
18.06      18.15      18.23      18.32      18.40 
18.92      19.01      19.10      19.18      19.27 
19.80     19.89      19.98      20.07     20.16 

4.5 
4.6 
4.7 
4.8 
4.9 

20.25      20.34     20.43     20.52     20.61 
21.16      21.25      21.34      21.44      21.53 
22.09      22.18      22.28      22.37      22.47 
23.04     23.14      23.23      23.33      23-13 
24.01      24.11      24.21      24.30     24.40 

20.70     20.79     20.88     20.98     21.07 
'Jl.t.2      21.72      21.81      21.90      22.00 
22.56     22.06     22.75     22.86     22.94 
23.52      li3.G2      23.72      23.81      23.91 
24.60      24.60      24.70      24.80      24.  IM) 

10 
10 
10 

5.0 
5.1 
5.2 
5.3 
5.4 

25.00     25.10     25.20     25.30      25.40 
26.01      26.11      26.21      2«;.;;j 
27.04      27.14      27.25      27.35      L'T.Jt, 
28.09     28.20     28.30     2841      L'K.W 
29.16      29.27      29.38      29.48      2H.5U 

25J50     25.60     25.70     25.81      25.91 
•;<      26.73      26.83      26.94 
27.56      27.67      L'7.77      27.88     27.98 
-i     K.M     ir.Mtf 

29.70      29.81      29.92      30.03      30.14 

10 
10 
11 
11 
11 

n 

01234 
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TABLE  III.— Squares  of  Numbers. 


T\jff 

5.5 
6.6 
5.7 
58 
6.9 

30.25     30.36     30.47     30.58      30.69 
31.70      31.81 
32.4!»      :<•_'.<;<)      32.72      32X5      32.K5 
•;'      34.11 
HM      :.4.!«3     35.05     35.16     35.28 

30.80     30.91      31.02      31.14     31.25 

B      32.:> 

34.22      M                        34.57     34.69 
35.40     36.52                         li     35.88 

11 
11 
12 

12 
12 

6.0 
(i.l 
SJ 

0.4 

36.00     36.12      36.24     36.36     36.48 

17.21     :;:.33    37.45    37.58    37.70 

88.41      38.56      38.69      38.81 
3i».ti9      3S»>2      :w.«J4      40.07      40.20 
40.96     41.0U      41.22     41.34     41.47 

36.60     36.72     36.84     36.97      37.09 
37.82     37.95      38.07     38.19     38.32 
39.06     «i».m     ::n.:;i     ;'/.'.  n 
40.32     40.45     40.58     40.70     40.83 
41.60      41.73      41.86      41.99      42.12 

12 
12 
13 
13 

13 

6.5 

6.7 

ti.8 

6.9 

42.25     42.38     42.51      42.64     42.77 

•i:;.:>r,     4::.i::i     4:i.*2     43.'.»t>     44.09 
44.89      45.02     45.16      4f,.29     45.43 
4<i.L'4      4r..3S      4I5.M      4H.65      4«.71» 
47.01      47.75      47.89     48.02     48.16 

42.90     43.03     43.16     43.30     43.43 
44.22      44.36      44  4'.l      44.62      44.7»i 
45.56     45.70     45.83     45.97      40.  Id 
46.92      47.06      47.20      47.  .S3      47.47 
48.30     48.44     48.58     48.72     48.86 

13 
13 
14 
14 
14 

7.0 
7.1 

7.2 
7.3 
7.4 

49.00     49.14     49.28     49.42     49.56 
50.41      50.55     50.69     50.84     50.98 
1.86      .V_'.i:{      52.27      52.42 
53.29     53.44     53.58      53.73     53.88 
54.76      54.91      55.06     55.20     55.35 

49.70     49.84     49.98      50.13     50.27 
51.12     51.27      51.41      51.55     51.70 
52.56      52.71      52.85     53.00     63.14 
54.02     54.17      54.32     54.46     54.61 
55.60     55.65     65.80     55.95     50.10 

14 
14 
15 
15 
15 

7.5 

7.6 

7.7 
7.8 
7.U 

56.25     56.40     56.55     56.70     56.85 
67.76     57.91      58.06      68.22     58.37 
59.29     59.44     59.60     59.75     59.91 
60.84     61.00     61.15     61.31      61.47 
62.41      62.57     62.73     62.88     63.04 

57.00     57.15     57.30     57.46     57.61 
58.52      58.68     58.83     58.9* 
60.06     60.22     60.37     60.53      t.ii.r.x 
61.62     61.78     t.l.'.-l      i-j.oy     62.25 
03.20      63.36     63.52     63.68     63.84 

15 
15 
16 
16 
16 

8.0 
8.1 
8.2 
8.3 
8.4 

64.00     64.16     64.32     64.48      64.64 
C5.61      65.77     65.93     66.10     66.26 
67.24     67.40     67.57      67.73     67.90 
68.89      09.06      69.22      69.39      69.56 
70.56      70.73     70.90      71.06     71.23 

64.80     64.96     65.12     65.29     65.45 
66.42     66.59     6675     66.91      67.08 
68.06     68.23     68.39     68.56     68.72 
69.72     69.89     70.06     70.22      70.39 
71.40     71.67      71.74      71.91      72.08 

16 
M 
17 
17 
17 

8.5 
8.6 
8.7 
88 
8.9 

72.25     72.42     72.59     72.76      72.93 
73.96      74.13     74.30     74.48      74.C5 
75.69      75.86      76.04     76.21      76.39 
77.44      77.62     77.79      77.97      78.15 
79.21      79.39     79.57      79.74      79.92 

73.10      73.27      73.44     73.62     73.79 
74.82     75.00     75.17      75.34     75.52 
7H.r.»;      7C.74      7t;.!'l      77.09      77.26 
78.32      78.50     78.68     78.85     79.03 
80.10     80.28     80.46     80.64     80.82 

17 

17 
18 

18 

18  . 

9.0 
9.1 
9.2 
9.3 
94 

81.00     81.18     81.36     81.54      81.72 
82.81      82.99     83.17     83.36      83.54 
84.64     8482      85.01      85.19     85.38 
86.49     86.68     86.86     87.05     87.24 
88.36     88.55     88.74      88.92     89.11 

81.90     82.08     82.26     82.45     82.63 
83.72     83.91      84.09      M.21 
85.56      85.75     85.93     86.12      86.30 
87.42     87.61      87.80     87.98      88.17 
89.30     89.49     89.68     89.87     90.06 

18 
II 

li 

in 

19 

o.r, 

9.6 
9.7 
9.8 
9.9 

80.25     90.44     90.63     90.82     91.01 

'.i-j.lt;      ;ij.;;r,      <c'.:.4      92.74      'X'.itt 
94.09      94.28      94.48      94.67      94.87 
'.it;.  i>4      !H>.24     96.43     96.63     96.83 
98.01      98.21      98.41      98.60      98.80 

91.20     91.39     91.58     91.78     91.97 
'.•::.  1-'                                 93.70     93.90 
95.06      U5.-.«;      n:,.4.r»     95.65     96.84 
•J7.H2      S«7.1.'2      H7.42      1'7.»  ! 
99.00     99.20     99.40     99.60     99.80 

19 
19 
20 
20 
20 

n 

01              234 
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